UNCLASSIFIED 


AD  NUMBER 

AD870195 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  NOV  1969.  Other  requests  shall 
be  referred  to  Commanding  Officer,  Naval 
Training  Device  Center,  Attn:  Code  423, 
Orlando,  FL  32813. 

AUTHORITY 

USNTDC  ltr,  9  Dec  1971 


THIS  PAGE  IS  UNCLASSIFIED 


Jto 

Oi 

rH 

O 

t- 

00 

Q 


Technical  Report:  NAVTRADEVCEN  68-C-0050-1 

ADVANCED  SUBMARINE  SYSTEMS  EQUATIONS  STUDY 

(Final) 


Brenton  R.  Groves,  Ph.  D. 
James  T.  Dorsey 
Dennis  Tucker 


Goodvear  Aerospace  Corporation 
Akron,  Ohio 

Contract  N61339-68-C-0050 


November  1969 


DoD  Distribution  Statement 

This  document  is  subject  to  special  export 
controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be 
made  only  with  prior  approval  of  the 
Commanding  Officer,  Naval  Training  Device 
Center,  (428),  Orlando,  Florida  32813 


D  C 

IrMSEETOlff 

IUJ  JUN  16  1270 


Accessiany 

CfSTI 

GDC 

UKANMMCEB 

JUSItfiGJJiOK 


WHITE  SEttlO* 
(UKF  SECtlOX 


lEttlO*  □  S' 
iECTlOK  IJ  ' 

D 


aiSiS'miM/MiusiUTt  wars 
3i£T.  AuR.  -T  SfECiAl 


Technical  Report?  NAVTRADEVCEN  68-C-0050-1 
ADVANCED  SUBMARINE  SYSTEMS  EQUATIONS  STUDY 


ABSTRACT 


- fhe^goal  of  this  study  was  to  determine  the  optimum  sets  of  equations  to 

be  used  with  two  general  types  of  -submarine  trainers.  The  equations  of 
motion  used  were  those  developed  by  the  Naval  Ship  Research  anti  Develop¬ 
ment  Center.  They  were  reduced  into  two  sets  for  use  in  training  simu¬ 
lators,  one  for  a  basic  submerged  control  maneuvering  trainer  and  one 
for  an  advanced  submerged  control  simulator  to  provide  highly  realistic 
ship  control  training  through  the  full  range  of  normal  and  casualty 
conditions. 


The  report  outlines  a  general-purpose  digital  computer  program,  fallowing 
the  N5RDC  standard  equations,  written  in  FORTRAN.  Integration  methods  for 
digital  simulators  are  discussed.  A  number  of  programs  for  testing  the 
degree  of  simulation  of  a  digital  simulator  program  are  given. 


Use  of  these  programs  as  applied  to  submarine  simulation  is  shown  with 
two  sets  of  equations  which  eliminate  74  of  the  131  coefficients  used 
in  the  original  NSRDC  equations.  The  mathematical  model  for  a  submarine 
simulator  using  a  very  small  computer  is  developed.  A  discussion  of 
near-surface  operation  and  wave-generation  is  followed  by  the  general 
requirements  for  determination  of  the.  simulation  requirements  for  training 
as  opposed  to  research. 
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FOREWORD 

The  NSRDC  (Naval  Ship  Research  and  Development  Gaiter)  developed  a  set  of 
equations  which  describes  submarine  motion  in  six  degrees  of  freedom.  The 
Naval  Training  Device  Center  recognized  this  set  of  equations  as  the  "standard" 
set  which  should  be  utilized  in  submarine  control  trainers  and  on  all  trainer 
procurements  required  that  the  complete  set  be  utilized.  A  medium  scale  digital 
computer  is  required  to  implement  the  entire  set  of  NSRDC  equations  in  a  submarine 
trainer.  Since  the  coefficients  of  many  of  the  equations  terms  are  zero  or  have 
little  effect  on  submarine  trajectories,  the  possibility  of  reducing  the  equations 
to  a  size  which  permits  small  scale  digital  computer  implementation  in  submarine 
trainers  was  recognized  by  NAVTRADEVCEN. 

This  report,  along  with  its  coup  anion  reports,  NAVTRADEVCEN  68-C-0050-2  and 
NAVTRADEVCEN  68-0-00^0-3,  is  the  culmination  of  a  study  by  Goodyear  Aerospace 
corporation  to  determine  abbreviated,  or  optimum,  equations  to  be  used  in 
conjunction  with  submarine  simulators. 

This  report  should  be  an  invaluable  aid  to  simulator  personnel  during  development 
of  submarine  simulators.  It  presents  a  discussion  of  several  factors  involved  in 
simulation,  such  as,  integration  methods,  wave  generation,  small  computer  utiliza¬ 
tion  and  verification.  It  not  only  presents  an  abbreviated  set  of  equations  wherein 
a  significant  number  of  the  original  coefficients  are  eliminated,  but  presents  a 
tool,  in  the  form  of  computer  programs  which  permit  further  coefficient  reduction 
investigations.  Since  the  required  complexity  of  equations  are  directly  pro¬ 
portional  to  the  training  device  requirements,  the  programs  are  written  so  that 
the  investigator  may  taylor  the  coefficients  to  reach  any  desired  level  of  simula¬ 
tion. 

NAVTRADEVCEN  63-C-00£0-2  presents  descriptions,  flow  charts  and  listings  of  each 
computer  program  and  NAVTRADEVCEN  68-C-0050-3  presents  results  of  the  computer 
programs  using  the  SS(N)5>9U  submarine  as  the  demonstration  model. 
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SECTION  I 
INTRODUCTION 


There  are  two  general  types  of  submarine  .trainers:  (1)  The  basic 
high  speed  submarine  trainer j  and  (2)  Tne  advanced  submerged  control 
trainer.  The  first  type,  which  does  not  necessarily  duplicate  a  specific 
class  of  submarine,  is  used  to  provide  initial  submarine  control  indoctri¬ 
nation  to  potential  pianesraen  and  diving  officers.  The  second  type  provides 
highly-realiatic  ship  control  training,  which  extends  over  the  operating 
ranga  of  the  submarine  to  ship  control  personnel.  Both  types  of  trainers 
utilize  computers,  programed  with  equations  of  motion,  to  simulate  sub¬ 
marines  under  various  operating  conditions. 

In  the  paet,  the  equations  used  in  these  programs  were  somewhat 
arbitrarily  chosen  by  manufacturers  of  training  devices.  Then  an  attempt 
was  made  by  the  Naval  Ship  Research  and  Development  Center1-*  toward 
standardizing  the  equations. 

This  brings  us  to  the  purpose  of  this  study* 

1.  To  determine  an  optimum  set  of  equatio/is  of  motion  to  be  used 
with  each  of  the  two  general  types  of  submarine  trainers. 

2,  To  write  the  two  sets  of  equations  as  compactly  as  possible  in 
order  that  they  may  be  used  with  a  small,  general-purpose  computer. 

Two  additional  reports  resulted  from  this  contract  and  are  described 
as  follows: 

Report  Number  Remarks 

68-C-00S>0*  3  Report  is  CONFIDENTIAL  and  titled,  "Advanced 

Submarine  Systems  Equations  Test  Data  (U)B. 

68-C-0050-2  Ret  ,  .is  UNCLASSIFIED  and  titled,  "Advanced 

Submarine  Systems  Programming''1. 

♦Superscript  numbers  indicate  references. 
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SECTION  II 

STATEMENT  OF  THE  PROBLEM 


The  determination  of  optimum  sets  of  equations  to  be  used  with  the 
two  general  types  of  submarine  control  trainers  forms  the  problem.  The 
tasks  required  to  solve  this  problem  are  the  determination  of  the  simu¬ 
lation  requirements.  Application  of  these  requirements  in  an  analysis 
of  the  equations  of  motion,  programing  the  equations  of  motion  on  a 
digital  computer,  and  then  verification  that  any  changes  to  optimize  them 
do  not  exceed  the  limits  specified  for  good  simulation. 

A.  PROGRAMING  THE  GENERAL  NSRDC  EQUATIONS 


The  problem  analysis  is  based  on  the  assumption  that  the  NSRDC 
equations  are  accurate  for  the  simulation  of  operational  submarine  motion. 
Therefore,  the  first  task  is  to  program  these  equations  so  that  a  digital 
program  representing  the  complete  equations  of  motion  is  available.  The 
most  practical  programing  is  achieved  with  the  use  of  a  problem-oriented 
language  (FORTRAN)  so  that  the  programs  are  independent  of  the  computer 
used.  This  will  allow  the  output  of  the  study  to  be  used  by  any  researcher 
interested  in  submarine  simulation  on-  any  computer  equipped  with  a  FORTRAN 
compiler. 

Because  this  program  is  going  to  be  a  research  tool,  it  is  necessary 
to  build  in  the  greatest  degree  of  flexibility  possible.  Therefore, 
consideration  must  be  given  to  allowing  any  changes  in  parameters,  variables 
and  analytical  methods  within  the  framework  of  the  basic  equations  without 
having  to  recompile  subroutines  at  each  run. 

The  final  result  should  be  a  program  that  will  input  coefficients 
easily,-  produce  any  outputs  desired,  and  will  make  a  permanent  record 
of  any  changes.  Also,  it  should  be  written  so  that  many  runs  can  be  made 
without  having  to  start  the  entire  program  over  again. 

The  next  problem  is  to  verify  the  research  program  against  the  avail¬ 
able  data  on  submarine  performance  to  determine  the  accuracy  of  the  program. 
A  complete  set  of  runs,  following  established  submarine  test  practice, 
are  run  using  the  NSRDC  coefficients,  and  these  are  compared  with  the 
NSRDC  reports  on  a  specific  demonstration  submarine.  These  runs  verify 
the  program  and  establish  proper  rerponse  in  various  input  conditions  and 
also  establish  a  data  base  with  which  the  reduced  programs  can  be  compared. 
In  effect,  a  test  guide  is  completed,  and  any  approximations  made  to  the 
program  must  oe  justified  against  the  base  program. 

B.  MODIFIED  AND  BASIC  EQUATIONS 

After  the  research  program  is  operational,  it  will  be  necessary  to 
establish  the  requirements  for  the  basic  and  modified  programs.  The 
first  step  will  be  to  analyse  the  existing  performance  data  on  submarines 
and  the  use  of  this  data  in  present  submarine  simulators*  The  main  sources 
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of  data  will  be  the  NSRDC  reports  and  the  dynamic  reports  on  present  day 
trainers.  Also,  a  number  of  studies  made  in  the  past  will  have  'to  be 
reviewed, 

A  reduction  to  a  basic  equation  set  can  be  made  after  the  simulation 
limits  are  set.  The  coefficients  can  be  divided  into  broad  categories 
and  the  effect  of  eliminating  certain  ones  can  be  studied.  Computer 
runs  are  used  to  aid  the  simplifications  of  the  coefficient  terms.  For 
example,  if  it  is  desired  to  know  the  effect  of  a  cross- coupling  term, 
runs  will  be  performed  with  and  without  the  coupling  term  in  question* 

The  importance  of  the  term  can  easily  be  established  when  the  two  runs 
are  compared. 

The  above  procedure  can  indicate  the  effect  of  any  change^  but  it 
does  not  indicate  the  effect  of  response  on  the  trainee  when  he  is  part 
of  the  closed-loop  controlling  the  simulator.  This  is  one  of  the  hardest 
areas  to  cover  in  a  study  of  this  type  because  closed-loop  real- time 
operation  cannot  be  run  on  a  typical  batch  process  digital  conqouter  oper¬ 
ation.  An  autopilot  will  be  programed  as  part  of  the  longitudinal  loop 
to  keep  the  submarine  on  a  level  course  when  turns  axe  made*  Measurement 
of  the  longitudinal  parameters  over  time  as  the  autopilot  corrects  for 
effects  of  lateral  displacement  can  be  used  to  determine  the  nature  and 
quality  of  dynamic  responses.  This  procedure  is  far  more  effective  than 
inserting  step  inputs  and  measuring  the  output  after  a  time  delay. 

After  a  complete  analysis,  the  two  sets  of  equations  will  cover  a 
basic  submerged  control  trainer  and  a  complete  simulation  of  a  particular 
submarine.  These  equations  will  serve  as  both  a  demonstration  of  data 
application  to  a  particular-  simulator  and  as  a  reference  for  any  submarine 
trainer  or  simulator  in  the  future. 

C.  ANALYSIS  OF  SURFACE  CONDITIONS 

The  coefficients  used  for  the  analysis  are  deep-surface  coefficients 
and  are  not  applicable  to  a  submarine  moving  close  to  the  surface.  The 
change  in  the  coefficients  due  to  the  presence  of  a  free  surface  can  be 
calculated  by  hydrodynamic  methods,  and  such  a  set  of  coefficients  has 
been  generated  by  NSRDC.  This  set  is  used  in  the  research  program  with 
values  adjusted  for  different  depths  and  the  result ing  given  frequencies 
compared  with  the  deep  submergence  set  to  determine  the  surface  effect 
on  the  submarine.  An  analysis  will  indicate  the  importance  of  any  change. 

The  biggest  disturbing  effect  of  a  suomarine  operating  close  to  the 
surface  is  that  of  waves  breaking  overhead.  Near-surface  operation  is 
becoming  more  important,  so  a  good  wave  simulation  is  essential  for  a 
modern  submarine  simulator. 

Wave  action  can  bs  described  only  on  a  statistical  basis  because 
the  sea  surface  is  never  regular,  Tt  is  a  jumble  of  different  waves 
of  different  heights  all  changing  and  all  moving  in  different  directions. 

This  action  can  be  described  by  means  of  spectrum  analysis  and  the  problem 
to  be  solved  i3  how  to  apply  this  technique  to  a  pratical  submarine  simulator 
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The  height  of  waves  above  the  submarine  is  the  first  thing  that 
has  to  be  known.  This  is  usually  specified  by  a  one-dimensional  spectrum, 
particularly  that  of  Neumann  which  was  presented  in  1953.  At  least  three 
other  spectral  densities  have  been  used  since  then,  so  the  first  task  of 
this  part  of  the  problem  was  to  determine  which  of  the  various  spectral 
densities  hove  the  most  suitable  form  for  both  submarine  response  and  for 
ease  of  computation  with  a  small  digital  computer.  Enough  background  will 
be  given  to  show  how  any  spectrum  can  be  used  in  a  wave  generator,  and 
one  will  be  chosen  and  programed  to  illustrate  the  method  on  a  practical 
basis. 

The  addition  of  a  direction  function  to  the  wave  spectrum  is  neces¬ 
sary  because  the  difference  between  direction  of  ship  travel  and  the  direc¬ 
tion  of  wave  propagation  has  a  direct  bearing  on  the  forces  felt  by  a 
submerged  submarine.  Therefore,  part  of  the  solution  will  be  to  develop 
this  directional  function  and  show  how  it  is  applied.  The  calculated 
forces  and  moments  will  be  determined  to  give  a  realistic  motion  under 
all  conditions  of  wave  inputs.  A  very  large  number  of  papers  and  ref¬ 
erenced  have  been  published  in  this  area  in  the  past  few  years.  These 
will  be  surveyed  and  a  list  drawn  up  that  will  enable  future  workers  to 
bring  themselves  up  to  date  quickly  on  the  present  state-of-the-art  in 
wave  generation. 

D.  SMALL  COMPUTER  ANALYSIS 

The  last  problem  to  be  investigated  once  the  program  limits  have 
been  determined  is  that  of  a  small  computer  analysis. 

Wien  the  modified  and  basic  program  outlines  have  been  set  and  a 
«»ximutn  reduction  has  been  made,  it  is  possible  to  determine  the  neces¬ 
sary  core  size  to  hold  each  of  the  programs .  The  analysis  should  try 
to  reduce  this  factor  as  much  as  possible  because  it  is  directly  propor¬ 
tional  to  the  cost  of  the  computer  being  used. 

The  programs  are  then  analyzed  from  the  total  number  of  tasks  to 
be  performed.  This  includes  not  only  the  arithmetic  operations  but  also 
all  input- output  subroutines  and  any  other  bookkeeping  operations  neces¬ 
sary  to  make  up  the  complete  cycle.  All  operations  are  sorted  by  category 
and  time  to  determine  the  necessary  computer  speed  for  real-time  operation. 
Finally,  the  reduced  programs  will  be  run  on  a  16-bit  computer  and  compared 
against  the  programs  run  on  the  IBM  360  32-bit  computer. 

A  small  computer  surrey  and  state-of-the-art  review  will  be  made  as 
of  the  date  of  this  study.  This  review  will  be  out  of  date  soon  at  the 
rate  at  which  the  computer  industry  is  advancing,  but  it  will  serve  to 
indicate  an  idea  of  what  can  be  done  in  future  submarine  simulators. 
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SECTION  in 
METHODS  AND  PROCEDURES 

A.  PROGRAMING  THE  NSRDC  EQUATIONS 


The  NSRDC  standard  equations  of  motion  for  submarine  simulation 
cover  all  phases  of  submarine  simulation  in  six-degrees-of -freedom, 
including  emergency  recoveries  after  casualties.  These  equations  come 
from  NSRDC  Report  2510  which  contains  a  brief  history,  defines  the  mathe¬ 
matical  model,  discusses  the  coefficients  required,  and  sets  a  standard 
to  be  used  in  the  simulation  of  submarines.  These  equations  were  pro¬ 
gramed  in  two  forms.  They  are  given  in  program  EB920  -  Submarine  Simu¬ 
lation,  and  ZC7?0  -  Submarine  Simulation,  Longitudinal  Freedom  Only 
found  in  the  programing  report,  NAVTRADEVCEN  68-C-0050-3. 

EB?20  is  a  research  submarine  simulation  program  that  allows  the 
simulation  of  any  submarine  on  a  large  general-purpose  digital  computer. 

It  is  complete  and  flexible  in  that  all  NSRDC  coefficients  are  programed, 
and  changes  can  be  made  easily  between  runs. 

Figure  1  is  a  block  diagram  of  the  program  with  all  subroutines 
shown.  The  program  operates  as  follows. 

The  main  submarine  simulation  program  ties  the  various  subroutines 
together  and  determines  the  operating  limits  placed  on  each  run.  The 
first  call  is  to  INPUT.  This  subroutine  reads  all  coefficients  from 
a  data  deck  along  with  controls  to  set  the  output,  determines  which 
integration  subroutine  to  use,  and  which  control  input  schedule  should 
be  used.  The  main  program  then  writes  out  all  coefficients  so  that  the 
input  to  the  run  has  a  permanent  record  printed  at  the  start.  An  integra¬ 
tion  switch  then  calls  either  the  KUTTA  or  INTEG  subroutine.  KUTTA  is  a 
very  accurate  method  of  integration  that  is  used  as  a  standard  for  compar¬ 
ing  changes  between  runs.  INTEG  has  three  short  methods  of  integration 
built  in  and  can  be  easily  adapted  so  that  any  integration  method  can  be 
used  without  changing  the  rest  of  the  program. 

Depending  on  which  integration  subroutine  is  used,  EVIL  or  EVAL  1 
is  called.  They  both  compute  the  equations  of  motion  as  given  in  equations 
(1)  through  (7)  from  NSRDC  Report  2510  .  EVAL  has  matrix  invert  and 
multiply  subroutines  so  that  the  equations  can  be  solved  exactly  in  each 
cycle.  EVAL  1  uses  values  computed  on  the  last  cycle  in  some  cases. 

Both  subroutines  call  CONTR  which  contains  schedules  for  a  number  of 
inputs  representing  the  deflection  of  the  controls  on  the  actual  submarine. 
A  return  to  INTEG  integrates  all  parameters  and  the  cycle  is  finished. 

The  main  program  then  prints  out  the  twelve  variables  of  interest  on 
demand  and  checks  the  RUN  TIME  to  see  if  the  program  is  finished. 

Plotting  calls  have  been  added  to  the  main  program  so  that  a  plot 
can  be  made  of  an?  of  the  twelve  output  variables  and  surface  control 
deflections  versus  time.  These  plots  can  be  made  in  any  order  and  the 
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Figure  1.  Research  Simulation  Program  Block  Diagram 
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cor vp.  t  labeling  will  be  added  to  the  plot.  However,  program  runs  do 
noi  -'•epend  on  the  plotting  subroutines,  so  if  a  plotter  is  not  avail¬ 
able,  the  printed  output  serves  in  place  of  the  plots. 

Almost  all  variables,  in  addition  to  the  coefficients,  are  located 
in  the  COMMON  area  of  the  program.  This  allows  a  degree  of  flexibility 
that  could  not  be  obtained  otherwise  and  the  main  program  takes  advantage 
of  this  to  maxe  changes  between  runs.  After  a  run  is  finished,  all  values 
are  restored  to  their  original  values  at  the  start  of  the  run.  The  next 
cards  in  the  input  deck  are  then  checked  and  a  coding  system  allows  any 
variable  in  COMMON  to  be  changed  to  a  new  value.  Integration  methods 
and  intervals,  input  schedules,  and  output  requirements  as  well  as  the 
coefficients  can  be  changed  between  each  run.  Because  of  the  heavy  use 
of  COMMON,  details  are  given  in  the  Programing  Report. 

ZC790  is  identical  to  EB920  except  that  only  three- dsgrees-of -freedom 
are  simulated.  Only  the  longitudinal  mode  is  active  so  all  terms  multi¬ 
plied  by  p, v ,  or  r  are  set  to  zero  and  removed.  The  block  diagram  is 
shown  In  figure  3.  Only  one  method  of  integration  is  supplied.  The 
CONTR  subroutine  has  removed  all  provisions  for  rudder  motion.  Therefore, 
only  meander,  overshoot,  and  similar  runs  can  be  calculated..  This  program 
was  written  to  save  computer  time  because  it  runs  in  one-third  the  time  of 
EB920.  This  is  a  real  advantage  where  a  large  number  of  longitudinal 
runs  are  required.  About  three-fourths  of  the  total  runs  made  during 
the  course  of  this  study  were  of  this  type.  The  mathematical  model  for 
program  ZC790  is  given  in  equations  (8)  through  (11). 

The  mathematical  model  for  the  simulation  program  used  two  different 
sets  of  notation  because  of  the  FORTRAN  restriction  that  required  the 
use  of  capital  letters  only.  Therefore,  the  following  notation  is  adapted 
and  carried  throughout  this  study.  The  letters  X,  Y,  Z,  K,  M,  and  N 
stand  for  coefficients  in  the  axial,  etc.,  channels  respectively.  U,  V, 

W,  P,  Q,  and  R  stand  for  the  velocities,  and  whenever  an  absolute  value 
is  stated,  the  velocity  will  be  proceeded  by  an  A.  Accelerations  are 
velocities  followed  by  a  D,  The  surfaces  are  proceeded  by  a  D  for  delta 
and  eta  terms  are  followed  by  an  E.  Inertias  and  other  terns  are  similarly 
identified.  For  example, 


»  XRP 

Ylrlt-r  * 

z  vV  *  ZWD 

V,  Msb  96  mdb 

N  ivivhrv  **  navave 
IXY*  1X1 


B.  INTEGRATION  METHODS 


This  section  discusses  soma  of  the  many  integration  methods  that 
can  be  used  for  digital  computation,  why  they  were  chosen  to  implement 
the  equations  of  motion  for  submarines,  and  how  they  are  utilised  in 
the  submarine  simulation.  The  optimum  integration  step-size  and  a  stabil¬ 
ity  criterion  for  any  integration  method  is  discussed  also. 
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Figure  3.  Research  Simulation  Program,  Longitudinal  Freedom 

Block  Diagram 
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Thera  are  a  great  number  of  integration  methods  that  can  be  used 
in  digital  computation.  Some  of  these  methods  have  been  in  use  for  many 
years.  They  were  developed  to  integrate  functions  numerically  that  could 
not  be  solved  in  closed  form.  Others  have  been  developed  fairly  recently 
to  take  advantage  of  the  digital  computer’s  high  speed  and  ability  to  handle 
large  bloeke  of  numbers.  They  range  all  the  way  from  the  simple  Euler’s 
method  to  a  fourteen-step  predictor-corrector  developed  by  NASA  to  calculate 
the  position  of  the  planets  with  great  accuracy. 

The  choice  of  any  integration  method  is  determined  by  considering  the 
trade-off  between  accuracy,  stability,  and  speed.  The  accuracy  is  depend¬ 
ent  on  the  amount  of  truncation  and  round-off  errors  present  when  using  a 
particular  method.  Truncation  error  results  from  the  fact  that  digital 
integration  is  a  discrete  process  using  a  finite  number  of  equally  spaced 
points  rather  than  a  series  with  an  infinite  number  of  terms.  Accuracy 
can  be  imp  roved  by  using  more  points  as  a  past  history  or  by  making  the 
integration  step  interval  smaller.  Either  method  buys  increased  accuracy 
at  the  cost  of  speed.  Round-off  error  results  whui  generating  a  sequence 
of  numbers  to  represent  the  integrated  function.  Each  number  is  limited 
by  the  size  of  the  computer  word  being  used.  However,  there  are  methods 
for  determining  an  upper  bound  on  the  magnitude  of  round-off  errors  ex¬ 
pected  for  the  numerical  accuracy  being  used.  Studies  have  shown  that 
for  five  decimal  figure  accuracy  within  the  limits  of  a  sixteen-bit 
maohine,  round-off  errors  are  of  no  oonsequenoe  compared  to  the  expected 
truncation  errors, z 

Stability  of  an  integration  method  is  defined  as  the  ability  to 
control  propagated  errors  so  that  they  do  not  increase  as  timo  passes. 
Propagated  errors  are  the  sum  of  both  truncation  and  round-off  errors 
over  all  integration  steps  calculated’ since  the  problem  started.  The 
amount  of  propagated  error  is  the  single  most  important  criteria  in  the 
choice  of  any  integration  method  for  simulation.  Fortunately,  there  is 
a  comparatively  easy  procedure  for  determining  if  a  particular  integration 
method  can  be  used  with  a  particular  simulation  problem.  This  procedure 
will  be  explained  below  using  the  submarine  simulation  as  an  example. 

Speed  depends  on  the  complexity  of  the  method  used  and  the  size  of 
the  integration  interval.  Real-time  simulation  requires  fitting  all 
computations  for  one  cycle  within  the  integration  interval  while  still 
keeping  the  desired  accuracy.  Therefore,  this  trade-off  las  to  be  made 
for  each  simulation.  In  the  case  of  this  study  it  was  proven  that  the 
problem  was  not  critical  for  submarines. 

Integration  methods  can  be  divided  into  CLOSED  and  OPEN  classes. 

If  a  rate  at  some  specified  tins  is  required  to  calculate  a  position, 
the  scheme  is  said  to  be  a  CLOSED  (or  corrector)  formula.  If  the  rate 
is  not  required,  the  scheme  is  an  OPEN  (or  predictor)  one.  Closed 
methods  cannot  be  used  to  integrate  the  equations  of  motion  because  the 
rates  are  a  function  of  the  position  at  that  time.  Since  the  positions 
are  the  output  of  the  calculations  integrated  from  the  rates,  the  rates 
are  not  known.  On  the  other  hand,  open  formulas  use  rates  calculated  in 
prior  integration  cycles  so  that  the  present  rates  are  not  required.  The 
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last  rates  are  used  in  the  open  integration  formulas  to  give  the  present 
position  and  the  new  positions  are  used  to  calculate  the  present  rates 
from  the  equations  of  motion.  In  order  to  keep  the  different  methods 
straight,  a  notation  has  been  developed  in  vhieh  the  open  and  closed  form¬ 
ulas  are  specified  by  Onm  and  Cnm  respectively,  and  n  is  the  number  of 
former  positions  used  and  m  is  the  number  of  former  rates  used.  The 
difference  equations  expressing  the  integrated  position  as  a  function 
of  ratss  and  position  is 

Xu  — J  +  h 

Pi  {*o 

where  i  is  ths  discrete  time  at  vhich  a  calculation  has  been  made  in  the 
past  and  n  ie  the  present  time.  The  a's  and  b*s  are  parameters  vhich 
determine  the  particular  integration  method,  H  is  the  integration  interval 
or  step  size  and  has  tfaa  dimension  of  time  in  order  to  keep  the,  dimensions 
consistent  serose  the  equation.  Note  that  if  b0  zero,  then  Sin  ie  not 
needed  to  calculate  X»  and  the  system  is  an  open  or  predictor  method. 
Therefore,  for  simulation  work,  be  should  be  zero. 


M 

2  »>!**-{ 


For  an  exaa^xLe,  let  us  look  at  the  2nd  order  Adams  method  of  inte¬ 
gration  which  is  an  Oi&  type*  In  this  case,  3|  ■  l,bo  =  0,  b)  *3/2, 
and  ba  *  -i*  The  difference  equations  would  be 


n-i 


Many  different  integration  formulas  have  been  extensively  studied  - 
a»*d  almost  any  text  on  numerical  methods  goes  into  them  in  some  detail. 


The  problem  was  to  determine  vhich  of  these  methods  should  be  used 
for  a  submarine  simulator  and  vhich  one  should  be  used  for  the  research 
program.  Finding  an  integration  method  for  the  research  program  was 
comparatively  easy.  Since  the  program  did  not  have  to  run  in  real-time* , 
the  only  criterion  that  had  to  be  met  was  to  have  the  most  accurate 
program  available. 

The  method  chosen  is  a  four  step  Rungs-Kutta  predictor-corrector. 

It  was  stated  that  corrector  (closed)  methods  could  not  be  used  in  simu¬ 
lation,  but  in  this  case  two  difference  equations  are  used.  An  open 
equation  is  used  to  predict  a  position  and  the  closed  equation  is  used 
to  check  this  calculation.  The  sequence  is  as  follows.  A  new  position 
is  integrated  from  the  prior  known  rates  using  an  open  formula.  This 
new  position  is  used  to  calculate  a  new  rate  from  lie  equations  of  motion. 
The  new  position  and  new  rate  can  now  be  used  in  a  closed  formula  to  inte¬ 
grate  a  new  corrected  position.  The  new  corrected  position  is  used  in 
the  equations  of  motion  to  calculate  a  new  corrected  rate.  This  process 
is  then  repeated.  This  is  enough  to  reduce  the  error  between  two  suc¬ 
cessive  position  calculations  to  a  very  small  amount. 

The  exact  difference  equations  are  not  listed  here  but  are  given 
in  the  programming  import  under  the  KUTTA  subroutine. 

The  next  task  was  to  determine  which  integration  method  should  be 
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used  in  a  submarine  simulator.  It  is  obvious  that  spaed  and  simplicity 
are  desired,  provided  that  accuracy  can  be  maintained.  A  number  of 
studies  have  been  performed  using  different  techniques  for  integration 
in  flight  simulation  and,  based  on  these  reports,  two  different  methods 
were  ohosen  for  examination  against  the  Runge-Kutta  method. 


The  simplest  method  of  all  is  an  On  due  to  Euler  with  the  following 
difference  equation. 

Xu®  Xn-I  +  H  %T\-2 

The  last  calculated  rate  is  multiplied  by  the  integration  step  size 
and  is  then  added  to  the  position  known  up  to  that  tine  to  give  the  new 
position.  This  difference  equation  has  a  bad  reputation.  In  fact,  Jennings  3 
says  that  "it  is  a  naive  method  of  very  limited  accuracy  and  not  recommended 
in  practice*.  However,  it  has  only  one  multiplication,  one  addition,  and 
no  storage  is  required  to  hold  former  values.  This  makes  it  the  shortest, 
and  therefore  the  fastest,  of  ary  method  and  since  one  of  the  objectives 
of  this  study  is  to  find  the  smallest  total  program  this  method  was 
examined. 


In  1965  Goodyear  Aerospace  Corp.  performed  an  i  -house  research  and 
development  study-*,  utilizing  their  simulation  laboratory  digital  computer 
facility,  to  investigate  a  number  of  complex  integration  algorithms  used 
for  digital  computers  in  general  to  see  if  such  complexity  was  needed 
for  flight  simulation.  At  that  time,  it  was  proposed  that  digital  flight 
simulators  use  such  formulas  as  Hammings. Cji or  the  o»  CjjMod  Gurk 
developed  at  the  University  of  Pennsylvania,  The  cot  elusions  of  this 
study  were  that  excellent  results,  well  within  the  tolerance  allowed  in 
the  simulation  of  aircraft  systems,  can  be  obtained  with  much  simpler 
algorithms.  It  was  found  that  the  smallest  errors  occurred  among  the 
simpler  routines  with  the  2nd  order  Adams  difference  equation 

X»  «  X-yi-j  +  jfaX-n-i  “in-*) 

This  requires  an  additional  subtraction  and  a  storage  for  each 
parameter  to  be  Integrated.  Since  this  one  has  been  used  very  success¬ 
fully  at  Goodyear  Aerospace,  and  Euler  integration  has  larger  errors, 

2nd  order  Adams  was  also  implemented  in  the  subroutines. 


The  next  step  was  to  look  at  the  way  integrations  are  handled  in 
the  NSRDC  equations  of  motion.  The  six- de gre e- of-f re  edom  problem 
requires  that  twelve  integrations  be  performed,  because  the  second  order 
equations  have  been  expressed  as  a  system  of  simultaneous  first  order 
equations.  Using  matrix  notation 

[M]R]=  M 

where  [*]  are  the  rates  and  M  are  the  coefficient  summations.  The 
matrix  £Jw|l  comes  from  the  fact  that  the  rates  are  functions  of  other, 
rates  and,  therefore  have  to  be  solved  /or  simultaneously.  If  t^JsLS 
'W.Vf-r  ,  then  this  matrix  is 
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O  o  o  niZq  -TnYg 

O  Tti-f^Y*  o  -f^rp-inZc,  O 

o  o  m-lfe  It) Yq  -|^ZvHlXs  o 

o  -f£V^ws  ix-f^p  -i*Y  -rXz-|jerK^ 

O  “Ixy  ~Iy£ 

"'YwYig  ^  “IxYJ^p  'Iyz 

The  reason  thatx>V|i,9/  $  and  ¥  do  not  appear  in  the  matrix  is 
that  they  are  not  dependent  on  Cl  etc.,  and  therefore  all  of  the  cross 
terms  are  zero.  This  reduces  the  matrix  to  six  levels  instead  of  twelve. 
The  solution  to  the  equation  is 

fM]  (*]  =  [X] 

M"  [m]  [*]  =  [M-][X] 

[I]  [*]  =  [M-]  [X] 


[*]  =  [M]-'  [X] 

where  is  the  inversion  of  the  £m3  matrix.  As  long  as  the  coef¬ 

ficients  that  make  up  the  [mJ  matrix  do  not  change  between  each  inte¬ 
gration  cycle,  the  matrix  can  be  inverted  once  when  the  program  is 
originally  written  and  the  inverted  matrix  used  every  time  on  the  rigit 
hand  side  of  the  equation.  Unfortunately,  X<j>  Y«  and  are  elements 
in  this  matrix  and  they  represent  the  amount  of  water  in  the  ballast  tanks 
of  a  submarine.  Therefore,  it  is  necessary  for  an  exact  computation  to 
calculate  the  elements  of  the  matrix  on  each  cycle.  The  matrix  is  then 
inverted  with  the  aid  of  a  subroutine.  Each  of  the  rates  are  calculated 
without  the  terms  that  are  functions  of  the  other  rates.  These  particular 
rates  are  multiplied  by  the  inverted  matrix  to  give  the  true  accelerations 
at  that  integration  cycle  time.  The  open  integration  method  then  integrates 
the  rate  to  produce  the  next  set  of  positions.  The  whole  process  is  then 
repeated. 

The  submarine  research  program  has  a  eontrol  card  that  can  be  used 
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to  romnvfl  the  evaluation  of  the  original  matrix  after  the  first  cycle 
if  none  of  the  elements  are  going  to  change  during  the  course  of  the 
program.  This  is  the  usual  case  if  the  tanks  are  not  going  to  be  changed. 


Matrix  inversion  and  multiplication  take  up  both  room  and  time  in 
the  computer  program  as  anyone  who  has  done  matrix  arithmetic  knows. 
Therefore,  a  programing  technique  was  borrowed  from  aircraft  flight  simu¬ 
lation  to  avoid  this  problem.  From  the  NSRDC  equations t 


a 

❖ 

w 

* 


* 


requires 

requires 

requires 

requires 

requires 

requires 


r  i 
*  > 


* 

c 

e 

* 

v 


W  > 

*  > 
V  > 


X 

♦ 

a 


The  usual  case  is  for  Yq  and  the  cross  inertias  to  be  zero  due 
to  symmetry  of  the  submarine/  If  the  matrix  is  not  used, 

«m  us.  (f  n  vU1  uoe  v-n 

v-n  “se?)H)rfr,yh9n<^  winua.  livjjWy, 

Wyi  1,111  m  Ul,,  will  w««  U.y, 


The  other  terms  have  zero  for  a  coefficient.  The  effect  of  this 
procedure  is  a  choice  of  either  having  the  translational  accelerations 
or  the  rotational  accelerations  slightly  in  error,  depending  on  which 
set  is  calculated  first. 


The  next  parameter  to  be  determined  was  the  integration  step  size. 

This  is  a  critical  value  because  it  determines  the  time  available  in 
real-time  to  compute  the  entire  problem  and  service  all  I/O  busses  before 
the  integration  interval  is  used  up.  If  this  step  size  (the  variable  H 
in  all  computer  programs)  is  too  small,  the  computer  will  not  finish  the 
computations  required  in  real-time .  If  H  is  too  large,  the  accuracy 
suffers  and  finally,  if  H  is  increased  beyond  a  certain  limit,  the  inte¬ 
gration  method  used  will  run  into  stability  problems.  Therefore,  the 
trade-off  is  to  make  H  as  large  as  possible,  consistent  with  both  accuracy 
and  stability. 

The  Runge-Kutta  integration  method  supplied  with  the  submarine 
research  program  has  a  feature  that  allows  the  step  size  to  be  adjusted 
between  each  cycle  dependent  on  the  accuracy  required.  Both  the  upper 
and  lower  bounds  on  H  can  be  specified. 

A  typical  program  was  run  with  the  demonstration  submarine  in 
which  the  rudder  was  deflected  to  3 5  degrees  and  stemplanes  were 
deflected  to  fifteen  degrees,  A  time  interval  of  200  seconds  was  considered 
sufficient  to  determine  the  effect  of  the  variable  step  size,  The  maximum 
values  of  the  twelve  parameters  were  printed  out  and  0.1$  of  this  value 
was  used  as  the  maximum  allowable  error  for  any  step.  The  program  was 
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re-run  with  H  allowed  to  vary  between  0.05  and  2,0  seconds  subject  to  the 
accuracy  requirements.  The  run  shewed  that  the  interval  for  0.1$  accuracy 
varied  from  0,25  to  2  seconds.  Many  steps  were  one  second  or  longer, 
indicating  that  the  use  of  one  second  intervals  is  probably  adequate  for 
most  maneuvers,  but  not  necessarily  good  enough  for  disaster  or  emergency 
conditions.  The  smallest  step  size  of  0,25  seconds  occurred  with  step 
inputs  and  this  indicated  that  this  time  interval  was  sufficient  for  any 
purpose.  Therefore,  the  studies  were  run  with  an  H  of  0.25  seconds  in  all 
cases.  However,  any  fixed  step  size  can  be  chosen  with  any  of  the  inte¬ 
gration  methods  used  in  the  subroutines. 

An  additional  study  was  made  into  the  largest  possible  step  size 
baaed  on  a  stability  criterion  that  assumes  the  local  growth  of  the 
propagated  error  must  be  kept  within  bounds.  The  criterion  assumes  that 
each  integration  method  has  a  range  that  is  a  function  of  the  step  size 
and  the  partial  differential  of  the  function  with  respect  to  itself  that 
will  keep  the  propagated  error  bounded.  The  system  may  not  be  stable, 
even  if  this  criterion  ie  mat,  but  it  will  certainly  be  unstable  if  it 
is  not. 

Assume  that  the  function 

i  =  Xs,  etc.) 


is  to  be  integrated  by  means  of  the  2nd  order  Adams  method.  The  numerical 
quantity 

h 

n  an 

has  to  have  a  value  between  0  and  -1  for  the  integration  to  be  stable. 
Looking  at  for  example, 


\where  U.  =  speed 
\J  JL  -  length 

*  YT)  *  101188 

■p  *  density 


If  the  specified  values  for  vhe  demonstration  submarine  are  used, 
the  value  of  the  above  function  at  high  speed  indicates  that  a  step  size 
of  2.1  seconds  could  be  used.  The  other  eleven  integrations  were  checked 
and  they  give  equally  large  maximum  step  sizes.  The  stability  region  for 
some  common  integration  methods  are  listed  in  table  1.  The  theory  for 
calculating  the  stability  region  for  any  method  is  given  in  Nigro^, 


Hie  submarine  research  program  has  a  control  option  that  allows  a 
program  to  be  run  either  under  the  Runge-Kutta  integration  method,  along 
with  the  matrix  inversion  and  multiply,  or  to  choose  any  one  of  a  number 
of  simpler  methods  in  which  the  matrix  is  not  used.  Three  are  coded  into 
the  subroutine  (IHTEG);  namely  Euler,  2nd  order  Adams,  and  the  best 
method  based  on  stability  alone.  Table  If  in  the  programing  report  has  22 
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TABLE  1.  INTEGRATION  STABILITY  REGIONS 


* 


Method 

Stability  Region 

Euler 

-2,  0 

2n<*  Order  Adams 

-1,  0 

Trapezoidal 

-  ,  0 

°30g31 

-.82,  0 

O33  Mod  Gurk 

-.82,  0 

Best  2n{*  Order 

-li,  0 

different  equations  listed  so  that  the  subroutine  can  be  recoded  to 
utilise  most  of  the  popular  numerical  integration  techniques.  This 
makes  it  possible  to  have  a  whole  library  of  integration  methods  and 
to  use  any  one  of  them  by  means  of  a  control  card  at  the  start  of  a  run. 


A  number  of  the  references  express  various  integration  methods 
as  a  Z- transform  transfer  function.  The  Z- transform  is  an  extension  of 
the  traditional  Laplace  transform  methods  with  the  transform  of  a  sampled 
function  written  in  a  simpler  notation.  The  substitution  Z  ■ 0ST  is  made 
where  T  is  the  sample  step  size  and  corresponds  to  the  integration  interval 
H.  The  Z- transform  is  a  handy  way  of  examining  integration  methods  ana¬ 
lytically  because  only  knowledge  of  Laplace  transforms  is  required  to 
use  them,  but  they  can  be  confusing  to  program  on  a  digital  computer. 
Therefore,  a  method  is  given  to  convert  any  Z- transform  transfer  function 
into  a  difference  equation.  It  can  then  be  programed  directly  in  FORTRAN. 


The  key  to  this  procedure  is  that  Z  or  0sTis  a  shifting  function 
in  the  time  domain.  Therefore,  any  negative  power  of  Z  times  a  parameter 
merely  indicates  that  the  value  of  the  parameter  so  many  cycles  back 
from  the  present  time  should  be  used.  An  example  will  make  this  clear. 

Baxter^,  for  example,  lists  several  integration  methods  I(Z).  For 
the  2nd  order  Adams, 

X  (z)  =  (b)  (3 Z-‘  - Z- *)/  (I -2" ') 

The  first  step  is  to  write  I{Z)  as  a  transfer  function  of  Y(Z)  to 

*(Z). 


~  I  -Z-‘ 
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Cross -multiply  to  clear  the  fraction 

yCz)(i-z-,)  =  't'(z)“5-  (3z-‘  -z-a) 
Y(Z)~  Y(z)Z'f  =  (3 1  (Z)  Z- '  -  Y  (Z)  Z’a  ) 

Y(z;  =Y(z;Z''+-|  C3Y(Z)Z'’-Y<Z)Z'2) 


This  procedure  assumes  that  the  Z-transform  is  a  linear  operator 
which  is  to  be  expected  since  the  Laplace  transform  is  linear.  The 
difference  equation  can  now  be  written  directly  by  making  the  substi¬ 
tution 


Y(Z)=Yv,,  Y(Z)Z-'=Yn.t,  etc. 


which  is  the  same  difference  equation  stated  earlier  for  2nd  order  Adams. 
Any  integration  methods  expressed  as  a  Z-transform  transfer  function  can 
be  used  in  subroutine  IN  TEG  by  coding  a  few  lines  of  FORTRAN  and  recom¬ 
piling  this  subroutine  only. 

A  number  of  different  runs  were  made  which  included  both  step  inputs 
and  realistic  surface  deflections.  From  these  runs  it  is  possible  to 
draw  some  definite  conclusions  as  to  step  size  and.  integration  methods 
to  he  used  in  a  submarine  simulator. 

The  Runge-Kutta  method  is  more  accurate  than  is  necessary  even  for 
a  research  program.  However,  it  serves  as  a  useful  check  on  other  methods 
and,  with  a  step  size  of  0.25  second,  is  the  reference  standard. 

It  is  not  necessary  to  use  matrix  multiplication  to  solve  for  simul¬ 
taneous  rates.  Using  the  rates  calculated  in  the  last  integration  cycle 
does  not  affect  the  accuracy  of  the  rune. 

The  minimum  step  size  for  best  accuracy  is  0.25  second  based  on  the 
variable  step  size  program.  However,  there  is  almost  no  difference  when 
the  program  is  run  at  0.5  second.  Whan  it  is  run  at  1.0  second,  all 
parameters  ream  in  the  same  except  for  p.  If  the  bank  angle  output  is 
examined  closely,  it  will  be  seen  that  the  high  frequency  wiggles  are 
out  of  phase  with  the  standard  and  take  a  little  longer  to  die  out.  At 
1.5  seconds,  p  goes  into  oscillation.  Therefore,  a  step  size  of  0.5 
second  should  be  used  for  a  casualty  trainer  where  high  rates  will  be 
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encountered,  A  step  size  of  1,0  second  is  satisfactory  for  a  maneuvering 
trainer. 

It  was  hoped  that  Euler* s  method  could  be  used  for  the  maneuvering 
trainer  because  it  gave  good  results  st  0.25  second  except  for  p.  However, 
at  1.0  second  most  of  the  parameters  diverged,  thus  justifying  its  bad 
reputation.  Therefore,  2nd  order  Adams  will  be  used  in  the  maneuvering 
trainer  so  that  a  large  step  size  can  be  vised. 

C.  COEFFICIENT  REDUCTION 

1.  LONGITUDINAL  IMPULSE  TESTING 

The  equations  of  motion  can  be  represented  as  a  transfer  function 
that  takes  a  certain  set  of  input  conditions  and  transforms  them  into  a 
certain  set  of  outputs.  The  coefficients  of  the  differential  equations 
determine  the  characteristics  of  this  transfer  function,  and  as  long  as 
the  system  under  consideration  is  linear,  it  is  possible  to  express  the 
response  of  the  system  independently  of  any  set  of  inputs  by  means  of 
the  impulse  function. 

The  impulse  function  is  a  spike  that  occurs  at  time  "t*o"  and  shocks 
a  system  into  oscillation.  The  system  will  oscillate  because  energy  has 
been  added.  These  oscillations  will  correspond  to  the  natural  frequencies 
and  denying  ratios  because  all  inputs  are  zero  immediately  after  *t«o". 

The  time  outputs  can  be  analyzed  to  determine  the  coefficients  necessary 
to  describe  a  linear  differential  equation  for  the  system.  These  will  be 
related  to  the  coefficients  required  in  the  equations  of  motion.  The 
techniques  for  doing  this  are  well  known  in  Laplace  transform  theory 
applied  to  servo  design?. 

A  procedure  somewhat  analogous  to  impulse  testing  is  used  in  the 
determination  of  handling  characteristics  during  the  full  scale  test  of 
the  actual  submarine.  The  test  is  known  as  a  "meander  maneuver"  and  one 
can  obtain  a  direct  measure  of  the  inherent  longitudinal  stability  of 
the  submarine  with  the  following  procedure.  Essentially,  the  meander  test 
is  conducted  by  disturbing  the  submarine  from  its  neutral  pitch  angle  by 
deflecting  the  stemplanes  a  fixed  amount  in  either  rise  or  drive.  After 
a  predetermined  pitch  angle  is  reached,  the  planes  are  returned  to  their 
neutral  position.  The  resulting  pitch  angle  trajectory  provides  both 
visual  and  quantitative  evaluations  of  the  longitudinal  stability  of  the 
submarine. 

This  comes  us  close  as  possible  to  impulse  testing  when  applied 
to  the  real-world  vehicle.  Mathematically,  an  impulse  occurs  in  zero 
time,  but  in  the  real-world  caas  the  surfaces  take  time  to  move  and  the 
controls  have  to  be  held  in  the  deflected  position  so  that  the  system 
is  disturbed  enough  to  give  a  good  reading  of  pitch  angle  change.  How¬ 
ever,  when  the  surface  movement  is  plotted  on  the  same  time  scale  as  the 
meander  maneuver,  it  forms  a  very  narrow  pulse  that  approaches  an  impulse. 
Once  the  pitch  angle  has  been  plotted  for  any  onset  speed,  it  can  be  ana¬ 
lyzed  for  such  numerical  maneuvers  as  damping  ratio,  undampened  natural 
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frequency,  stability  roots,  and  trajectoiy  envelope  decay.  These  are 
given  in  the  evaluation  reports  of  the  handling  qualities  for  each  sub¬ 
marine  class. 

It  was  hoped  that  impulse  testing  of  the  equations  of  motion  could 
serve  two  purposes s  (l)  A  series  of  runs  would  be  made  to  determine 
the  natural  frequencies  and  damping  ratios  at  all  speeds:  and  (2)  these 
figures  could  then  be  compared  with  the  actual  figures  cn  the  submarine 
to  check  the  accuracy  of  the  complete  set  of  coefficients.  Once  this 
had  been  done,  the  figures  could  be  used  to  generate  a  new  set  in  which 
the  equations  are  assumed  to  be  linear  and  therefore  much  simpler  than 
the  original  set.  The  first  purpose  was  successful,  the  second  was  not. 

The  following  procedure  is  used  to  generate  an  inpulse  run.  Com¬ 
puter  program  ZC790,  Submarine  Simulation,  Longitudinal  is  used  to  gen¬ 
erate  a  run  with  the  control  subroutine  set  to  generate  an  input  impulse 
and  the  output  option  set  to  punch  cards  with  the  tine  and  pitch  angle 
recorded  every  two  seconds.  The  output  on  these  cards  is  assumed  to  follow 
the  form 


'(=a)e-st  +ajLea't  (cos fit  +  y>) 

where  y  and  t  are  on  the  cards  and  ,  and  ifj  are  con¬ 

stants  that  are  to  be  determined. 


The  reason  for  choosing  this  form  as  an  equation  comes  from  the 
form  of  the  characteristic  equation  for  the  longitude  channel  when  the 
equations  of  motion  are  linearized.  The  Laplace  transform  of  y  is 


^DfCt)]= 


4-  K|S  -f-  Ko 

(S+  X)  (CS-f  +3X) 


,  and  'ft  are  the  natural  frequencies  of  the  system  and  remain 
the  same  in  either  the  time  or  complex  domain  and  G.,,  Q-j.,  and  V 
are  functions  of  these  and  |^,  K,  ,  and  JC0. 

If  we  consider  the  equations  of  motion  as  a  linearized  set  in  the 
complex  domain,  they  can  be  written  as 

CyCs)]  =  [eCs)][c(s;] 

where  Y(S)  is  the  output  matrix,  E(S)  are  the  equations  of  motion  and 
C(S)  are  the  input  conditions.  An  impulse  is  a  constant  in  the  complex 
domain  so  Y(S)  is  equal  to  E(S).  Therefore,  X(t)  in  the  time  domain  in 
response  to  an  impulse  can  be  determined  by  taking  the  inverse  Laplace 
transform  of  the  equations  of  motion  for  the  particular  input  (the  stem- 
planes  in  this  case),  and  the  particular  output  (pitch  angle).  The 
details  of  setting  up  the  equation  of  motion  for  the  longitudinal  channel 
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are  described  in  the  programing  report  in  program  ECUO  and  will  not  be 
repeated  hers.  From  this  analysis  ths  Cwsrscfc6i\Law.c  souation  is 


where 


A3S3+Aa.  S2  +  A|  S  +A0  =  O 

**  -  m)  ***)(r*i) 

M<j) 

(• pD4  x  v  7 oJ)3  .  pQ$ 

T  -r  Uo Z?)-(x  UoM^Xf-  z£> 

.4,- T«o  M f  M?* w-fn)(B  -Z%) 

-  (^^o^vj)(^-atn  +  "•£“  UoZ 

Ao-  (*§  Mo  Zw)CB  *  Z8) 


where  all  the  symbols  are  defined  in  the  equation  of  motion  report. 
Actually,  computer  program  EdUO  solves  this  equation  in  direct  matrix 
form  and  allows  for  the  variation  in  speed. 

A  compensation  is  made  for  the  non-linear  terms  in  the  longitudinal 
loop  by  modifying  the  basic  coefficients.  For  example, 

2wM  •  w(iwz^v9x  becomes  Zwi^j  *7  -  w 

where  w  is  an  average  value  of  w  used  as  a  constant.  It  is  determined 
by  running  a  number  of  full  set  runs  and  scanning  the  output  for  a  typical 
value  of  w  at  steady  state  conditions.  This  allows  ZW  and  ZWAW  to  be 
combined  Into  a  single  term,  A  similar  process  is  used  with  }  being 
set  to£.  “ 

It  is  obvious  that  the  characteristic  equation  i3  a  cubic  with  at' 
least  one  real  root.  It  can  be  assumed  that  the  other  two  roots  are  a 
complex  pair  from  the  behavior  of  the  submarine  in  a  meander  maneuver. 
Hierefore,  if  the  equations  were  linear,  the  frequencies  could  be  deter¬ 
mined  directly  from  the  coefficients  and  conversely,  the  coefficients 
could  be  calculated  from  the  measured  frequency  response. 

When  the  cards  have  been  punched  with  the  time  and  pitch  angle; 
o<  ,/$  t  and  can  be  determined  by  two  other  programs.  EC310  is  a 
program  that  performs  a  ’'least  squares  fit”  to  any  specified  equation 
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format  from  a  series  of  X-X  points.  Thi3  program  requires  that  fairly 
close  initial  guesses  be  supplied  for  the  parameters  that  are  to  be 
calculated.  Program  EC330  calculates  a  set  of  numerical  values  to  be 
used  as  starting  conditions  for  EC310.  Details  are  left  to  the  programing 
report, 

A  summary  of  the  situation  will  be  helpful  at  this  point. 

1.  The  handling  qualities  report  for  a  particular  submarine  gives 
the  results  of  meander  tests  taken  during  full-scale  trials  and  expressed 
in  terms  of  frequencies  and  delays, 

2.  The  stability  and  control  report  gives  the  same  information 
calculated  from  the  coefficients  as  determined  from  the  model  tests. 

3.  Program  EGlItO  uses  the  NSRDC  coefficients  to  calculate  the  fre¬ 
quencies  and  damping  ratios.  The  values  are  identical  to  the  report  if 
the  coefficients  given  in  the  stability  and  control  reports  are  used  in 
this  program. 

U.  Computer  program  ZC790  can  be  used  to  produce  a  time  output 
which  can  be  analyzed  to  findo<  ,  and  y ,  The  natural  frequency 
and  damping  ratio  can  then  be  calculated. 

It  would  be  convenient  if  all  four  sources  agreed  exactly,  but  they 
do  not.  Therefore,  it  was  deemed  advisable  to  determine  the  validity 
of  the  supplied  coefficients  and  the  programing  of  ZC790,  the  main 
simulation  program. 

The  coefficients  were  taken  as  supplied  by  NSRDC  and  run  in  ZC790 
to  produce  a  full  set  of  meander  curves.  These  curves  were  compared 
with  the  new  data  acquired  on  the  demonstration  submarine  and  stored  at 
NSRDC,  and  they  matched  over  the  entire  speed  range.  Therefore,  we  can 
be  assured  that  th9  longitudinal  coefficients  used  and  simulation  program 
Z(?790  (and  E6920)  will'  accurately  simulate  the  demonstration  submarine,  ^ 
The  iangitudinal  and  lateral  channels  of  ££920  were  checked”out  in  another 
manner  in  order  to  give  further  support  to  this  conclusion.  The  handling 
qualities  report  for  the  demonstration  submarine  gives  a  time  history  of 
a  20-knot  vertical-plane  overshoot  maneuver  and  a  time  history  of  a  typical 
2h-knot  submerged  turn.  The  input  conditions  were  duplicated  as  closely 
as  possible  and  the  output  curves  matched  the  actual  output  data  taken 
during  sea-trials  to  a  much  closer  degree  than  required  for  an  accurate 
simulator. 

We  can  now  establish  a  data  platform  from  the  approved  criteria 
for  any  submarine  simulator.  The  coefficients  as  received  from  NSRDC 
are  inserted  into  ZG790  for  longitudinal  tests  or  EC920  for  six-degrees- 
of-freedom  and  any  specified  set  of  inputs  can  be  run  off.  This  will  not 
only  serve  to  check  the  reduction  mde  as  a  part  of  this  study  but  more 
importantly  will  serve  to  check  the  operation  of  any  simulator  using  the 
standard  set  of  equations  of  motion.  Even  before  a  simulator  is  delivered 
to  the  user,  a  complete  set  of  performance  specifications  can  be  computed. 
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An  investigation  was  made  into  why  the  handling  qualities  report 
and  the  stability  and  control  report  should  have  slightly  different 
frequency  and  damping  ratio  parameters  for  what  is  essentially  the 
same  vehicle.  It  was  concluded  that,  the  processing  of  the  data  was  re¬ 
sponsible  for  this  difference.  The  frequencies  in  the  stability  report 
are  calculated  from  the  coefficients  as  determined  by  model  tests  using 
the  linear  model  described  in  ECliiO.  This  neglects  such  coefficients  as 
Mw  w  which  were  proven  to  be  important  by  the  investigation.  The  fre¬ 
quencies  in  the  handling  qualities  report  are  determined  from  a  time- 
history  on  the  assumption  that  a  single  pair  of  complex  poles  are  present 
to  represent  the  system  rather  than  a  real  pole  and  a  pair  of  complex 
ones.  This  is  a  standard  treatment  in  servo  texts  and  is  used  to  calcu¬ 
late  damping  ratios  in  the  verification  section  of  this  study.  It  can 
lead  to  different  numbers  depending  on  which  points  of  the  time-history 
are  used.  Authorities  are  not  in  agreement  as  to  the  way  to  measure  the 
damping  ratio  of  a  higher  order  function,  and  most  of  them  use  only  the 
single  complex  pair  case  as  an  example. 

The  method  of  calculating  the  frequencies  from  the  linear  coefficients 
is  possible  lead  to  an  idea  for  .calculating  a  simplified  set  of  longitudinal 
coefficients  that  seemed  extremely  attractive  in  theory.  The  full  set  of 
coefficients  are  used  to  run  off  a  set  of  impulse  curves  at  all  speeds. 

These  curves  are  reduced  to  numerical  values  by  the  technique  described, 
and  the  reverse  equations  contained  in  EClitO  are  used  to  generate  a  new 
set  of  coefficients  for  all  speeds.  These  would  be  linear,  but  would 
contain  all  the  variations  described  by  the  complete  set  because  of  this 
change  with  speed.  Unfortunately,  this  attempt  proved  to  be  unsuccessful 
due  to  the  difficulty  of  solving  the  reverse  set  on  a  "least  squares*  fit. 
The  approach  is  outlined  here  because  the  problem  is  in  the  computing 
methods  available  and  not  in  the  theoretical  basis  behind  the  problem. 

The  unknowns  are  the  eight  basic  coefficients: 

^  m  f  M  ^  y  Z”  ••£!  jM  \v  >  M  ^ 

which  are  used  to  form  the  characteristic  equation: 

Ay  S*+  AaS*+A($  +  Ao  =0 

If  we  assume  that  the  variation  of  the  longitudinal  stability  indices 
with  speed  are  known,  (these  are  the  art  ,  ft  and  f  's  from  either  the 
impulse  runs  or  the  stability  and  control  report)  the  characteristic 
equation  takes  the  form 

(s+  irXs  +  «)1'-t-/St)  =  O 

Therefore,  we  can  set  up  a  one-to-one  correspondence  with  the 
stability  indices  and  the  A*s, 
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^Va3  -  z°<  +  ir 
a‘/A3  -  z  tr*  +  cx*  +  /9 

A/43  =  r(**+/3*) 

and  then  solve  for  the  basic  coefficients.  If  tests  are  run  at  more 
than  three  speeds,  more  knovns  than  unknowns  are  generated  and  a  best 
fit  of  runs  at  one-knot  intervals,  for  example,  would  give  us  a  basic 
set  that  could  be  used  to  linearize  its  equation  of  motion, 

A  number  of  different  numerical  methods  were  tried  and  discarded 
because  they  all  broke  down  on  the  fact  that  the  slopes  of  the  change 
of  coefficients  with  respect  to  speed  is  so  shallow  that  not  even  double 
precision  subroutines  were  able  to  prevent  the  problem  from  blowing  up. 

Most  of  the  non-linear  equation  solvers  require  an  initial  guess  and  the 
shallower  the  slope,  the  closer  the  guess  has  to  be  to  the  true  answer 
for  all  eight  coefficients.  Usually,  the  original  values  are  suitable  as 
initial  guesses,  but  in  this  case  they  do  not  satisfy  the  requirements. 

The  mathematical  difficulties  are  covered  in  detail  in  Freeburg8.  Never¬ 
theless,  this  approach  might  prove  useful  in  further  work  in  this  field. 

2.  BASIC  COEFFICIENT  ANALYSIS 

The  equations  of  motion  of  a  body  can  be  developed  in  terms  of  force 
equations  derived  from  d'Alembert's  principle.  This  principle  states 
that  for  any  body  in  motion,  the  applied  forces  must  be  in  equilibrium 
with  the  inertia  reactions  of  the  body.  3y  assuming  that  the  inertia 
reactions  have  the  dimensions  of  force,  the  dynamic  problem  can  be  reduced 
to  a  static  one  in  which  the  only  requirement  is  that  of  “'sum  of  forces  a 
0*.  Once  the  force  equations  are  developed  and  expressed  in  terms  of 
coefficients  times  functions  of  velocities,  .. ,  ies,  etc,,  the  accelerations 
can  be  computed.  Known  velocities  with  the  external  forces  (rudder  angle, 
plane  angles,  and  thrust)  are  used  as  input  conditions.  This  assumes 
that  tho  mass  and  inertias  are  known.  The  accelerations  are  then  inte¬ 
grated  to  give  a  new  set  of  velocities  which  are  then  used  to  recompute 
a  new  set  of  accelerations. 

The  coefficients  can  be  divided  into  three  categories  to  facilitate 
the  investigation  of  their  effects  on  the  acceleration  computation  for 
any  set  of  inputs^.  They  are  longitudinal,  lateral,  and  cross-coupling 
coefficients. 

In  the  general  case,  there  are  six  simultaneous  differential  equations 
to  be  solved  for  the  six-degreea-of-fraedom  for  a  body  that  can  rotate 
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And  translate*  However^  IT  certain  conditions  arc  pieced  on  ihc  inputs , 
the  problem  can  be  greatly  simplified.  It  is  helpful  to  start  the  inte¬ 
gration  with  all  forces  exactly  balanced  by  the  control  surface  deflec¬ 
tions  so  that  the  six  accelerations  are  zero.  This  is  fairly  simple  to 
achieve  in  a  computer  program  by  setting  the  steady-state  inputs  and 
velocities  to  the  correct  precomputed  values.  When  the  program  begins 
to  run,  tbs  accelerations  are  zero  and  the  velocities  do  not  change  until 
a  surface  is  deflected. 

The  first  set  to  be  investigated  are  the  longitudinal  coefficients. 
The  program  starts  with  the  rudder  at  neutral  and  the  lateral  axis  veloc¬ 
ities  at  zero.  Since  the  longitudinal  axis  velocities  do  not  couple  into 
the  lateral  channel,  the  lateral  velocities  remain  at  zero  and  any  coef¬ 
ficients  multiplied  by  these  velocities  can  be  ignored  if  only  the  diving 
planes  and  thrust  are  varied.  The  three- degree -of -freedom  research 
program  not  only  assumes  that  this  condition  holds,  but  it  eliminates 
the  lateral  channel  entirely  along  with  any  coefficients  using  lateral 
velocities  in  the  longitudinal  channel. 

Therefore,  only  those  coefficients  multiplied  by  the  longitudinal 
velocities  are  of  interest  in  the  first  category.  They  are  listed  in 
table  2. 


TABLE  2.  LONGITUDINAL  COEFFICIENTS 


Axial  Force 

Normal  Force 

Pitching  Moment 

XQQ 

ZDQ 

XUD 

ZWD 

M3AQ 

XWQ 

ZQ 

MWD 

xou 

ZAQDS 

MQ 

xww. 

ZWAQ 

MAQDS 

XBSBS 

ZSTR 

MAWQ 

XDBDB 

ZW 

MSTR 

Thrust  coefficients 

ZWAW 

MW 

XWWE 

ZAW 

MWAW 

ZDSDSE 

zww 

MAW 

ZDS 

MWW 

ZD3 

MDS 

ZQE 

MDB 

ZWE 

MQE 

ZWAWS 

MWE 

ZDSE 

MWAWE 

MDSE 

The  second  set  of  coefficients  are  those  that  govern  the  behavior 
of  the  lateral  channel.  The  assumption  is  made  that  an  autopilot  is 
attempting  to  keep  the  submarine  level  so  that  q  is  zero  and  the  value 
of  w  remains  a  small  constant.  There  will  be  some  coupling  through  the 
change  of  forward  speed  as  the  vehicle  begins  to  turn,  but  it  is  assumed 
that  this  is  small  colored  to  the  effect  of  the  rudder.  The  lateral 
coefficients  are  given  in  table  3. 
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TABLE  3.  LATERAL  COEFFICIENTS 


Lateral  Force 

Rolling  Moment 

Yawing  Moment 

YRD 

KPD 

NRD 

YPD 

KRD 

NPD 

YPAP 

KPAP 

NRAR 

YVD 

KP 

NVD 

YR 

KR 

NP 

YP 

KVD 

NR 

YARDR 

KSTR 

NARDR 

YVAR 

KV 

NATO 

T5TR 

KVAV 

NSTR 

Y7 

KDR 

NV 

YVAV 

KSTRE 

NVAV 

mt 

NDR 

IRE 

NRE 

YVE 

NVE 

YVAVE 

NVAVE 

YDRE 

NDHB 

Finally,  there  are  the  cross-coupling  terns  that  couple  one  channel 
into  the  other  and  only  cone  into  prominence  during  a  climbing  or  diving 
turn.  The  cross-coupling  terms  are  given  in  table  k* 

TABLE  lu  CROSS -CCUPLIKJ  TERMS 


Axial 

Lateral 

Normal 

Rolling 

Pitching 

Yawing 

XRR 

YPQ 

2PP 

KQR 

MPP 

NPQ 

XRP 

YQR 

ZRR 

K?Q 

MRR 

NQR 

XTO 

YVQ 

ZRP 

KVQ 

MRP 

NWR 

XW 

YWP 

ZTO 

KWP 

MTO 

NWP 

XDRDR 

YWR 

ZVP 

KWR 

M TO 

NVQ 

IWE 

XDRBRE 

YVW 

ZTO 

KVW 

MTO 

NVW 

It  should  be  noticed  that  the  cross-coupling  terms  in  the  lateral 
channels  always  have  one  lateral  velocity  as  a  multipling  factor,  but 
the  longitudinal  terms  are  always  multiplied  by  two  lateral  velocities. 
This  is  why  longitudinal  variations  do  not  affect  the  lateral  channel 
when  the  lateral  velocities  are  zero,  but  any  turning  rate  'will  affect 
the  performance  of  the  longitudinal  channel. 

It  was  established  above  that  the  NSRDC  equations  are  accurate  for 
the  simulation,  of  operational  submarine  motion  when  they  are  used  with 
the  full  set  of  coefficients  supplied  by  this  facility,  Therefore,  a 
certain  amount  of  program  complexity  has  already  been  eliminated  during 
the  analysis  performed  at  NSRDC.  The  terms  which  are  omitted  can  be  of 
three  types.  The  first  type  are  those  for  which  the  coefficient  is  zero 
because  of  geometrical  symmetry  of  the  hull  or  mathematical  symmetry  of 
the  equations  of  the  hydrodynamic  reactions  in  Taylor’s  theorem  expansion. 
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XUU  YQR  ZWW  KPQ  MJAQ  NSTR 

YSTR  ZAW  KSTR  MAW 
XPAP  KVW  WWW 

The  second  type  are  those  which  can  exist  but  have  been  determined 
to  be  too  small  for  submarines  in  general.  They  include* 

YVQ  ZPP  KDR  MPP  NQR 

IWR  KR  NRAR 

KVQ  NWR 

KWR 

finally,  there  are  those  coefficients  which  have  been  set  to  zero 
because  they  have  not  been  investigated  and  the  effect  of  their  presence 
and  their  magnitude  is  unknown. 

The  terms  are  included  in  this  category.  They  cover  the  effect 
of  the  propeller  disturbing  the  flow  in  addition  to  producing  thrust  in 
the  body  of  the  submarine.  As  long  as  a  propeller  is  subjected  to  uniform 
axial  flow,  the  stream  velocity  is  constant  for  all  blade  positions  around 
its  path  of  rotation.  As  a  result,  the  angLe  of  attack  of  the  blade,  the 
thrust,  and  torque  are  constant,  and  the  ^  coefficients  can  be  set  to 
zero. 


As  a  rule,  stem  propellers  must  operate  in  non-axial  flow  conditions. 
This  is  true  ever  for  an  axially  mounted  propeller  on  a  body  of  revolution 
where  the  flow  is  symmetrical  and  converges  around  the  afterbody.  Then 
too,  the  body  axis  may  not  be  parallel  to  the  direction  of  the  stream 
velocity,  resulting  in  further  deflection  of  the  flow.  The  blade  angle 
of  attack  is  altered  with  respect  to  flow  direction,  thus  a  greater  thrust 
is  produced  at  one  side  resulting  in  an  upward  force  being  applied  to  the 
shaft.  Any  change  in  the  ship's  direction  of  motion  will  be  reflected  by 
the  hydrodynamic  forces  and  moments.  A  change  in  the  propeller  RPM  will 
result  in  a  changed  force  diagram.  The  propeller  blades  will  be  acting  at 
different  angles  of  attack  with  respect  to  the  f  luid  velocity  vectors 
until  steady-state  conditions  again  prevail  after  the  propeller  velocity 
is  changed. 

The  demonstration  submarine  is  known  to  have  certain  non-linear 
effects  that  are  not  reflected  in  the  equations  of  motion  with  the 
present  set  of  coefficients.  These  effects  could  be  provided  through  the 
addition  of  the  terras,  but  their  magnitudes  are  not  known  at  the 
present  tis».  However,  the  NSRDC  equations  of  motion  report  states  that 
for  the  moderate  changes  in  ahead  speed  involved  in  most  normal  maneuvers, 
all  of  the  -1)  terras  usually  can  be  neglected.  Therefore,  the  following 
coefficients  are  set  to  zero  under  the  assumption  that  axial  flow  exists 
at  the  propeller  at  all  times. 
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XWE 

IRE 

ZQE 

KSTRE  MQE 

NRE 

XWWE 

YVE 

ZWE 

MWE 

NVE 

XDRDRE 

YVAVE 

TWAWE 

MWAWE 

NVAVE 

XDSDSE 

YDRE 

ZDSE 

MESE 

NDRE 

However,  all  o.f  these  terms  are  Included  in  the  programing  so  that 
if  the  necessity  of  using  them  should  arise,  they  can  be  added  merely  by 
placing  a  number  on  an  input  card. 

A  survey  was  made  of  four  different  submarine  simulations,  and  in 
each  case  the  coefficients  listed  above  vers  not  used.  This  survey  was 
correlated  with  interviews  of  personnel  who  had  received  training  in  the 
simulators  in  question  and  had  served  on  the  submarine  class  being  simu¬ 
lated.  The  conclusion  reached  was  that  the  simulation  was  a  good  repre¬ 
sentation  under  normal  training  maneuvers  even  with  these  terms  missing. 

3.  INERTIAL  TERMS 

It  was  stated  that  the  equations  of  motion  can  be  derived  from 
d’Alembert's  principle  that  the  applied  forces  must  be  in  equilibrium 
with  the  inertial  reactions  of  the  body.  This  section  examines  these 
inertial  reactions  of  the  body  with  respect  to  submarine  simulation. 

The  inertial  terms  are  those  terms  that  result  from  the  mass  of  the 
vehicle  rather  than  from  any  outside  forces  on  the  vehicle.  As  such, 
they  are  written  on  the  left  hand  side  of  the  equation  in  the  standard 
equations  of  motion  published  by  NSRDC.  These  terms  are  of  two  types? 
those  resulting  from  mass  and  inertia  of  the  hull,  and  those  that  appear 
in  the  equations  because  the  center  of  gravity  does  not  remain  at  the  center 
of  the  coordinate  system  chosen  to  specify  the  coefficients. 

Generally,  the  coordinate  system  used  in  aircraft  simulation  is 
coincidental  with  the  center  of  gravity  and  moves  along  with  it.  This 
change  does  not  affect  the  coefficient  because  the  shift  is  small  and  is 
mainly  due  to  depletion  of  on-board  fuel  located  in  tanks  close  to  the  CG. 

An  entirely  different  case  is  found  in  submarine  simulation.  Large 
weights  of  water  are  taken  on  board  in  tanks  far  removed  from  the  origin 
of  the  reference  point  used  for  the  inertial  frame.  Therefore,  the 
parameters  Xq  ,  »  and  2<S  representing  the  shift  of  the  CG  off  the  origin 

are  included  in  the  simulation. 

The  first  step  was  to  eliminate  functions  that  are  zero  because  of 
the  nature  of  the  shape  of  a  submarine.  It  car.  be  assumed  that  Yq  can 
be  set  to  zero  because  the  CG  is  not  going  to  shift  far  in  the  lateral 
direction  as  the  ballast  tanks  are  filled.  These  tanks  are  symmetrical 
with  reference  to  the  X-Z  place.  Iy  has  the  same  value  as  and  the 
cross  inertial  2*y  ,I*z  >  andl^are  zero  on  all  submarines  studied.  This 
reduces  the  inertial  terms  of  the  equations  to; 

Axial  ->n(d-vr  +  w^-  Xq  (?**ra)+Z<5(p'r+ 
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Lateral  -  Jn  ($~Wp  -+  U.r  +  Lqfy'r-p)  -f*  ^qOfP^  Y-)) 

Normal  -  *m  -  U  v  p>  ~ZqCT,i+v2)4  Xq(yp* 

Ron  w^+u  r) 

Pitch  -X^  +  (lx-Iz)yp+>}(Z<j(u~Yr+w£)-X(;(w-U2~*vp)) 

Yaw  -Ixf  +  (i'f  ~Tx)7°  f  +  ^Xq($-w  p-fr  uy) 

The  various  values  of  cross  velocities  were  kept  in  the  reduced 
equations  of  motion  without  investigating  their  magnitude  of  effect. 

This  is  because  all  of  these  terms  are  needed  to  calculate  the  forces 
from  the  coefficients,  and  therefore  they  can  be  included  at  the  cost 
of  a  single  addition.  Also,  any  one  of  these  terms  can  become  quite 
large  under  casualty  conditions  and  will  therefore  be  needed  for  good 
recovery  simulation. 

If  Xq  and  Zq  can  be  held  at  zero,  as  in  an  aircraft  simulation,  the 
inertial  terms  could  be  greatly  reduced  in  complexity.  Therefore,  an 
investigation  was  made  into  the  effect  of  Xc  and  J_c. on  a  typioal  simu¬ 
lation.  M 

The  first  task  is  to  determine  the  range  of  Xq  and  Zq  •  The  equa¬ 
tions  for  these  parameters  are  ^  v 


where  the  Vs  are  the  weights  of  the  submarine  and  water  in  the  tanks  and 
the  2*3  and  X*a  are  the  coordinates  of  the  respective  water  ballast  tanks. 
The  figures  used  for  the  investigation  are  contained  in  the  report 
"Measures  of  Performance  Based  on  Motion  Simulation"  for  the  demonstration 
submarine.  In  the  case  of  flooding,  a  compartment  filled  with  water  can 
be  treated  as  an  additional  tank. 

Program  EGU30  calculated  X5  ,  Yq  ,  Z q  ,  and  W  for  any  tank  arrange¬ 
ment.  The  data  on  the  submarine  and  the  individual  tanks  is  used  as  an 
input,  and  any  tank  arrangement  is  controlled  by  an  input  card  following 
the  data.  In  one  of  the  submarines  studies,  there  were  values  of 
associated  with  individual  tanks,  but  the  shift  in  Yq  was  very  small. 

Also,  where  tanks  are  split  laterally,  the  two  sides  are  always  filled 
equally,  reducing  y^  to  zero. 
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A  number  of  different  tank  arrangements  were  tried  to  cover  a  range 
of  Xg  and  z<*._  The  final  values  used  were2<2  closest  to  the  center  of 
buoyancy,  a  normal  operating  condition  (W=B  an a  equal  distribution  in  the 
trim  tanks),  and  a  worst  case  inX^  with  major  Hooding  at  the  forward 
end  of  the  submarine.  Tests  were  made  at  5  and  25  knots  to  cover  the 
operating  range  under  these  conditions. 

A  normal  overshoot  maneuver  is  used  to  check  the  effects  of  varying 
X<^  and  Zfi .  Programs  EB920  or  ZC790  are  used  with  the  different  values 
and  u,  theta,  and  z  are  recorded,  The  results  are  given  in  figure  ii. 


T  T  T 


l  ~Ref.>Xq  O 

2- Z<j 

3-  Mown* l 

4- X<s  Maxrmwm 

.Figure  k.  Overshoot  Response  to  Tankage  Variation 

It  should  be  noted  that  these  curves  do  not  represent  the  true 
response  of  the  submarine  from  the  change  in  the  center  of  gravity  due 
to  distribution  of  ballast.  This  is  because  onlyX<^  andjTfi  is  changed. 
The  weight  still  equals  the  buoyancy. 


It  can  be  seen  that  at  25  knots  there  is  very  little  change  in 
response  to  changing  the  two  parameters.  However,  at  5  knots,  the  change 
is  drastic?  in  fact,  the  pitch  angle  effect  reverses  with  the  change  in 
Xq  .  The  change  in/^is  not  as  great. 

Program  ZC691  was  used  to  find  the  amount  of  error  in  three  oases. 
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TABLE  5.  CHANGE  OF  PERFORMANCE  WITH  TANKAGE  CHANGE 


PERCENT  CHANGE  OF  VARIABLE  { 

AU 

A  6 

AZ 

5  KNOTS 

ZG  MAXIMUM 

“1.3 

-6.6 

-23.4 

normal 

-372.6 

-681.2 

-  1307*7 

Xg  maximum 

-692.1 

-856.8 

25  KNOTS 

-  1  545.0 

Z$  MAXIMUM 

"2.9 

O.fc 

7.7 

NORMAL 

- 1 1*6 

-16.4 

— 4l  *4 

*6  MAXIMUM 

—16.0 

-20.7 

-  53.5 

Z  <5  could  be  sat  to  zero  if  it  ware  not  for  the  23. change  in 
maximum  depth  at  5  knots.  The  direction  of  change  of  z  remains  the  same. 

A  consideration  of  the  results  leads  to  the  conclusion  thatX^  and 
Z<^  should  be  kept  in  the  modified  equations  of  motion  set,  and  Zq  can  be 
set  to  a  ccr  w-.nt  in  the  basic  set. 

The  degree  of  variation  in  andZ^  can  be  used  to  determine  how 
accurate  the  tank  representation  shculd  be  in  a  submarine  casualty  control 
trainer.  The  parameters  were  varied  and  changes  measured  at  the  end  of 
one  minute.  A  15%  change  inZ^  resulted  in  no  change  in  u,  less  than  l£ 
change  in  the  response  uf  theta,  and  about  a  2%  change  in  maximum  depth 
change.  A  change  in  resulted  in  as  great  a  change  in  u,  and  a  3% 
change  in  theta.  Z  was  changed  very  little.  Therefore,  a  very  close 
simulation  of  in  the  tank  representation  is  unnecessary  even  in  a 
casualty  trainer,  but  the  factors  that  go  to  make  up  Xa should  be  computed 
to  within  of  the  actual  reeponse  in  both  time  and  weight. 

The  standard  equations  of  motion  allow  for  the  center  of  uuoyancy 
(GB)  to  be  displaced  from  the  center  of  the  coordinate  axis  as  well  as 
the  center  of  gravity.  These  parameters  are  expressed  as  functions  of 
Xg  ,Yg  , Zb  »  which  are  the  distances  in  each  dimension  of  the  dis¬ 
placement.  The  CB  is  fixed  in  the  submarine,  and  the  coefficient  co¬ 
ordinate  system  is  always  chosen  so  thatXg  and  Yg  are  zero.  Zg  has 
to  have  a  definite  value  that  puts  the  CB  above  the  CG  when  the  submarine 
is  submerged.  This  causes  the  mass  to  act  as  a  stable  pendulum  which 
returns  to  neutral  at  the  natural  frequency  of  the  system  when  the  outside 
disturbing  forces  are  removed.  The  restoring  force  is  the  force  of  gravity 
and  the  equations  of  motion  become 

Axial  -  (W~  B)  SIN  © 

Lateral  (w-8)  COS  O  SIN  <{> 

Normal  (W-B)  cos  ©  COS$ 


Lo 
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Roll  *“  (Zq  Zg  8)  CoS  0  S  iKi  ^ 

Pitch  — X^W  cos  0  co  s  (j)  -CZ<$  W-  ZgB)  SIN0 

law  X^  W  cos  ©  SIN  4* 

with  the  Xb  >Yq  t  an<?-Y<5  *®t  to  zero.  It  can  be  seen  thwt  time  will  be 
saved  if  terms  such  as  McosGsin0B  are  computed  only  once  and  saved  for 
use  in  the  valuation  of  the  accelerations.  This  also  applies  to  «feW fyd". 
The  multiplication  does  not  have  to  be  made  when  separate  tanks  are 
being  simulated, 

ZqW=  WfX*  +-^Wt  1  *  *T  i  .  w 

-  WoXo  +  £  Wt^  * Y [_ 

and  the  latter  teirm  can  be  used  before  the  division  by  W  is  made, 
h.  LONGITUDINAL  COEFFICIENT  REDUCTION 

Although  different  sets  of  operations  were  followed  for  each  group 
of  coefficients,  a  general  reduction  procedure  was  followed  that  was 
common  to  all  groups,  Hie  task  was  broken  up  into  three  parts:  first, 
an  analysis  to  determine  the  affect  of  each  coefficient  on  the  quality 
of  simulationj  second,  the  determination  of  the  effect  of  removing  groups 
of  coefficients;  and  finally,  the  degree  to  which  the  remaining  coefficients 
should  be  changed  in  order  to  improve  the  output  response  of  the  remaining 
set  of  coefficients  after  the  reduction  is  made, 

A  data  platform  is  established  by  treating  a  set  of  input  conditions 
that  will  vary  those  parameters  that  are  multiplied  by  the  coefficients 
of  interest.  The  full  set  of  coefficients  is  used  to  run  a  program  over 
a  sufficiently  long  period  of  time  to  set  the  trend  of  the  outputs  to 
tills  set  of  input  conditions.  Points  of  interest  occur  where  pitch  changes 
direction  or  where  most  of  the  variables  have  reached  a  steady-state 
condition.  Most  of  these  runs  are  about  one  or  two  minutes  simulated 
time  in  length.  The  variables  used  in  the  verification  nans  should  be 
those  which  the  trainee  will  observe  during  trainer  operation  such  as 
pitch  angle,  bank  angle,  and  forward  velocity. 

Hie  simulation  program  is  run  a  number  of  times  with  the  same  input 
conditions  and  the  coefficients  are  set  to  zero  on  a  term  by  term  basis. 

The  percentage  change  is  calculated  at  the  point  of  interest  and  a  table 
is  sat  up  so  that  the  effect  of  removing  any  one  coefficient  can  be  seen 
at  a  glance. 

Some  of  tiie  coefficients  are  not  tested  because  the  nature  of  the 
program  requires  a  numerical  value  or  else  the  output  will  not  resemble 
any  vehicle,  much  less  a  submarine. 

Hie  progr-  *  $  rerun  with  all  coefficients  having  small  percentage 
errors  set  to  WvO,  Ine  usual  effect  is  to  have  all  errors  increased  by 


lil 


iri  ts-**  4  /  O  n  ■** 

ns  vrmiur<voc*fl  du'^w^u** 


more  than  the  amount  of  the  separate  runs  added  together  because  of  the 
interactions  contained  in  the  program.  Terras  are  added  until  an  acceptable 
level  of  error  is  reached  or  remaining  terms  are  used  to  compensate  for 
the  changes. 

Compensation  by  adjusting  the  remaining  coefficients  seems  attractive 
at  first  glance,  but  it  cannot  be  carried  too  far.  This  is  logically 
obvious  because  of  the  amount  of  work  that  has  been  done  on  the  original 
equations  of  motion.  Any  term  that  could  be  combined  with  another  term 
over  the  full  range  of  operation  of  the  vehicle  would  have  already  been 
combined  with  that  term  when  the  coefficients  were  supplied  by  NSRDC, 
However,  the  requirements  for  a  training  simulator  are  not  as  rigorous 
as  those  needed  for  research,  so  simplifications  are  possible. 

In  the  sections  that  follow,  the  coefficients  will  be  divided  into 
natural  categories  and  the  method  taken  to  reach  a  successful  conclusion 
will  be  outlined.  What  cannot  be  shown  in  a  report  is  that  coefficient 
reduction  is  an  art  rather  than  a  science.  The  basic  assumptions  and 
conclusions  will  hold  true  for  any  submarine  chosen  for  simulation,  but 
a  cut  and  try  procedure  will  be  necessary  before  it  is  proven  that  the 
exact  reduction  used  is  effective  in  any  particular  case. 

The  longitudinal  coefficients  can  be  divided  into  two  classes  after 
the  coefficients  that  have  been  previously  set  to  zero  are  eliminated. 

They  are  those  which  are  necessary  for  the  simulation  to  operate  and 
those  which  are  to  be  investigated.  The  two  classes  are  given  in  table  6. 


TABLE  6.  NON-ZERO  LONGITUDINAL  COEFFICIENTS 


Axial 

Normal 

Pitch 

Necessary 

XUD 

ZQD 

MQD 

XDSDS 

ZWD 

MWD 

XD3DB 

ZQ 

MQ 

ZSTR 

MS7R 

ZW 

KW 

ZDS 

MDS 

ZDB 

WEB 

To  be  investigated 

XQQ 

2AQDS 

HAQDS 

XWQ 

ZWAQ 

MAWQ 

xww 

ZWAW 

MWAW 

The  velocity  and  acceleration  coefficients  are  required  to  give 
stability  to  the  differential  equation  solutions,  ZSTR  and  M5TR  are 
required  to  give  correct  solutions  to  steady-state  conditions  when  the 
submarine  is  in  level  flight.  The  effect  of  eliminating  them  is  discussed 
later.  The  rest  are  required  so  that  operation  of  the  planes  will  result 
in  response  of  vehicle  movement  in  the  proper  direction. 

This  does  not  seem  to  leave  many  coefficients,  but  if  the  remaining 
ones  are  examined  it  will  be  seen  that  eliminating  them  will  reduce  the 
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complexity  of  the  equation  to  a  rauoh  greater  degree  than  their  number 
would  indicate.  Let  ua  look  at  the  ZWAQ  as  a  horrible  example.  It  is 
multiplied  by  \(v* which  requires  five  multiplications, 

a  division,  a  square  root,  and  three  absolute  values* 

A  set  of  initial  conditions  has  to  be  established  before  investigation 
of  the  longitudinal  loop  can  begin.  The  natural  set  would  be  to  establish 
leve?.  flight  before  each  rur..  This  is  not  as  easy  as  it  sounds,  and  a 
side  investigation  was  made  inco  the  subject  of  critical  speeds  and  neutral 
angles  for  a  submarine.  Neutral  angles  are  defined  as  those  steady  pitch 
angles,  sailplane  angles,  and  sternplane  angles  which,  when  combined, 
produce  a  condition  of  steady  level  flight  or  zero  path  angle  at  a  particular 
speed*  This  requires  that  the  #angla  of  attack  equal  the  pitch  angle,  and 
that  the  accelerations  w  and  equal  zero.  Reference  is  made  to  angle 
of  attack  in  many  texts,  including  the  NSRDC  reports,  but  the  parameter 
does  not  appear  in  the  equations  of  motion*  This  is  because  a  body  axis 
system  of  coordinates  is  used  instead  of  a  flight  path  set.  The  angle  of 
attack  is  the  angle  between  the  flight  path  and  the  body  when  roll  angle 
is  zero. 


Figure  5,  Angle  of  Attack 

The  angle  of  attack,  and  in  the  case  of  steady  level  flight  the 
pitch  angle,  can  be  calculated  from  the  relationship 

U*  =  U?  +  uu  a 

5IN  ©  3 1 hi  “  ‘Qp  —  y  X 

as  seen  in  figure  5*  This  removes  the  parameter  0  by  replacing  it  with 
a  function  of  w.  For  any  speed  u,  the  equations  of  motion  can  be  solved 
for  level  flight  in  terms  of  w  and  or  by  setting  w,  and  q  to  zero. 

From  No-Tael  Equation  (3) 

Zw  |wl  *W|W|+Zwvv -W  2+Zw  •  U-vu  +ZM  •  ‘V  M+Z*  •  uSZf-uSS  *Q 

From  Pitch  Equation  (5) 


U3 


P  =  density  of  water,  1  ■  length  of  submarine 

&  can  be  either  $y  or  5 5  provided  that  the  correct  coefficients 

Zfb  ,M$fcor7$s>MSs  are  used  as  inputs. 

Program  EGU70  was  written  to  solve  the  equations  for  w  and  £  , 

This  program  utilizes  all  coefficients  of  interest,  so  that  it  can  be 
used  with  any  set  of  coefficients  applicable  to  the  NSRDC  equations. 

The  pitch  equation  is  solved  for  £  with  a  trial  value  of  w,  and  this 
value  is  used  in  the  normal  equation  via  Newton's  method.  Once  w  and 
5  have  been  calculated,  the  program  calculates 

©  —  til 71"'  <—£  radians 

©«»=  57-3  ©  degrees 

<fj=  57.3  6  degrees 

and  prints  out  u  in  ft. /sec.  and  knots  as  well. 

This  program  serves  another  purpose  in  addition  to  determining  the 
initial  conditions  for  longitudinal  runs.  The  values  for  S  are  determined 
for  the  operational  speed  range  of  the  submarine.  At  some  point,  the 
values  become  very  large.  This  speed  is  known  as  the  critical  speed 
because  at  this  point  the  controls  are  ineffective  in  controlling  the 
pitch  attitude  of  the  submarine.  Fortunately,  the  critical  speed  is 
different  for  the  stemplanes  and  the  sailplanes.  The  sailplanes  are 
effective  at  low  speeds,  below  ten  knots.  The  sternplanes  have  a  very 
interesting  effect  that  must  be  simulated  correctly  in  order  not  to 
introduce  negative  training  into  the  simulator.  At  a  certain  low  speed, 
the  neutral  angles  suddenly  become  very  large.  Pitch  angle  and  stemplane 
deflection  reverse  direction  over  a  very  narrow  range.  Just  bclcw  this 
speed,  the  angles  come  back  to  approximately  their  former  position..  The 
"Model  Investigation  of  the  Stability  and  Control  Characteristics"  report 
for  the  demonstration  submarine  has  a  graph  of  this  behavior.  It  was 
used  to  check  the  steady-state  angles  generated  by  the  simulation  program 
when  using  the  original  coefficients.  There  was  a  perfect  fit  between 
program  EC^70  output  and  the  NSRDC  data  for  both  6b  and  $S  over  the 
operational  range  of  the  demonstration  submarine. 
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The  low  speed  critical  area  is  vital  to  effective  submarine  training. 
Therefore,  if  any  coefficient  changes  affeot  this  determination  to  any 
great  degree,  they  must  keep  their  original  value  in  the  modified  simulation. 
The  basic  simulation  is  not  as  critioal  since  a  particular  submarine  is  not 
being  simulated.  This  is  -eason  for  keeping  ZSTR  and  MSTR  on  the  list  of 
necessary  coefficients.  ZWa.  nd  MWAW  should  be  cheoked  closely  because 
they  have  a  definite  effect  on  the  values  produced  by  ECii70.  ZW,  ZAW,  MWW, 
and  MAW  have  already  boon  sot  to  zero  in  the  original  analysis  by  N5RDC. 


It  was  hoped  that  a  new  set  of  linear  coefficients  could  bo  generated 
by  means  of  impulse  testing,  but  this  did  not  work  in  the  case  of  the 
longitudinal  loop.  Therefore,  a  series  of  overshoot  runs  were  made  to 
determine  the  effect  of  setting  all  coefficients  in  the  second  class  to 
zero.  This  input  was  the  standard  overshoot  schedule. 


Figure  6.  Comparison  Criteria 

The  comparison  procedure  used  for  the  longitudinal  channel  uses  the 
maximum  change  in  pitch  angle  and  depth  as  the  reference  points.  In 
order  to  remove  any  bias  due  to  the  particular  overshoot  maneuver  used, 
all  percentages  are  calculated  from  change  of  the  change  at  the  reference. 
If  A0  -full  set  and &Q  -modified  set  are  compared,  the  percentage 
change  is 


°/o  A© 


a©fs 


X  IOO 


The  same  holds  true  for  changes  in  z.  Program  simulation  times  are  set  to 
be  large  enough  so  that  the  reference  points  are  always  reached  for  all 
runs  as.  shown  in  figure  6. 

A  large  number  of  runs  were  made.  Table  5  gives  several  represent¬ 
ative  cases. 

ZWAQ  is  an  extremely  desirable  term  to  remove,  and  it  seems  to 
compensate  nicely  for  MAQDS  and  ZAQDS.  Therefore,  these  are  set  to  zero 
in  both  programs. 
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TABLE  7.  PERCENT  CHANGE  IN  OVERSHOOT,  COEFFICIENTS  SET  TO  ZERO 


Coefficient  *  0 

% 

$ 

ZAQDS 

2.3 

O.lj 

ZWAO 

2.3 

0.0 

ZWAW 

ii.6 

0.1 

MAQDS 

-2.3 

-1.8 

MAWQ 

U.6 

0.3 

MWAW 

7.0 

3.3 

ZAQDS 

O.li 

ZWAQ  =  0 

MAQDS 

0.0 

All  zero 

20,9 

13.6 

All  zero 

ZW  ♦  75$ 

0 

0.2 

MW  -  11$ 

of  ZWAW,  MAWQ,  and  MWAW  and  that  ZWAW  and  MWAW  are  two  of  the  coefficients 
that  determine  the  critical  speed,  it  was  decided  that  these  three  would  be 
kept  in  the  modified  set  of  equations.  This  does  not  increase  the  amount 
of  computation  by  a  large  amount  once  the  factor  has  been 

calculated.  Since  it  was  found  necessary  to  calculate  the  square  root  for 
another  ‘reason,  the  most  these  three  terms  will  add  is  six  multiplications, 
three  additions,  and  three  data  storage  locations  for  the  coefficients. 

There  is  another  extremely  important  reason  for  keeping  ZWAW,  MWAW, 
and  MAWQ  in  the  modified  simulation.  These  terms  govern  the  action  of  the 
submarine  if  it  is  going  straight  up.  This  would  not  happen  under  ordinary 
maneuvering  procedures,  but  during  casualty  conditions  it  is  possible  to 
blow  the  main  ballast  tank  while  sitting  dead  in  the  water.  This  causes 
an  upward  movement  in  which  w  is  the  only  body  velocity  to  have  any  magni¬ 
tude.  The  conditions  are  equivalent  to  an  angle  of  attack  of  ninety  degrees 
at  this  time. 

The  impulse  runs  did  show  that  it  is  possible  to  produce  the  same 
output  with  a  linear  set  of  coefficients  even  if  they  could  not  be  deter¬ 
mined.  Therefore,  all  the  terms  were  set  to  zero  in  spite  of  the  large 
errors,  and  the  remaining  coefficients  were  varied  on  a  cut  and  try  basis 
to  see  if  the  errors  could  be  reduced.  Runs  showed  that  the  acceleration 
coefficients  ZWD,  ZQD,  etc,  had  very  little  effect  on  the  magnitude  of 
pitch  angle  and  depth  and  only  varied  the  frequency  response,  MQ  and  ZQ 
were  a  little  more  sensitive  and  the  greatest  response  was  to  changing 
MW  and  ZW.  Seven  runs  determined  that  ZW  should  be  increased  by  75$  and 
MW  was  decreased  by  11$.  The  magnitude  of  the  response  was  excellent  and 
the  timing  was  not  greatly  altered.  Additional  runs  shoved  that  XQO,  XWQ, 
and  XWW  did  not  affect  either  longitudinal  loop  or  the  acceleration-deac- 
celeration  response  whan  thrust  was  changed.  Therefore,  these  coefficients 
were  also  set  to  zero. 
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5.  LAtkkaL  COEFrJLuIEra  REDUCTION 

The  lateral  coefficients  can  be  divided  in  the  same  manner  as  the 
longitudinal  coefficients. 


TABLE  8.  NON-ZERO  LATERAL  COEFFICIENTS 


Lateral 

Roll 

Yaw 

Necessary 

YRD 

KRD 

Nttl) 

YPD 

KPD 

NPD 

YVD 

KVD 

NVD 

YR 

KP 

NR 

YP 

KV 

NP 

YV 

NV 

YDR 

NOR 

To  be  investigated 

YARDR 

KPAP 

NARDP. 

YVAR 

KVAV 

NAVR 

YVAV 

NVAV 

Again,  the  number  of  coefficients  that  can  be  eliminated  seems  to 
be  small,  but  removing  YVAR  would  help  in  particular  with  simplifying 
the  program. 


The  lateral  loop  is  more  complex  than  the  longitudinal  loop  because 
three  directions  are  interacting  rather  than  just  two.  In  spite  of  this, 
impulse  testing  provided  information  of  value  unlike  the  other  case. 

The  procedure  is  as  follows,  A  run  is  made  using  the  full  simulation 
program  EB920.  A  new  set  of  linear  coefficients  is  generated  from  the 
assumption  that  the  non-linear  terms  can  be  replaced  by  adjustments  in 
the  remaining  coefficients.  The  EBpPO  run,  made  under  autopilot  control 
with’  a  rudder  impulse  at  t=0,  is  used  to  determine  an  average  constant 
value  of  v  and  r  under  th8  assumption  that  w  is  zero.  These  values  are 
used  in  program  EC320,  lateral  characteristic  equation  calculator,  to 
solve  for  the  adjusted  coefficients  and  natural  frequencies  using  the 
principles  of  linear  servo  theory.  The  process  is  identical  to  that 
explained  in  the  section  on  longitudinal  impulse  testing,  except  that 
the  lateral  characteristics  equation  has  four  roots  instead  of  three. 

Two  of  the  roots  are  real  so  that  time  response  to  an  impulse  becomes 


<t> (t)  =  a. ,e  cos (fi t  +  a,  e~ st 


The  additional  term  over  the  longitudinal  case  is  a  very  slowly 
converging  spiral  mode  that  is  present  in  almost  all  six-degree-of- 
freedom  vehicles.  Pilots  are  generally  unaware  that  this  mode  exists, 
because  they  automatically  correct  bank  angle  for  it  while  piloting  the 
vehicle.  In  aircraft,  this  mode  often  diverges  slowly  and  will  cause 
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entrance  to  e  spin  if  the  pilot's  wind  wanders  for  any  length  of  tiros* 

Cards  are  punched  with  the  values  of  bank  angle  and  time  and  a 
convergence  run  is  produced  so  ths  actual  frequency,  rather  than  the 
estimated  one,  can  be  used  to  plot  the  output.  The  values  of  estimated 
frequencies  produced  by  SC32C  along  with  the  cards  are  used  by  pregram 
EC120,  tins  response  coefficient  estimator,  to  calculate  input  values 
for  Brown's  convergence  routine  (program  EC310)  for  eight  variables  instead 
of  six,  EC120  will  take  o<  j  /$  j  2T  -  and  £  calculated  by  EC320 
and  output  j  ^  1  j  and  ip  as  'defined  by  the  time  response 
equation.  Generally,  Brown's  routine  is  unreliable  when  eight  variables 
are  to  be  optimized  so  the  estimated  values  for  ^3  and  £  are  held 
constant  for  the  first  convergence  run  while  optimizing  the  other  six, 

EC320  also  produces  adjusted  values  for  YV,  YP,  KV,  NN,  and  NR, 

These  values  are  used  in  the  coraple te  simulation  program  with  YARDF, 

YVAR,  YVAV,  KPAP,  KVAV,  RARDR,  NAVR,  and  N7AV  set  to  zero.  Impulse  runs 
were  made  at  5,  10,  and  25  knots,  and  the  results  were  plotted  against 
the  full-set  runs  determined  above.  The  ourves  show  excellent  agreement 
between  the  two  impulse  runs  at  5  and  10  knots  and  good  agreement  at  25 
knots.  This  states  that  ths  non-linear  coefficients  can  be  set  to  zero 
and  the  lateral  loop  will  still  have  the  same  natural  frequencies  and 
dancing  characteristics  as  with  the  full  set  of  coefficients.  Therefore, 
it  was  assumed  that  this  procedure  would  lead  to  the  basic  set. 

It  was  discovered  later  on  that  this  approximation  was  too  drastic, 
when  all  the  modifications  were  made  to  the  modified  and  basic  programs. 
Additional  tests  were  made  under  closed  loop  control  to  simulate  the 
response  when  the  pilot  controls  the  vehicle. 

The  trouble  with  making  this  type  of  investigation  on  a  general 
purpose  digital  computer  is  that  it  is  not  a  simulator.  There  is  no  way 
to  have  a  man  monitor  the  output  and  change  £s  »  for  example,  as  the 
program  proceeds.  A  set  of  input  conditions  is  established  and  the 
program  calculates  to  the  end  of  a  set  time.  In  order  to  simulate  closed 
loop  response,  an  autopilot  was  designed  to  keep  the  submarine  in  level 
flight  at  a  constant  depth.  Otherwise,  the  longitudinal  loop  was  free 
to  move  to  any  position  during  the  course  of  a  steady-state  turn.  The 
equation  used  for  the  autopilot  was 


D S  -  DB=  .008(zC-Z)+5.5©  +.01  sm©-  cos©)+2 .0^ 


DS  and  DB  are  the  stern  and  sailplanes  respectively  in  radians,  ZG  is 
the  command  depth  in  feet,  and  the  other  variables  are  the  same#as  in  the 
simulation  program.  The  term  (u  sin  9  -  w  cos  9)  was  used  for  z  because 
1,  was  not  available  in  the  CONTROL  subroutine.  The  autopilot  programing 
limited  both  plane  deflections  to  35°  in  either  direction. 

Closed  loop  control  was  found  to  be  extremely  necessary  to  determine 
that  approximations  made  to  reduce  the  coefficients  were  still  valid 
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when  tha  simulation  was  tested.  For  example,  when  the  longitudinal  loop 
was  frozen  (w  =  0,  6  a  o)  the  lateral  loop  responded  veiy  closely  to  the 
full  set  of  coefficients  to  constant  rudder  commands  when  the  non-linear 
terms  were  eliminated.  This  was  not  the  case  when  closed  loop  control 
was  used.  Therefore,  all  tests  were  made  with  both  the  lateral  and  cross¬ 
coupling  terms  us^ng  the  autopilot  with  a  steady  turn  requirement*  Even 
though  the  submarine  was  being  held  in  level  flight,  the  behavior  of  the 
autopilot,  as  reflected  by  z  and  «,  can  be  used  tc  determine  the  effective¬ 
ness  of  coefficients  in  both  the  longitudinal  and  lateral  loops. 

The  degree  of  maximum  bank  angle  change  and  steady-state  bank  angle 
late  in  a  turn  were  used  as  a  standard  to  compare  the  effect  of  the 
various  coefficients  as  shown  in  figure  7. 


figure  7,  Bank  Angle  Comparison 

Table  9  shows  the  relationship  of  $  max  and  (j^  for  several  coef¬ 
ficients  to  the  rudder  deflection  that  was  used. 


TABLE  9.  BANK  ANGLE  COMPARISON 


Coefficient  -  0 

A<P  max? 

YARDR 

10c  0 

11.0 

YVAR 

11.8 

11.8 

YVAV 

13.7 

lh.7 

KPAP 

l.C 

3.2 

KVAV 

-65.0 

-70.0 

NARBR 

-  h.o 

5.o 

NAVR 

3.5 

-  3,5 

NVAV 

3.1 

-  2.0 

It  is  obvious  that  KVAV  is  necessary  if  the  maximum  and  steady  roll 
angles  are  to  resemble  the  original  effects  at  all.  KVAV  is  multiplied 
by  v*(v*4  w  %  as  in  YVAV  which  has  the  second  longest  magnitude 
change,,  Therefore,  XVAV  can  be  utilized  with  only  one  additional  multi¬ 
plication  and  addition  along  with  a  data  storage  location  for  YVAV.  NVAV 
was  kept  for  the  same  reason  although  its  affect  is  minor.  This  is  because 
it  can  be  added  with  very  little  effort  and  its  sign  is  such  that  it  aids 
in  stabilizing  YV  when  computing  v  *  This  could  be  important  during 
casualty  conditions  when  the  equations  are  making  large  changes. 

Therefore,  the  lateral  loop  is  reduced  by  setting  YARDR,  YVAR,  K?AP, 
NARSR,  and  NAVR  to  zero.  The  effects  of  this  did  not  cancel  out,  so 
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several  runs  were  made  to  adjust  KVAV  to  bring  maximum  bank  angle  and 
steady-state  bank  angles  into  line.  In  the  case  of  the  demonstration 
submarine,  KVAV  was  reduced  by  50#.  It  was  found  that  heading  angle 
was  increasing  at  a  faster  rate  than  originally  when  all  the  modifications 
were  made.  Therefore,  NDR  was  adjusted  on  a  cut  and  try  basi3  to  bring 
heading  rate  back  to  its  original  value  without  affecting  bank  angle, 

NDR  was  reduced  by  20#  to  achieve  this  effect.. 

Since  KVAV  is  necessary  and  is  to  be  adjusted  no  matter  what  happens 
to  the  other  coefficients,  it  was  concluded  that  this  lateral  set  should 
be  used  for  both  the  basic  and  modified  set  of  equations. 

6.  CROSS-COUPLING  COEFFICIENT  REDUCTION 

The  cross-coupling  coefficients  are  those  that  are  used  only  when 
all  six-degrees-of -freedom  are  energized.  They  are  listed  in  table  10. 


TABLE  10.  NON-ZERO  CROSS-COUPLING  COEFFICIENTS 


Axial 

Lateral 

Normal 

Roll 

Pitch 

Yaw 

XRR 

YPQ 

ZRR 

KQR 

MRR 

NPQ 

XRP 

YWP 

ZVR 

KWP 

MVR 

NWP 

XVR 

YVW 

ZW 

MW 

NVW 

XW 

ZRP 

MVP 

ZVP 

MRP 

All  of  these  were  set  to  zero  and  a  rumber  of  steady-state  turns 
were  made  under  autopilot  control.  It  was  soon  evident  that  something  was 
wrong  when  the  autopilot  could  not  hold  depth  to  within  several  hundred 
feet.  When  the  original  coefficients  were  used,  the  depth  change  wa3 
about  four  feot.  A  term-by-term  check  showed  what  the  difficulty  vast 


TABLE  *.  i .  DEPTH  CHANGK  FOR  STEADY  TURN 


Coefficient  =  0 

Change  in  2? ,  ft. 

None 

Ii,l6 

XRR 

5.29 

XRP 

5.21 

XVR 

9.31 

XW 

M3 

ZRR 

18.00 

ZVR 

12.38 

ZW 

-33.00 

ZRP 

5.21 

ZVP 

M3 

MRR 

-25.1*0 

MVR 

21.27 

MW 

21.26 

MVP 

h.28 

MRR 

3.85 

_  » 
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ZRR,  ZVR,  ZW,  MRR,  KVR,  and  MW  are  required  to  make  the  equation* 
behave  liku  a  submarine  during  steady-state  turns,  the  other  longitudinal 
cross-coupling  coefficients  and  all  of  the  lateral  0-C  coefficients  can 
be  set  to  zero.  Any  unbalance  in  the  lateral  loop  can  be  taken  oare  of 
with  the  KVAV  adjustment. 

7.  THRUST  TERMS 

thrust  is  computed  as  a  function  of  the  engine  parameters  in  the 
general  equations  of  simulation.  This  quantity  is  compared  to  the  coef¬ 
ficient  of  drag,  computed  from  aerodynamic  parameters,  in  the  axial  force 
equation.  When  thrust  equals  drag,  u.  is  zero,  and  the  forward  body 
speed  is  in  a  steady-state  condition,  A  slightly  different  method  is 
used  in  the  N5RDC  equations  of  motion  because  the  coefficient  of  drag 
does  not  appear  explicitly  in  the  axial  equations. 

Drag  coefficients  are  computed  as  functions  of  rotational  and  side 
velocities  through  the  water  and  as  functions  of  the  control  surfaces, 
but  all  of  these  terras  go  to  zero  when  the  forward  velocity  lies  along 
the  X-body  axis  (©4  :  0),  Tne  coefficient Xhm.  is  not  the  main  drag 
coefficient,  but  represents  small  departures  of  the  drag  coefficient 
from  an  equilibrium  value.  It  is  set  to  zero  in  the  NSRDC  equations. 

Therefore,  the  main  computation  for  the  comparison  between  thrust 
and  drag  falls  on  the  term 

i  />-£*(aiiC*  +  b;.U.U.e  +  CiU*) 


where  P 

SL 

OL 


at>bl)cl 


■  water  density 
s  submarine  length 
=  forward  speed 
«  command  speed 
s  thrust  coefficients 


Q.*  are  a  set  of  constants  that  are  used  to  represent 
propeller  thrust  and  body  drag.  They  are  computed  from  both  the  steady- 
state  thrust  conditions  and  the  acceleration-deceleration  values  for  the 
submarine.  Three  sets  are  supplied  with  the  submarine  coefficients.  One 
set  is  used  when  the  commanded  speed  is  half  of  the  actual  speed  or  greater. 
A  second  set  is  for  below  this  value,  and  a  third  one  is  for  full  reverse. 
The  three  coefficients  in  each  set  must  add  up  to  zero.  This  is  so  that 
when  commanded  speed  equals  actual  speed,  thiust  equals  drag  and  the  output 
of  the  above  terra  is  zero. 


Accaleration  and  deceleration  data  was  not  available  for  the  demon¬ 
stration  submarine.  Therefore,  it  was  assumed  that  the  NSRDC  coefficients 
accurately  represented  thrust  and  they  were  programed  as  indicated  in  the 
aquations  of  motion.  This  means  that  the  modified  and  basic  equation  sets 
are  compared  against  the  full  set  in  acceleration-deceleration  runs  instead 
of  using  actual  submarine  data  for  verification. 
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In  order  to  check  out  thrust  response,  program  2C300  was  written 
to  calculate  both  thrust  in  pounds  and  acceleration  in  feet  per  second 
squared  for  any  condition  of  actual  and  commanded  speed,  the  equation 
of  interest  is 

+  hi  UUC  +  c-k  *1) 

where  a  is  the  mass  of  the  submarine.  Tnu.  has  the  dimensions  of  pounds 
(English  units)  which  is  to  be  expected  from  F  =  ma.  Therefore,  the 
calculation  of  the  term  on  the  right-hand  side  of  the  equation  can  be 
printed  out  directly  as  pounds  of  thrust  and,  when  divided  by  the  mass 
cf  the  submarine,  as  the  acceleration  in  feet  per  second. 

The  equations  of  motion  use  the  parameter  **v,  tc  indicate  which 
set  of  coefficients  to  use  in  the  thrust  equation.  This  is  merely  the 
ratio  of  the  commanded  speed  divided  by  the  actual  speed  and  indicates 
the  breakpoints  on  the  curves  used  to  determine  the  coefficients.  The 
thrust  program  allows  limits  to  be  specified  during  input. 

The  commanded  speed  Uc  is  related  to  the  RPM  of  the  propeller  by 
a  linear  constant.  This  constant  can  be  determined  from  curves  which  are 
found  in  the  NSRDC  report  "Resistance  and  Propulsion  Characteristics  of 
the  .  Submarine  as  Predicted  from  Model  Tests*  for  the  particular 

submarine  of  interest. 

It  should  be  noted  that  Uc  can  go  negative  and  this  can  cause  **1 
to  be  negative.  This'  happens  when  backing  down  during  emergencies,  and 
the  coefficients  are  computed  correctly  for  this.  However,  that  does  not 
mean  that  U.  can  go  negative.  All  of  the  coefficients  assume  that  U.  is 
always  positive  and  the  thrust  terms  are  calculated  on  this  basis, 

Uc  can  be  calculated  from  RPM  times  a  constant,  but  some  additional 
subroutines  would  be  required  in  a  submarine  simulator.  The  time  history 
of  an  RPM  change  has  to  allow  for  communication  lag; ,  human  response  time, 
and  engine  response  lags.  This  could  be  simulated  by  a  delay  followed 
by  an  exponential  change  to  the  new  value  of  RPM. 

8.  KINEMATIC  RELATIONS 

The  kinematic  relations  can  be  considered  as  a  pair  of  transformations 
that  will  transform  the  body  axis  rates  into  Ruler  angle  rates  and  the  body 
axis  velocities  into  earth  axis  velocities.  Part  of  these  transformations 
are  given  on  page  12  in  the  NSRDC  Report  2 $10,  but  not  in  a  form  that  can 
be  used  conveniently  for  digital  computation. 

All  forces  and  accelerations  are  calculated  in  the  body  axis  system. 
This  system  lies  along  the  physical  axis  of  the  submarine  body  and  moves 
along  with  it.  These  are  not  the  axes  about  which  bank,  pitch,  and  yaw 
rotate,  and  therefore  a  transformation  is  needed  to  calculate  the  angular 
rates.  It  is  known  as  an  Euler  transformation  and  the  development  of 
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this  particular  set  is  given  in  many  text  books3”*  NSHDC  Report  2510 
~ives  the  following  expression  for  bank,  pitch,  and  yaw  rates. 


<J>  —  P  +  SJN  © 

X  _  *1+^  co$©  StKlC t> 

O  -  C^slj) 

*  _  Q  ^in  <1> 

^  cos  ©cos  $ 


The  trouble  with  this  set  of  equations  is  that  yaw  rate  has  to  be 
known  to  calculate  pitching  rate  and  pitching  rate  has  to  be  known  to 
calculate  yaw  rate.  If  the  body  axis  angular  rates  are  placed  on  one 
side  of  the  equation,  the  Euler  rates  on  tho  other,  and  matrix  notation 
is  used,  the  difficulty  is  obvious. 


X  O  SIN  © 

— 

—  — 

p 

O  CoS  <t>  COS©  SIN  0 

e 

% 

O  —SIN  0  COS©coS0 

y  _ 

'r 

The  above  equation  assumes  that  the  Euler  rates  are  known  and  the 
body  axis  rates  desired  when  the  reverse  is  true. 

The  inverse  matrix  is 


& 

4 

SIN  ©SIN  C|> 

SiN©coS0 

i 

CoS  © 

cos  © 

0 

O 

cos  0 

—  SIN  0 

SIN  <t> 

COS  4) 

o 

COS  © 

CoS  © 
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With  this  matrix,  p,  q,  and  r  are  calculated,  and  the  transformation 
is  aade  to  solve  for  the  Euler  rates.  An  additional  saving  in  complexity 
can  be  made  if  the  new  rates  are  calculated  in  the  proper  order.  Expand¬ 
ing  the  matrix  gives  ths  following. 


4>  --p  + 


=  7°  + 


SIN  ©  SIN  <t> 

—  ?  + 


SIN  ©  cos 

cos  ^  ^ 


_  !%  sin<!>  +  r  cos  <t>\ 

V - ST© - / 


—  <t>  V  cos(t>  _  %  SIN  <j>4y  CCS  4> 

^  "  CoS  O  COS  © 


cos© 


Therefore,  if  1ft  is  calculated  from  q  and  r  before  the  calculation 
for  0  ,  the  three  equations  can  be  written  as  follows. 


©  =2  «|C©S  <P~Y  SIN  <$) 

/  —  *is|n4>  +  >  c°s  4) 

•  COS© 


<f>  =  7O  +  SIN  0 


The  order  of  the  calculations  is  important  sincq  a  digital  computer 
can  evaluate  only  one  term  at  a  time.  Notice  that  0  is  the  same  as 
the  HSitDC  set,  but  that  the  other  taro  equations  can  be  calculated  inde¬ 
pendently.  Che  problem  with  any  Euler  transformation  is  that  singular 
points  can  occur.  In  the  second  equation,  the  determinate  of  the  matrix 
is  cosO,  ao  the  transformation  is  not  valid  at  0  =TT/2.  The  assumption 
is  made  that  although  large  pitch  angles  (0  >  60°)  can  occur,  a  submarine 
cannot  get  into  a  true  vertical  position.  Therefore,  the  singularity 
should  not  give  any  trouble.  As  a  practical  matter,  such  a  condition 
would  give  either  the  largest  possible  yaw  rate,  or  the  problem  would 
stop  due  to  a  zero  divide  check  error  depending  on  the  computer. 

Generally,  a  maneuvering  trainer  is  concerned  with  angles  and  depth, 
but  not  with  position  in  the  X-Y  plane.  Because  of  this,  the  NSRDC  report 
gives  an  equation  for  vertical  speed  only  related  to  the  earth  axis,  and 
ignores  the  equations  needed  for  X  and  Y  positions.  These  equations  were 
added  to  the  research  submarine  program  so  that  the  outputs  would  be 
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complete.  The  only  reason  for  including  them  in  a  submarine  simulator 
program  would  be  if  it  were  to  be  used  in  a  ballistic  missile  sub-trainer 
where  position  with  reference  to  distant  points  is  important.  The 
f enaction  used  changes  the  body  axis  velocities  to  earth  axis  velc*  uties 
with  x  along  the  north  axis,  y  along  the  east  axis,  and  s  downward.  The 
equations  are 


Z  =  —U.  SIN©  4“  V  coS0  SIN  0  4*  NMCoS  ©  C.OS  ^ 

X~  ut  cos  ©  cos  (p  4*  v(sin(J)ssw©  CoS  (^-coS^SM  (ft) 

4*  w(ccs  CpSJKiOcos  + 

Y  =  a  COS  ©  SIN  (ft  4*  V  (si  N  0  SIN©  SIN  (p  4*  c©S  ( {>CoS  (ft ) 

4-  w(«foS  4>  SIN  ©SIN  Ip  *“  $IN 


It  can  be  seen  that  time  and  space  will  be  saved  in  the  computer 
if  x  and  y  do  not  have  to  be  calculated. 

There  exist  fairly  simple  algorithms  for  calculating  sine  and  cosine 
functions.  The  simplest  one  of  all  is  suggested  for  the  basic  trainer 
if  core  space  problems  require  that  the  final  program  be  as  small  as 
possible.  They  are 

sin©  -  0.971  0  radians 

cos©  =  1.02  -  0.28  8  radians 

These  functions  have  less  than  $%  error  up  to  h2°.  Submarine  pitch 
angles  can  exoeed  this  limit,  but  it  should  not  occur  too  often  in  a 
basic  maneuvering  trainer. 

A  more  satisfactory  algorithm  for  the  modified  simulation  is 

s,N  ©  *  ©  ~  § 

©x  ©4 
cose  =  1  -  g  +  24 
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They  hare  lass  than  3%  error  at  75°  which  should  serve  for  any 
purpose.  Of  course,  both  work  only  over  the  range  of  Tf  /2  to  “TT/2 
so  sin  yu  and  cos  will  have  tc  be  calculated  by  a  more  coaplar.  sub¬ 

routine  if  X  and  Y  are  to  be  calculated, 

9.  NEAR  SURFACE  OPERATION 

Two  affects  come  into  play  when  a  submarine  is  operating  close  to 
the  surface.  The  coefficients  of  the  equations  of  motion  are  changed 
due  to  the  presence  of  the  surface  and  forces  and  momenta  are  generated 
by  waves  on  the  surface. 

in  investigation  was  aaue  into  the  effects  of  coefficient  changes 
on  the  demonstration  submarine.  These  coefficients  were  supplied  by 
MSKDC  as  a  result  of  depthkeeping  studies  under  a  seaway  and  are  assumed 
to  be  correct.  The  reasons  for  these  changes  are  that  a  body  moving 
under  the  free  surface  of  a  liquid  produces  a  flow  pattern  which  in  turn 
affects  the  body  itself.  In  addition,  the  nature  of  the  surface  in  the 
form  of  wind  waves  influences  the  body.  The  magnitude  of  all  of  these 
effects  is  related  to  the  depth  of  the  body  below  the  surface  and  its 
velocity, 

When  the  body  is  near  the  surface,  the  flew  pattern  produced  in  be¬ 
tween  the  body  and  the  surface  is  a  Bernoulli  contour  system.  The  positive 
pressures  at  the  bcw  and  stem  produce  wave  crests,  while  the  lower 
pressure  at  mid-body  produces  a  truugn.  The  reduced  deptn  over  the  mid¬ 
portion  of  the  body  reduces  the  streamline  flow  area  between  body  and 
surface.  The  establishment  of  the  Bernoulli  contour  system  increases 
the  positive  pressures  at  bow  and  stem  and  decreases  the  mid-body  pressure 
on  the  upper  surfaco  of  the  boefcr.  Since  this  effect  is  not  produced  on 
the  bottom  at  relatively  low  velocities,  the  low  pressures  act  over  a 
larger  area  on  the  top  surface  resulting  in  a  net  vertical  force  or.  the 
body  acting  to  push  the  body  toward  the  surface. 

As  the  velocity  of  tie  body  increases,  the  wave  pattern  on  the  surface 
lags  more  and  more,  resulting  in  changes  to  the  pressure  pattern  to  the 
extent  that  the  positive  pressure  in  the  bow  area  covers  more  and  more  of 
the  upper  area.  A  not  transverse  force  in  the  downward  direction  is  then 
produced.  Because  of  the  Bernoulli  contour  system  effect,  the  positive 
pressure  at  the  bow  is  greater  while  the  positive  pressure  at  the  stem 
is  lover.  This  condition  results  in  an  increase  in  pressure  drag, 
varying  with  body  velocity  and  depth,  and  known  generally  as  wave  resistance. 

Che  of  the  references  indicated  that  the  change  in  damping  ratio, 
in  particular,  could  be  substantial .under  near-surface  conditions®,  but 
another  one  indicated  the  opposite10.  Therefore,  impulse  runs  made 
using  coefficients  modified  for  the  above  effect  and  were  compt  -a  with 
the  deep  coefficients. 

Coefficients  were  available  at  6  and  12  knots,  so  impulse  runs 
were  run  at  these  speeds.  Pitch  angle  was  recorded  as  a  typical 
parameter.  The  two  runs  for  deep  and  surface  coefficients  were  plotted 
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against  each  other.  A  shift  in  the  vertical  axis  was  due  to  different 
initial  conditions  of  level  flight,  and  the  result  shows  vexy  little  dif¬ 
ference  in  the  response  at  either  6  or  12  knots.  In  each  <wjse,  the  "ini¬ 
tial  conditions"  program  ECltfO  was  used  to  generate  a  neutral  pitch  angle, 
normal  velocity,  and  stemplane  deflection  for  level  flight. 

Because  the  critical  speed  for  the  demonstration  submarine  is  below 
six  knots,  it  was  decided  to  repeat  an  impulse  run  at  two  knots  to  ensure 
that,  performance  did  not  change  at  very  low  speeds.  Coefficients  were 
not  available  at  this  speed,  so  the  values  for  6  and  12  knots  were  ex¬ 
trapolated  into  the  2  knot  area.  The  main  changes  were  i nZW  and  MW 
which  were  lowered  by  3%  and  15?  respectively.  This  run  had  a  minor 
phase  shift  that  would  not  be  noticeable  to  the  trainee,  but  was  c  iherwise 
identical  to  the  deep  coefficient  case  at  two  knots. 

The  results  of  the  impulse  runs  indicate  that  the  handling  qualities 
of  the  demonstration  submarine  do  not  change  appreciatively  d”e  to  the 
presence  of  a  calm  near-surface.  A  further  check  was  made  using  overshoot 
runs  as  a  criteria.  The  change  in  forward  speed,  pitch  angle,  and  depth 
at  values  of  maximum  change  between  the  two  sets  of  coefficients  are  given 
in  table  12. 


TABLE  12.  .'3RC3HI  CHANGE  OF  CHANGE,  DEEP  TO  SURFACE  COEFFICIENTS 


6  Knots 

12  Knots 

AU 

7.7 

13.0 

A  6 

-1.1 

3.7 

AZ 

-ii.3 

10.9 

Because  the  change  in  u  is  very  small  during  thase  overshoot  runs, 
the  13?  difference  represents  an  absolute  change  of  less  than  one  tenth 
of  a  knot  at  twelve  knots.  A  comparison  of  the  curves  of  the  runs  show3 
that  very  little  change  occurs  when  surface  coefficients  are  used. 


The  one  thing  that  might  bother  the  trainee  is  that  the  conditions 
for  level  flight  are  different  close  to  the  surface.  For  example,  the 
neutral  pitch  angle  becomes  slightly  negative  as  the  surface  is  approached. 
It  was  also  felt  that  much  of  the  change  in  overshoot  response  was  due  to 
s  difference,  so  a  run  was  made  at  twelve  knots  and  used  surface 
coefficients  except  for  Z^*  ar.dM^.  The  changes  in  this  case  are  seen 
in  table  13. 


In  other  words,  most  ox'  the  change  during  an  overshoot  run  was  duo 
to  the  different  and  Zs|fr  determining  different  neutral  angles  ana 
not  due  to  the  use  of  other  surface  coefficients. 


Therefore,  it  will  not  be  necessary  to  u  ;e  surface  coefficient©  in 
e.ther  trainer.  If  very  accurate  near-surface  operation  is  required  as 
part  of  the  training  mission,  fyt^and  Zjfc  should  be  computed  as  a  function 
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PERCENT  ChASO?.  SUP.5 AGL'  jCEFFMEKTC 


Parameters  at  12  Knots 

1  %  Change  ! 

ALL 

3.6 

&& 

1.1 

AZ 

1.3 

oi  veloci ty  and  depth  in  the  modified  program. 


The  valuta  of  the  wave  coefficient  magnitudes  were  not  supplied  by 
NSRDC  aithov^!:  work  ha^  been  done  in  this  area*  An  investigation  was 
made  into  the  magnitude  of  these  coefficients,  during  a  Sea  State  $  sea 
using  the  classical  methods  of  hydrodynamics  . 


Calculations  for  the  wave  surface  amplitude  have*  been  calculated  in 
the  section  on  Wave  Generation.  The  forces  and  moments  are  resolved 
through  the  difference  between  the  heading  angle  of  the  submarine  and  the 
heading  angle  of  the  '~\ves.  Usually,  a  simplification  is  made  that  the 
wind  is  always  blowing  from  the  north,  but  a  variable  wind  direction  can 
be  added  easily.  The  forces  and  moments  are  then  varied  as  a  function 
of  the  depth,  reaching  a  negligible  value  at  about  five  hull  diameters 
below  the  surface. 


D.  MODIFIED  AND  BASIC  PROGRAM  SIMULATION 

The  modified  submarine  math  model  can  be  ^nltten  when  all  the 
approximations  to  th8  full-set  equations  of  motion  have  been  determined. 

The  first  step  is  to  set  all  negligible  coefficients  to  zero  and  remove 
the  parameters  not  being  used,  Th-  equations  are  solved  for  the  six 
accelerations,  and  any  set  of  terras  that  can  be  reduced  to  a  constant 
is  replaced  by  that  constant,  A  certain  amount  of  juggling  is  performed 
in  order  to  reduce  the  number  of  arithmetical  operations,  and  the  final 
result  is  known  as  the  coiqpact  submarine  simulation  program.  This  is  the 
program  that  would  be  used  in  a  submarine  simulator  irrespective  of  whether 
analog  or  digital  computation  was  used  and,  in  the  case  of  a  digital 
computer,  whether  or  not  FORTRAN  or  assembly  language  was  need.  The 
program  for  six-degrees-of -freedom  movement.  Compact  Equations  of  Motion 
(EC780),  are  given  in  equations  (12)  through  (17). 

Axial  (12) 


d~  (+  M  (Xq  (<£*+ -Zq  (pY  +  +  v  Y  ~  -*•  %) 

+  +Xte%*,Sl  +  X$bSb6l+  Z\) 

+  (b^  +  CC  uc)  -  (W- B)  S  IN  G)/(M  -  X  •  ) 
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Lateral 


(13) 


lioraal 


f\T  =  (  uCX^r 

*Y.f  +Y>f>  *  Sitin'*#')1 
4-  Vfc  vjf&y  -v  cos  8  sm 

(li*) 


*■  M  t£*(  *>**<£>  -  X^P'V 
+  +  +  *  *>«*flr* 

♦  l^UJ  (flT*  4.  *)* 

♦  IW-&)  cosecos<^)/(^"^) 


Roll 


(15) 


m,  -uip  ^urVfy  Kfr  *  ***• 

♦  u(Kvp  +  vc^at)  *•  AT(nrl*  *>%)k 

~  K*(2fcW  -£#B)  cos©*cos$ 


Pitch 


(16) 


1WII  . 

=  Hrfrp  4-tt.,a>  4-  a(j\%  * 

-  m  -rtvr  vva)^)  -  X^(vl>  -  ^ 

*  M^r*  4-  Mwr<*r  *  Hw<a«rl 

*  ♦*>')*  +  M««»u,^,'4U>t) 

4-  u1  {.M*  4-  Miw4%  +  Kttss) 

-  W^X^WcosGcos^  4  -  i.%B)  sm8) 
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Nr%P%  +  Nr#  +  •*•  «*'  Ntv6,. 

*  m  *.X4(<v  -uijp  ♦  urVN£  +  mx)x 

+  *4(Nfp  ♦  Krr  ♦  +  N^lX^W)  cot>0  vm  <► 


The  parameters  are  defined  in  equation  (18).  In  the  equations  that 
follow,  the  terms  on  the  left  side  of  the  equation  are  used  in  the  compact 
submarine  simulation  and  the  terms  on  the  right  side  are  the  original  NSRDC 
coefficients  as  defined  in  NSRDC  Report  2510.  The  only  exception  to  this 
is  Tjj,  Tcj,  etc.  which  are  used  on  both  sides  of  the  equation  once  the 
terms  for  the  inertial  mass  has  been  inverted.  These  values  are  different 
in  the  two  programs,  because  one  is  normalized  and  the  other  is  not,  but 
they  employ  the  same  symbology  because  they  occupy  the  same  locations  in 
the  two  equations  of  motion. 

Coefficient  Description  (EC780)  (18) 


X0£D(t 

s 

* 

£  (£y 
*  * 

X  OSDS' 

«■* 

= 

£  q}  x 

ZJL 

XOBDB* 

= 

a 

£  9}y 
*■  **  A*vw 

AH' 

2 

*'u 

a 

£.  0*a 
^  *  an 

Air 

2 

b‘n 

£  iTb 

AV3‘ 

a 

C4l 

a 

A  21* 

* 

aU 

= 

i 

•» 

bU 

- 

-  £*b 

Z  A* 

«• 

c’t. 

5 

£  l*c 

60 


NAVTRADSVCEN  68-C~O05O-X 


A^l' 

« * 

5 

- 

A*2‘ 

»  bU 

-  C  r  bl% 

A^‘ 

=  cl» 

- 

XUD' 

*  Xa' 

-  ^3Xa. 

YR* 

*  Y* 

r 

=  (|£*Y, -m) 

YRO’ 

-  V’ 

'  f 

»  4*% 

YPO1 

=  y'* 

*  irv. 

YP* 

*  y 

t 

YV 

*  Y' 

<r 

*  €  *‘Y„ 

YVAV1 

«*  Y*  , 

-  4**^ 

YOR* 

* 

•  f^Yw 

rvD 

*  Y* 

*Q’ 

"  £'l 

US3L 

-  4^ 

£RR’ 

=  2’ 
rr 

=  4 

KAVTHADEVCEH  68-0-00^0-1 


ZVR' 

*» 

Z‘ 

3 

i  !•  i„r 

ZSTR* 

tt 

« 

££lz* 

ZW* 

*• 

z’ 

*“44» 

3 

« J2l  Z„ 

ZWAW* 

V 

z' 

wM 

B 

- 1*? 

2  ^  E*wiu>i 

ZVV’ 

«* 

z*‘ 

CW 

5 

UX» 

EDS1 

e 

z’ 

Civ 

Z0&1 

« 

= 

iti* 

2W0 

«• 

ECr 

*3 

{&*£*■ 

3 

i/  u*  - 1  £ 

AKRO* 

at 

KV 

3 

1&X-T* 

AKP‘ 

St 

*V 

* 

£  £%  • T* 

AKv/O* 

«• 

«fr 

= 

•  T* 

AKV 

«o 

*C 

* 

AKVAV* 

* 

K* 

•'•rim  . 

- 

r^t-o 

f 

5 

f 

f* 

AK¥»o* 

«■» 

KV 

8 

62 


KAVTRADSVCBi  68-C-OOb'O-l 


T, 

i/a,  - 1 

AHRP* 

= 

- 

cx,  -I, +  f  J2.’m 

AMRR‘ 

= 

m' 

rir 

“ 

?  •  T, 

Amwd* 

Kb 

- 

i  *T» 

AHVR* 

S 

K.r 

= 

?  A*  •  T. 

AMQ' 

r 

{  jH*Mt  •  Ts 

AMAWQ* 

* 

= 

C  ?  T, 

AMSTR* 

m; 

= 

*  Ts 

AMW* 

- 

K. 

= 

£  £*  ft  *  j 

Z  *  ‘*o  *» 

Amwaw' 

= 

MU., 

r 

AMW’ 

- 

- 

l  f  M„  •  T. 

Amos’ 

- 

— 

1  £  Mu  -  T, 

AMOB* 

s 

«w 

s 

S4X.-T. 

AM<?0* 

= 

* 

c 

T, 

63 


NAVTRADEVCSN  68-C-0050-1 


anpq’ 

s 

N'n 

X 

AMP0‘ 

NV 

s 

AM  VO* 

•4 

a 

ANP‘ 

» 

"V 

s 

AMR* 

« 

»C 

B 

ANV1 

* 

m; 

* 

amvav' 

«  ‘ 

» 

ANDR* 

« 

•C 

«• 

ANRO* 

* 

Nf 

s 

s 

V. 

X 

i/a,  -  f  2’Nfi 

-:a*n,  *t4 

4  2  N*  •  Tk 

{ I'K  •  T„ 

.  C  ^^.Tk 

I  Ntr  *  T<. 

\ 

CB*> 


6k 


NAVTRADEVC15N  63-C-OOSO-l 


Since*  MRP  and  NPQ  can  be  added  to  the  cozqpact  simulation  at  no 
additional  cost,  these  terns  were  kept  even  though  their  necessity  is 
doubtful. 

Although  the  compact  coefficients  could  be  calculated  by  hand, 
it  is  far  easier  to  use  the  computer,  this  is  especially  true  consid¬ 
ering  that  the  input  deck  for  program  EB920  contains  all  of  the  necessary 
NSRDC  coefficients.  The  resulting  values  will  have  to  be  punched  on 
cards  as  input  to  the  compact  submarine  simulation  program. 

Program  EC790  -  Compact  Coefficient  Calculator  was  written  to  perform 
this  chore.  The  program  consists  of  a  main  program  and  Wo  subroutines 
called  INPUT  and  WRITE.  The  main  program  calls  INPUT  to  read  in  a  deck 
of  cards  and  then  calls  WRITE  to  write  them  out.  The  seventy  or  so  new 
coefficients  are  calculated  and  WRITE  is  called  again  to  write  out  the 
new  values.  These  new  values  are  then  punched  on  cards  in  a  format  t?  • 
can  be  used  as  an  input  to  the  compact  submarine  simulation  program. 

Subroutine  INPUT  was  taken  directly  from  the  main  program  of  EB92G. 

It  reads  in  the  same  data  deck  as  EB?20  including  the  values  that  are  not 
going  to  be  used.  This  saves  having  to  make  up  a  new  input  deck  to  use 
this  program  after  EB920  has  been  used  in  the  investigation  of  a  set  of 
coefficients. 

Subroutine  WRITE  was  also  taken  from  EB920.  It  uses  the  same  output 
format  for  printing  out  the  coefficients.  Whan  the  program  prints  out 
the  same  format  for  the  second  time,  the  values  have  become  those  on 
the  output  cards.  These  are  then  used  in  the  compact  simulation  pregram 
even  though  the  nomenclature  remains  the  same.  Values  not  used  by  the 
new  program  remain  the  same.  . 

Program  EC?80  -  Compact  Submarine  Simulation  is  a  FORTRAN  program 
using  the  mathematical  model  for  the  reduced  equations  of  motion.  It 
represents  the  chief  goal  of  the  study  in  that  it  is  as  small  as  possible 
while  still  simulating  the  path  of  a  submarine  accurately. 

This  program  can  be  run  in  any  digital  computer  that  has  a  FORTRAN 
compiler  and  enough  core.  It  has  bean  run  on  the  IBM  36O/I1O,  SDS  Sigma  5, 
and  an  SDS  Sigma  2  which  is  a  small  16-bit  machine.  In  every  case,  it 
gave  good  results  when  compared  against  EB?20  with  the  same  set  of  input 
conditions. 

The  MAIN  program  reads  in  the  input  cards  produced  by  EC790,  prints 
out  a  heading  called  UPDATE,  prints  out  the  twelve  parameters  of  velocity 
and  position  and  then  cheeks  the  tins  to  see  if  the  program  is  done. 

The  UPDATE  subroutine  is  equivalent  to  a  combined  EVAL  end  IN TEG 
subroutine  in  EB920.  It  calculates  velocities,  thrust,  trig  functions ; 
accelerations,  Siler  angles,  and  then  integrates  using  2nd  order  Adams. 
This  program  is  all  that  is  necessary  for  the  flight  portion  of  a  sub¬ 
marine  simulator. 
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Subroutine  CONTR  is  identical  to  the  sane  subroutine  in  EB920, 

It  Mas  added  so  that  the  compact  program  would  be  cheeked  out  against 
EB920  and  no  attempt' was  made  to  shorten  it*  It  allows  climbing  turns, 
flat  turns i  meanders ,  overshoots,  and  impulse  runs  in  the  same  manner  as 
the  original. 

Finally,  the  output  format  for  EC780  was  designed  so  that  the  twelve 
variables  could  be  printed  out  on  a  teletypewriter  rather  than  a  line 
printer.  This  is  so  that  a  very  small  computer  can  be  used  and  useful 
results  can  still  be  obtained. 

8.  WAVS  GBHSR1TICN 


The  analytical  treatment  of  wind-generated  ocean  waves  has  undergone 
a  revolution  in  the  last  five  years.  This  started  in  1953  when  Or,  Pierson 
and  Or,  Neumann  published  "Practical  Methods  for  Observing  and  forecasting 
Ocean  Waves  by  Means  of  Wave  Spectra  and  Statistics".  The  field  was 
entirely  shifted  over  from  a  determinate  one  to  a  statistical  basis  by 
1963  with  the  publication  of  "Ocean  Wave  Spectra"1?,  This  was  an  account 
of  proceedings  of  the  Fasten,  Md,  conference  in  that  year.  The  field  is 
a  complex  one  and,  rather  than  giving  a  superficial  treatment  here,  this 
section  will  attempt  to  outline  some  of  the  references  that  lead  to  an 
understanding  of  wave  spectra.  This  will  be  followed  lay  a  discussion  of 
two  different  methods  of  generating  wave  amplitude  with  a  digital  oomputer 
program. 

The  basic  reference  in  Chapter  8  of  "Wind  Wavoe*^  by  Blair  Kinsman, 
This  chapter  gives  a  short  history  of  the  development  of  wave  spectra 
theory,  discusses  some  of  the  stationary  Gaussian  process  mathematical 
models,  shows  the  development  of  the  Neumann  spectrum,  and  finally  indicates 
the  uses  to  which  the  Pierson  -  Neumann  theory  can  be  put.  Some  of  the 
mathematics  is  fairly  rough  going,  but  the  text  is  very  dear  and  under¬ 
standable. 

The  second  reference  that  is  a  "most*  is  "Sea  Spectra  Simplified"1^ 
by  Walker  H.  Michel,  This  report  discusses  the  concept  of  a  wave  spectrum 
from  a  physical  point  of  view  and  shows  why  a  frequency  spectrum  is  used, 
written  with  a  minimum  of  mathematics.  More  ijgjortantly,  it  gives  a 
review  of  the  different  spectra  that  have  been  developed  in  the  past  few 
years  and  outlines  the  notation  and  units  used  by  the  different  workers 
in  the  field.  Just  to  illustrate  one  source  of  confusion,  some  use  energy, 
some  use  wave  height,  and  some  use  amplitudes.  Finally,  ten  references 
are  given  that  oover  the  development  of  the  wave  speotra  theory, 

The  first  nine  pages  of  "Recent  Developments  in  Seakeeping  Research 
and  its  Application  to  Design"1?  give  a  brief  discussion  of  the  same 
areas  as  Michel's  paper.  The  end  of  the  paper  has  a  list  of  11*2  references 
covering  almost  every  paper  of  the  "statistical  era"  since  I960.  It  is 
available  as  a  reprint  from  SRAMS. 

1 

The  three  references  listed  give  a  good  overall  view,  but  they  do 
not  tie  down  an  actual  spectrum  that  could  be  used  throng  a  wave  generation 
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investigation*  This  is  done  in  "A  See  Spectrum  for  Model  Tests  and 
Long-Term  Ship  Prediction"16  by  J.R.  Scott.  A  particular  spectrua  is 
given  and  justified  along  with  the  transformation  that  determines  the 
energy  in  the  spectrum  proportional  to  either  the  average  wave  height 
(sea  state)  or  the  wind  speed  existing  at  the  place  where  the  waves  are 
assumed  to  exist. 

All  of  these  references  indicate  how  to  measure  waves,  determine 
the  spectrum  characteristics,  and  apply  the  theoretical  results  obtained! 
but  they  do  not  give  any  help  as  to  how  a  mathematical  model  might  be 
developed  to  generate  wave  amplitudes  with  a  computer.  The  two  different 
approaches  are  covered  in  "The  Analysis  and  Modeling  of  Irregular  Waves*1? 
and  "Mathematical  Generation  of  a  Realistic  Seanl8  by  Chen  et  al.  The 
second  report  is  out  of  print,  but  it  may  be  obtained  from  the  Defense 
Documentation  Center  under  number  AD609906  by  qualified  users.  A  review 
of  these  papers  and  references  will  give  a  complete  background  into  the 
subject  of  statistical  wave  generation  for  research  purposes. 

However,  no  one  in  the  field  seems  to  be  using  these  techniques 
with  a  digital  computer  to  produce  waves  with  the  desired  statistical 
properties  in  real-time.  The  closest  approach  is  in  the  wave-making 
machine  used  with  a  towing  tank  for  testing  uses.  In  this  case,  signals 
are  generated  using  an  analog  computer  and  recorded  in  analog  form  on  mag¬ 
netic  tape.  The  tape  is  played  back  when  a  wave  signal  is  desired.  The 
two  different  methods  were  examined,  a  mathematical  model  was  developed, 
and  a  computer  program  was  written  to  generate  a  random  wave  over  time. 

The  concept  of  a  wave  spectrum  is  based  on  the  energy  content  of 
the  wave,  which  is  proportional  to  the  sea  state,  being  divided  into 
different  frequencies  of  superimposed  sinus odial  waves.  The  energy 
divides  according  to  hydrodynamic  principles  into  a  particular  frequency 
curve  which  can  be  determined  by  making  many  measurements  of  waves  and 
analyzing  the  resulting  data.  The  irregular  surface  is  due  to  the  fact 
that  these  sinusoids  have  phase  angles  which  occur  in  a  completely  random 
manner  with  time.  Therefore,  they  add  and  subtract  to  present  the  typical 
irregular  surface  seen  on  the  ocean,  rather  than  the  regular  wave-form 
expected  from  a  trigonometric  Fourier  expansion  with  fixed  phases. 

The  energy  in  a  sine  wave  is  related  to  the  amplitude  by  the  expression 


where  p  is  the  water  density  and  h  is  the  wave  height. 

To  calculate  a  wave  height  at  a  particular  point,  a  sine  wavo  is 
calculated  for  each  of  a  number  of  frequencies.  Their  amplitude,  and 
therefore  iheir  energy,  at  each  frequency  is  determined  by  the  wave  spectra 
curve.  These  frequencies  are  shifted  in  random  phase  and  added  together 
for  each  time  point.  The  main  difference  between  the  two  methods  outlined 
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is  how  the  frequency  division  is  made.  The  older  of  the  two  methods  is 
the  addition  of  a  finite  number  of  cosine  waves  of  different  amplitudes 
and  frequencies.  The  following  mathematical  model  was  developed  from 
"Mathematical  Generation  of  a  Realistic  Sea"  and  section  8,3  of  "Wind 
Waves".  The  surface  elevation  can  be  represented  as  a  stationary  Gaus’sian 
process  In  three  dimensions.  This  can  be  reduced  to  a  function  of  time 
by  assuming  that  only  the  single  point  directly  over  the  submarine  is 
of  interest,  and  placing  the  coordinate  reference  system  origin  at  this 
point.  The  surface  elevation  then  becomes 

90  <jj> 

i?(t)  =X  4  CoS  +  €(u/lP)]  \/A2(^t/S) 


where 


u*>  s  frequency 


s  difference  between  wind  direction  and  wave 
direction 

€  (iUj  >3  )  a  random  phase  angle  (0-2W ) 

:  wave  spectrum 

The  reason  for  putting  the  differential  under  the  radical  sign  is 
covered  in  Chapter  8  of  "Wind  Waves".  In  order  to  program  this  integral 
on  a  digital  computer,  it  was  necessary  to  divide  it  into  discrete  elements 
so  that  they  c an  be  added  numerically.  This  is  a  double  integral  so  the 
division  has  to  be  made  in  both  w  and  ft  . 

The  first  task  is  to  divide  the  directional  dependence  energy  spectrum 
into  equal  bands  on  each  side  of  the  primary  wind  direction.  This  is 
usually  assumed  to  be 


o 


%  -  % 

oihcvvw  »%e 


This  converts  the  spectrum  into  a  function  of  frequency  only  and 
allows  separate  calculations  of  the/9  part  of  the  integral.  We  then 
can  divide  the  integral  according  to  the  example. 


X  cotZ/&  4  & 
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This  now  has  five  equal  energy  bands  off  the  main  wind  axio  by 
using  constants  evaluated  iri  table  lit  as  calculated  in  Chen's  report. 


TABLE  lit.  WAVE  ENERGY  BANDS 


A/6  SECTOR 

at  =f<os*j0  4/i\siCTOR 

-60° 

-90°  to  -52° 

0,0906 

-30° 

-32°  to  -l5.ii° 

0.1030 

0° 

-l3.1i°  to  +15.I10 

0.3325 

30° 

l5.ii°  to  32° 

0.1030 

60° 

32°  to  90° 

0.0906 

3L  (,  represents  the  energy  in  a  sector  which  contains  >5  dividing 
the  seotor  into  two  equal  areas  on  each  side  of  JB  . 

The  wave  spectrum  is  also  divided  into  equal  energy  areas.  The 
total  area  under  the  spectrum  curve  represents  total  energy  in  the  wave 
and  it  is  a  single  Job  to  solve  for  the  frequency  points  which  give  equal 
areas  after  integrating  the  function  for  A2(w)  .  If  each  of  these  points 
istu^  |  then 

u/j-w'cj-0 

i  z 
Awi =  ~ 


We  can  now  write  the  surface  elevation  as  a  summation  of  discrete 
points 


S  N 

c°s Cw t  +  £; j } Awj 


where  all  the  terms  have  been  defined. 

This  mathematical  model  was  programed  using  the  Pierson  -  Moskowitz 
spectrum.  This  is  a  spectrum  curve  that  is  being  used  as  the  Interim 
standard  in  both  analytical  and  experimental  studies.  Its  equation  is 
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* 


8.1  x  ICT^g^ 

UJ  5 


4 


where  U  is  the  wind  speed*  Hie  actual  equation  programed  was 


XICf*x  32»Zg> 

ir 


cos(W^  t+ 


W  |  had  ten  different  values  which  were  computed  from  the  spectrum 
and  ran  from  a  lower  limit  of  0.2  rad/soc.  to  1.8  rad/sec.  U  was  choson 
to  be  UO  feet  per  second  and  the  *s  were  read  in  from  cards.  A  random 
number  generator  was  called  fifty  times  and  scaled  for  0  to  2TT  to 
calculate  values  for ?  .  These  values  were  computed  and  stored  each 
time  the  program  was  used.  Figure  8  shows  iiOO  points  calculated  at  one 
second  intervals.  The  T  axis  is  in  per  cent  of  the  largest  point  calculated. 
This  would  correspond  to  a  wave  height  of  25  ft,  at  this  wind  speed. 

This  program  takes  about  one-tenth  of  a  second  to  calculate  one 
point  in  time  and  requires  quite  a  bit  of  core  space  to  hold  the  values 
used  to  make  up  the  auiaaation  (fifty  in  the  case  of  the  example),  so 
it  is  useful  only  for  the  research  program.  It  would  hava  to  be  greatly 
simplified  for  use  in  a  submarine  simulator. 

The  second  method  was  recently  borrowed  from  communications  and 
filter  theory.  It  consists  of  a  filter  network  that  has  a  frequency 
response  that  is  equivalent  to  the  desired  spectral  density  energized 
by  a  "white  noise*  random  excitation.  "White  noiae*  has  a  constant  spectral 
density  containing  all  frequencies  equally,  and  s:'noe  a  filter  in  the 
frequency  domain  can  be  considered  as  a  variable  gain  device,  the  output 
will  have  frequencies  statistically  equal  to  that  of  ocean  waves.  White 
noiae  cannot  be  realised  as  a  practical  matter,  but  ocean  wave  spectra 
are  limited  to  2  radians  per  second,  so  a  Hat  band-limited  signal  will 
do.  This  is  a  very  easy  process  when  using  an  analog  computer,  A  gas 
tube  is  used  as  a  white  noise  source  and  a  filter  is  designed  and  programed 
using  conventional  techniques.  The  job  is  a  little  harder  when  using  a 
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digital  computer. 


figure  8.  Typical  Wave  Output, 

Hia  "Analysis  and  Modeling  of  Irregular  Waves"  gives  a  mathematical 
model  that  can  be  used  as  a  filter  to  create  the  Pierson  -  Moskowits 
spectrum.  The  desired  frequency  response  for  the  filter  is 


H  (u,)  =  (a2  (w))z  = 


-.57* 


and  the  frequency  response  can  be  duplicated  by  means  of  a  low  pass 
filter  followed  bv  a  three  section  high  pass  filter.  The  constants 
were  derived  empirically  by  plotting  H(w)  and  comparing  the  curve  with 
the  filter  frequency  response  as  the  damping  factor  and  natural  frequency 
of  each  section  va3  varied.  The  filter  equation  is 


5  +2<fuWhuS+  \SM2SHWhHS+W 


* 
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TABLE  15.  WIND  CONSTANTS 


msss 

Sl 

J&kL. 

u  5 

0.39 

0,325 

0.21 

0.70? 

iiO 

0.i*3 

0.25 

0.21* 

0.707 

35 

0.1*8 

0.325 

0.27 

0.707 

The  factors  were  given  in  the  reference  and  were  used  in  a  digital 
computer  program  to  generate  waves  for  the  hO  knot  case*  This  program 
was  not  included  in  the  study  library  of  programs  because  it  uses  a 
noise  generator  that  is  both  IBM  36O/I1O  machine-dependent  and  GAC  DOS 
Operating  System  dependent,  and.  therefore  will  run  only  at  t.he  Goodyear 
Aerospace  facility.  A  white  noise  generator  that  operates  as  fast  as 
possible  and  has  a  Gaussian  probability  distribution  function  would 
be  written  far  the  particular  computer  used. 

The  programing  of  the  filter  is  straight  forward.  It  is  necessary 
to  find  a  set  of  difference  equations  having  a  system  function  R(z)  which 
is  equivalent  to  that  of  H(s).  A  method  due  to  Baxter”  is  given. 

If  a  transfer  function  is  written  in  terms  of  s,  the  Laplace  operator, 
the  function  1/s  can  be  considered  a3  an  integration.  Therefore,  if  an 
integration  z- transform  I(z)  is  known,  the  filter  z  function  can  be  written 
by  replacing  the  s  with  the  factor  l/l(z).  An  example  will  make  this 
clear.  Assume 


Filter 


H(S)  = 


sx  +  £Wh<fS  +  Wh* 

Euler  Integration  with  time  integral  h 

T/  \  -  hZ-1 

1(2)  -  rrpt 

Thus  H(z)  equals  B(l/I(z))  which  equals 
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or 


hlz-z 


1-(Z+  2Wh  S  h)  z=r^-  (1-2  W),  i h+ W(,1h*)Z'*' 


This  can  be  written  directly  as  a  difference  equation  through  the 
method  outlined  in  the  section  on  integration  aetheds  end  can  then  be 
programed  into 


1  „  =  (  % +  2  Wh  s  h)  V  J  -  (»-* 5  h  +WhV)  V*+  H’ij 


The  same  method  is  used  with  the  other  sections  of  the  filter. 

The  first  HI  h-3.  is  the  output  of  the  random  noise  generator  and  each 
H n  becomes  the  tjj  of  the  next  section. 

This  method  worked  as  well  as  the  cosine  series  addition  because 
there  were  inly  16  multiplication  and  additions  along  with  the  generation 
of  the  random  number.  It  should  be  noted  that  most  Gaussian  random-number 
generators  are  extremely  slow  and  are  not  suitable  for  real-time.  The 
one  used  at  GAC  is  written  in  IBM  machine  language  in  order  to  perform 
shifting  operations  that  are  difficult  to  do  in  FORTRAN.  This  is  not  a 
problem  with  the  cosine  generator  because  the  xmndom  numbers  are  generated 
only  when  the  subroutine  is  first  entered. 

Since  this  work  was  done,  another  reference  on  direct  programing 
of  Laplace  transfer  functions  was  found.  The  reference  "Digital  Filter 
Design  Techniques  in  the  Frequency  Dona  in"  by  Rader  and  Gold  gives  several 
methods.  One  that  has  been  utilized  for  filter  design  on  another  project 
at  Goodyear  Aerospace  is  as  follows.  Replace  the  natural  frequency  of  the 
filter  by  a  new  natural  frequency 


Wj.Cn€w#)  =  tsn 


w».h 


where  h  is  the  time  interval  of  the  difference  equation  to  be  used. 

Then  replace  s  by  (z-l)/(z+l)  in  E(s)  and  perform  the  algebra  necessary 
to  express  the  resulting  H(z)  as  a  ratio  of  polynomials.  The  difference 
equation  y  can  then  be  written  diroctly. 
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F.  SMALL  COMPUTER  UTILIZATION 


The  utilization  of  a  small  general  purpose  computer  was  investigated 
onee  the  characteristics  of  the  compact  submarine  simulation  program 
were  known.  This  was  accomplished  in  the  most  direct  manner  possible) 
namely,  by  running  the  compact  program  in  a  16-bit  machine  and  examining 
the  results. 

There  is  no  doubt  that  a  l6»bit  computer  can  be  used  in  a  submarine 
simulator.  Accuracy  is  not  a  problem  if  a  higher  level  language  such  as 
FORTRAN  is  used  because  of  the  way  a  l6-bit  machine  utilizes  floating 
point  numbers.  Two  words  are  used  to  store  data.  The  first  word  contains 
the  exponent  and  the  first  six  bits  of  the  fraction.  The  second  word 
contains  the  rest  of  the  fraction  so  that  a  2l*-bit  word  plus  exponent 
is  actually  being  used.  This  allows  ssven  figure  decimal  acouracy  to 
be  obtained  in  the  simulation  computations. 

The  computer  used  for  running  the  compact  sim  ulation  program  was 
an  SB8  Sigma  2.  It  has  a  2.25  v^aec  add  and  cycla  time  which  is  typical 
of  this  class  of  computers.  It  did  not  have  either  floating  point  hardware 
or  fixed  point  railtiply-divida  hardware,  so  any  programs  represent  a 
worst  case  condition  of  slowest  speed  and  maximum  core  size. 

The  cpmpact  program  was  run  on  the  IBM  360/UO  (32-bit  word)  to 
produce  a  check  program  against  the  Sigma  2.  The  output  format  in  each 
case  was  six  decimal  figures  plus  exponent, and  the  two  outputs  agreed 
in  the  sixth  place  in  all  twelve  outputs  for  every  run  made.  Therefore, 
the  accuracy  is  exactly  the  same  whether  a  16-bit  or  32-bit  computer  is 
used  if  floating  point  numbers  are  used. 

In  the  past,  the  main  problem  in  using  floating  point  arithmetic 
has  been  trust  the  subroutines  are  slow  and  floating  point  hardware  doubles 
the  cost  of  the  computer.  Therefore,  the  timing  of  the  coapaet  program 
was  measured  under  worst  case  condition  to  see  if  real-time  computation 
was  a  possibility.  The  compact  program  was  run  through  eighty  complete 
cycles  in  slightly  under  fifteen  seconds.  This  represents  1*0  seconds  of 
real-time  whan  the  integration  interval  is  one-half  second  and  60  seconds 
wham  &  is  one  second.  The  duty  cycle  is  38  JS  in  the  first  case  and  19% 
in  the  second.  This  allows  plenty  of  time  for  additional  system  simulation, 
wave  generation,  and  input/  output  subroutines. 

The  core  size  presented  a  problem.  The  basic  core  size  for  a  small 
computer  is  LO96  words  and  any  above  this  amount  requires  the  core  to 
be  extended  to  8192  words,  thus  increasing  the  cost.  An  analysis  of  the 
core  used  for  the  compact  program  is  given  in  table  16. 

The  simulation  program  needs  2632  words  without  I/O  and  system 
simulation.  The  FORTRAN  I/O  routines  were  were  used  as  an  Indication  of 
what  was  needed  ard  the  total  was  just  over  the  1*096  limit  without  in¬ 
cluding  wave  generation  and  tank  simulation.  It  seems  unlikely  that 
enough  cutting  could  be  done  to  get  the  entire  program  under  this  amount, 
and  therefore  an  8192  core  is  indicated. 
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TABLE  16.  SUBROUTINE  CORE  SIZE 


Subroutine 

Size 

Total 

UPDATE 

1253  words 

SQRT,  SIN/COS 

317 

FLOATING  POINT 

66h 

DATA  TRANSFER 

130 

CGJMON 

268 

2532 

2632 

MAIN 

hl7  words 

FORTRAN  I/O 

1698 

ERROR  CHECK 

215 

2330 

2330 

1I962  words 

Hand  coding  might  help,  but  then  all  of  the  advantages  of  making 
changes  easily,  not  having  to  worry  about  scaling,  aid  being  able  to 
work  with  the  simulation  program  without  having  to  learn  a  particular 
language  would  be  lost. 

Fixed  point  hardware  irultiply-divide  might  help  the  problem.  The 
compact  program  was  run  on  a  SDS  Sigma  5  with  this  option  and  the  core 
space  was  cut  in  half.  Hardware  V/to  is  about  one- third  as  expensive  as 
extending  the  core. 

Floating  point  hardware  doubles  the  cost  of  the  computer  and  is 
used  only  when  speed  requirements  make  it  necessary.  It  cannot  be  justified 
at  present,  but  the  cost  keeps  coming  down  as  more  computers  came  onto 
the  market. 

A  further  cost  is  a  real-time  clock  and  interrupt.  This  is  necessary 
to  tie  the  integration  interval  to  the  outside  world  through  the  I/O 
system.  A  teletypewriter  is  also  needed  to  communicate  with  the  computer. 
The  total  costs  are  estimated  below. 


Basic  Computer  100$ 

Teletypewriter  PP/PR  reader  30$ 

Interrupt  6% 

Real-time  Clock  8$ 

Simulator  I/O  100$ 

WO  Hardware  20$ 

or 

Extended  Memory  60$ 

Therefore,  the  total  hardware  cost  of  a  submarine  simulator,  ex¬ 


cluding  the  motion  platform,  will  run  between  265$  and  325$  of  the  basic 
computer. 

The  I/O  equipment  doubles  the  cost  of  the  computer.  This  section 
contains  the  D-A  and  A-D  converter  discrete  inputs  and  circuits  that 
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tie  the  computer  to  the  motion  platform.  Generally,  it  is  designed  for 
the  particular  job  and  this  is  the  reason  for  the  expanse. 

there  are  two  ways  in  which  this  cost  can  be  reduced.  One  is  to 
do  all  programing  at  a  central  location.  A  computer  is  chosen  whose 
manufacturer  has  software  that  simulates  the  program  assembly  on  a  much 
larger  computer.  A  program  is  written  in  FORTRAN  and  run  on  a  Sigma  7, 
for  example.  A  paper  tape  is  produced  that  contains  an  object  program 
for  the  small  machine.  This  tape  is  sent  into  the  field  and  loaded  with 
a  simple  tape  reader  at  the  user's  facility.  A  good  program  will  write 
the  most  compact  program  possible  using  the  advantages  of  the  large  com¬ 
puter.  This  saves  the  cost  of  the  teletypewriter  for  each  small  computer. 
More  importantly,  it  reduces  the  cost  of  human  programing* 

The  second  step  uses  a  read-only  memory  in  the  small  machine.  The 
program  is  written,  debugged,  and  tested  at  a  central  facility.  A  tape 
containing  the  complete  program  is  sent  to  the  computer  manufacturer  who 
returns  read-only  memories  containing  this  program  for  each  simulator. 

These  are  plugged  into  the  small  computer  at  each  location  and  no  programing 
or  testing  is  required  at  all.  This  could  be  the  most  economical  method 
of  all  If  a  number  of  simulators  were  going  to  be  produced. 

A  review  was  made  of  the  small  computers  now  on  the  market.  Up 
to  a  year  ago,  the  only  computers  on  the  market  were  16-bit  machines  in 
the  class  of  the  SEL810A  (used  in  the  20A62  Emergency  Ship  Handling 
Trainer),  the  PDP-8  and  DDP-516 .  These  machines  ran  from  $18,000  up. 

Since  then,  a  number  of  different  manufacturers  have  cone  onto  the  scene 
and  prices  are  dropping  while  speed  is  increasing.  The  list  below  includes 
same  of  the  computers  that  should  be  considered  and  represents  the  state- 
of-the-art  at  the  time  of  this  report.  It  cannot  be  complete  because  of 
the  rate  at  which  the  field  is  changing.  Price  includes  4K  of  memory 
and  a  teletypewriter. 


TABLE  17.  SMALL  COMPUTER  BATA 


Computer 

Manufacturer 

Price 

Data  6201 
DataMate  8 

HP  yn)|i 
Interdot*  3 

NOVA 

MAC  16' 

VarUn  Bata  Machines 

Gamo  Industries 
Hewlett-Packard 
interdata 

Bata  General 

Lockheed  Electronics 

*14,500 

*15,000* 

*12,000 

*11,200 

*10,000 

*12,000 

♦price  includes  8£  core 


It  would  be  advisable  to  run  the  compact  simulation  program  on  a 
particular  computer  before  the  final  decision  is  made. 

G.  VERIFICATION 

Section  C-l  of  Specification  3kl-101,  Study  of  Submarine  System 
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Equations,  requires  that  certain  tests  be  made  to  verify  that  programs 
having  reduced  sets  of  coefficients  are  compatible  with  the  full  set 
of  coefficients  as  determined  by  David  Taylor  Model  Basin.  These  tests 
shall  be  computer  runs  simulating  submarine  operation  through  a  varied  set 
of  controlled  maneuvers  such  as  meander,  overshoot,  etc.  Tests  will  be 
run  for  the  full  set  of  coefficients  and  each  different  set  of  reduced 
coefficients. 

Three  groups  of  runs  were  performed;  one  for  the  full  set  of  co¬ 
efficients,  one  for  the  modified  set,  and  one  for  the  basic  set.  The 
modified  set  had  twenty-one  coefficients  from  the  full  set  to  zero,  and 
two  other  coefficients  changed  to  help  moderate  the  effects  of  eliminating 
coefficients.  The  basic  set  has  twenty-four  coefficients  set  to  zero 
and  four  other  coefficients  changed  to  help  offset  the  effects  of  elim¬ 
inating  coefficients. 

The  submarine  maneuvers  that  were  performed  by  computer  to  verify 
the  compatibility  of  the  reduced  sets  are  described  below. 

1.  MEANDER 

This  manuever  demonstrates  the  capability  of  a  submarine  to  recover 
from  a  rise  or  a  dive  without  use  of  primary  control  surfaces  or  emergency 
ballast  changes.  It  provides  an  evaluation  of  the  longitudinal  stability 
of  the  submarine. 

The  submarine  is  operated  at  a  given  speed  on  a  given  course  until 
steady,  level  flight  is  achieved.  The  submarine  is  then  disturbed  from 
its  neutral  pitch  angle  by  deflecting  the  sternplanes  a  fixed  amount  in 
either  rise  or  dive.  After  the  submarine  reaches  a  preselected  pitch 
angle  called  execute  pitch  angle  (®* ),  the  planes  are  returned  to  their 
original,  neutral  position.  They  are  held  fjxad  while  the  submarine  is 
allowed  to  pull  out  and  resume  a  condition  of  steady,  level  flight.  The 
resulting  pitch  angle  trajectory  provides  both  visual  and  quantitative 
evaluations  of  the  longitudinal  stability  of  the  submarine. 

Meander  tests  were  run  on  the  digital  computer,  program  EB920, 
Submarine  Simulation  Program,  for  each  set  of  coefficients  at  five  differ¬ 
ent  command  speeds  (!>,  10,  15,  20,  and  2 £  knots).  Maximum  diving  plane 
angle  was  20  degrees,  and  the  execute  or  pull  out  angle  was  10  degrees 
on  all  meander  maneuvers. 

Computer-plotted  graphs  of  stemplane  angle  versus  time,  pitch 
angle  versus  time,  and  depth  versus  time  were  recorded  for  each  run. 

The  damping  factor  was  calculated  for  the  full  set  data  and  modified 
set  data  in  accordance  with  methods  detailed  in  BuAer  report,  "Dynamics 
of  the  Airframe". 

The  curves  for  the  basic  data  showed  greater  dancing  and  could  not 
be  calculated  by  the  same  method.  Therefore,  the  damping  ratio  was 
estimated  by  comparison  of  these  curves  to  other  sets  and  prior  general 
knowledge  of  the  relation  of  damping  factor  to  response  curves. 
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2,  VERTICAL  OVERSHOOT 

The  overshoot  is  a  specific  maneuver  utilized  to  provide  numerical 
values  of  the  inherent  effectiveness  of  the  diving  planes  in  initiating 
and  checking  changes  in  depth. 

In  this  maneuver,  the  submarine  is  operated  at  a  fixed  command 
speed  and  course  throughout.  The  controls  are  maintained  for  level 
flight  until  the  selected  speed  is  achieved.  The  stern  elevators ' are 
then  deflected  to  a  given  angle  and  are  held  until  the  submarine  reaches 
preselected  pitch  angle,  called  execute  pitch  angle  (O')*  At  that  time, 
the  elevators  are  reversed  an  equal  and  opposite  amount  from  the  original 
neutral  position  and  held  until  the  rate  of  change  of  depth  reverses. 

This  maneuver  was  performed  on  the  digital  computer  for  the  throe 
sets  of  coefficients,  at  five  different  command  speeds  (5,  10,  15,  20, 
and  25  knots).  Maximum  diving  plane  angles  of  10  and  20  degrees  were 
also  used  for  two  complete  groups  of  curves.  Execute  pitch  angles  of 
hf  8,  12,  and  16  degrees  were  also  selected. 

Computer-plotted  graphs  of  stern  elevator  position,  submarine 
pitch  angle,  and  depth  versus  time  were  prepared  for  all  runs.  Units 
given  are  radians,  feet,  and  time. 

From  the  data  of  these  rune,  a  number  of  graphs  of  parametric  curves 
were  prepared  and  machine-plotted  which  provide  relative  measures  for  the 
runs  using  the  three  different  sets  of  coefficients, 

3.  STEAM,  SUBMERGED  TURN 

This  is  a  test  to  investigate  the  submarine  response  to  action  of 
its  rudder  in  horizontal  plane  maneuvers;  in  particular,  loss  of  speed, 
steady  turning  diameter,  roll  angles,  and  time  to  roach  90-  and  180- 
degree  change  of  heading. 

In  this  maneuver,  the  submarine  is  operated  at  a  fixed  command  speed 
throughout.  The  computer  program  provides  automatic  pilot  operation 
during  this  maneuver  to  maintain  depth  as  constant  as  pdssible.  The 
controls  are  maintained  fixed  until  steady,  level  flight  at  the  selected 
speed  is  achieved.  The  rudder  is  then  deflected  to  a  given  angle  and 
held  steady.  During  the  maneuver,  computer-plotted  graphs  were  made  of 
rudder  angle  (DR),  stemplane  angle  (DS),  pitch  angle  (0),  roll  angle (0), 
heading  angle  (\|0,  normal  velocity  (w),  yaw  angular  velocity  (r),  forward 
speed  (u),  and  depth  (z).  The  angle  of  attack  (oc  )  was  described  by 
this  equation  for  small  values  of  w. 


°<  =tavr'  ^  *  —  (foT  «  <io°) 
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Therefore,  the  plotted  values  of  w  may  be  used  for  values  of  o<  . 

The  data  from  these  runs  was  used  to  prepare  graphs  of  parametric 
curves  which  provide  relative  measures  for  the  runs  of  full  set  coeffi¬ 
cients  and  runs  of  modified  set  coefficients.  Results  of  runs  with 
basic  set  coefficients  were  almost  identical  to  runs  with  the  modified 
coefficients, 

ll.  RPM  VERSUS  SPEED 

Since  the  thrust  coefficients  have  been  neither  changed  in  value 
nor  eliminated,  there  was  no  need  to  repeat  these  calculations « 

<.  acceleration/deceleration 

This  test  determines  the  submarine  response  to  change  in  command 
speed  values. 

The  autopilot  function  is  operating  throughout  this  test  in  attempting 
to  hold  a  constant  depth.  The  submarine  is  started  from  an  initial  at 
rest  or  zero  velocity  condition  at  800  feet  depth.  The  command  speed  is 
increased  in  5  knot  increments  up  to  a  maximum  of  25  knots  at  equal  time 
intervals  of  230  seconds.  The  command  speed  is  then  reduced  in  5  knot 
increments  to  zero  over  the  same  time  intervals. 

Runs  were  made  for  full,  modified,  and  basic  sets  of  coefficients, 

6,  MAXIMUM  ACCELERATION/DECELERATION 

This  test  is  an  examination  of  the  submarine  response  to  a  large 
change  in  command  speed  or  throttle  setting. 

The  autopilot  is  in  operation  during  this  complete  maneuver.  The 
submarine  is  maintained  in  steady,  level  position  at  800  feet  depth. 

The  command  spaed  is  then  immediately  increased  to  25  knots.  The  value 
is  held  for  5 00  seconds  and  then  abruptly  reduced  to  zero. 

During  this  maneuver,  computer-plotted  graphs  of  the  same  parameters 
as  used  in  Section  5  were  made  versus  time,  Sg,  0,  u,  v,  q,  and  z. 

Runs  were  made  for  full,  modified,  and  basic  sate  of  coefficients, 

7,  LONGITUDINAL  TRIM 

This  test  shows  the  submarine  response  to  a  small  change  in  ele¬ 
vator  deflection. 

The  submarine  is  held  at  level  flight  at  a  particular  speed  and 
depth,  the  s  tempi  an  e  is  then  deflected  a  small  amount,  and  the  changes 
in  certain  dynamic  parameters  are  noted.  A  charge  of  3  degrees  was  used 
in  this  case.  These  runs  were  performed  at  5,  10,  15,  20,  and  25  knots 
for  all  three  sets  of  coefficients  over  a  time  period  of  2b0  seconds. 
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SECTION  IV 
RESULTS 


A,  RESEARCH  PROGRAM  RESULTS 


The  research  program  simulates  the  equations  of  motion  as  expressed 
in  NSRDC  Report  2f>lCr.  The  use  of  this  program  was  described  in  the 
section  on  "Methods  and  Procedures"  of  this  report,  A  brief  description 
follows, 

EB920  -  Submarine  Simulation  -  This  program  calculated  the  dynamic 
changes  in  submarine  position,  Telocity,  and  attitude  as  a  function  of 
time  in  six-degrees-of-freedom.  It  consists  of  the  following  parts. 

Main  Program  -  This  program  ties  the  subroutines  together,  calcu¬ 
lates  a  few  constants,  initializes  all  arrays,  writes  out  the  hydrodynamic 
coefficients  used,  writes  out  the  twelve  rates  and  positions  as  the 
program  proceeds,  writes  the  plotting  tape  if  plot  subroutines  are  avail¬ 
able,  and  finally  checks  to  see  if  multiple  runs  are  required. 

Subroutine  INPUT-  -  Reads  in  a  deck  of  cards  with  all  coefficients 
and  initial  values,  stores  values  far  initializing  during  multiple  runs. 

Subroutine  KUTTA  -  Integrates  twelve  variables  with  a  Runge-Kutta 
It-step  integration  subroutine. 

Subroutine  INTEG  -  Does  the  same  job  as  KUTTA,  but  has  a  ehoioe  of 
three  simpler  integration  routines. 

Subroutine  EVAL  -  This  routine  simulates  the  actual  equations  of 
motion  for  integration  with  the  KUTTA  subroutine. 

Subroutine  EVA L  1  -  Sams  as  EVAL  except  uses  subroutine  INTEG. 

Subroutine  CONTR  -  This  subroutine  allows  the  researcher  to  pro¬ 
gram  an  input  which  follows  standardized  submarine  tests  for  climbing 
turns,  flat  turns,  meanders,  overshoots,  and  impulse  inputs. 

Subroutine  INVER2  and  Subroutine  MATHPI  -  Theee  subroutines  invert 
a  matrix  and  multiply  it  by  a  column  vector  for  solving  the  equations 
of  motion  whan  extreme  accuracy  is  desired, 

ZC790  •  Submarine  Simulation,  Longitudinal  Freedom  -  This  program 
is  identical  to  EE02O  except  that  it  has  only  three- degrees-of -freedom. 
Therefore,  only  meander,  overshoot,  and  similar  runs  can  be  calculated, 
^hla  program  was  written  to  save  computer  time  because  it  ran  in  one- 
jnird  the  time  of  EB?20.  This  will  be  a  real  advantage  where  a  large 
number  of  longitudinal  runs  are  required.  It  is  limited  to  2nd  order 
Adam's  integration,  so  less  core  space  is  required. 
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A  number  of  additional  programs  were  written  to  assist  the  researcher 
in  submarine  simulation, 

£Ci*30  -  This  program  uses  the  location  and  weight  of  all  tanks  cm 
the  submarine  to  calculate  the  center  of  gravity  shift  limits  and  weight 
changes  for  any  ballast  arrangement, 

EC&70  -  This  program  calculates  the  neutral  angles  for  any  spaed 
in  level  flight.  The  coefficients  are  used  as  an  input  and  the  level 
flight  values  of  plane  angles,  pitch,  and  normal  velocity  are  calculated. 

ZC300  -  This  program  calculates  the  amount  of  thrust  and  the  accel¬ 
erations  of  the  submarine  for  the  ratio  of  commanded  speed  to  actual 
speed, 

ZC590  and  ZO69I  -  These  programs  calculate  the  error  present  when 
values  from  two  different  runs  are  supplied  according  to  the  criteria 
set  up  to  evaluate  the  effect a  of  changes. 

The  next  group  of  programs  are  used  to  determine  the  exact  natural 
frequency  and  damping  ratio  of  either  the  longitudinal  or  lateral  channel, 
in  impulse  is  used  to  excite  the  simulation  program  and  the  values  of 
pitch  angle  or  bank  angle  are  punched  on  cards.  Those  cards  and  the  so¬ 
lution  to  the  linear  characteristic  equation  are  used  to  determine  the 
numerical  values  of  the  time  response  equation. 

EClhO  and  BC200  -  These  programs  solve  for  the  roots  of  the  char¬ 
acteristic  equation  of  the  system  in  longitudinal  and  lateral  channels 
respectively, 

EClf>0  and  EC330  •  These  programs  use  the  output  of  the  above  pro¬ 
grams  and  the  time  response  cards  to  solve  for  the  coefficients  of  the 
time  response  equation. 

EC310  -  This  program  uses  the  output  of  EC150  or  EC300  as  initial 
conditions  to  solve  for  a  "least  squares  fit"  against  the  submarine 
solution  output  cards.  The  fitted  and  actual  response  is  either  listed 
or  plotted  for  comparison  purposes. 

The  above  programs  will  allow  further  research  in  the  field  of 
submarine  simulation  on  a  general  purpose  digital  computer  equipped 
with  a  FORTRAN  compiler.  The  final  program  is  complete,  flexible,  and 
can  be  used  for  both  research  and  check-out. 

One  of  the  most  important  features  is  that  a  complete  test  guide 
can  be  run  for  any  set  of  coefficients  before  the  simulator  is  built. 

Any  approximation  can  be  checked  in  a  week's  time  against  the  full  set  of 
coefficients  as  supplied  by  NSRDC.  If  plotting  subroutines  ere  available, 
graphs  can  be  made  of  any  output  parameter  using  standard  submarine  test 
inputs  as  a  reference. 
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B.  NSRDC  COEFFICIENT  RELIABILITY 

■  I—W—B  f|l  I  ■■■■  IN———— 

The  equations  of  notion  can  represent  a  particular  submarine  only 
if  the  coefficients  supplied  by  NSRDC  truly  represent  that  submarine. 
For  the  demonstration  submarine,  these  coexfieients  did  an  excellent 
job  of  simulating  the  correct  output,  in  spite  of  some  differences  be¬ 
tween  the  various  NSRDC  reports  supporting  these  coefficients. 

C.  NSRDC  EQUATION  ANALYSIS 

1.  COEFFICIENTS  -  MODIFIED  EQUATIONS 

The  standard  equations  of  motion  were  analyzed  using  a  demonstra¬ 
tion  submarine  as  an  example.  Coefficients  were  divided  into  longitu¬ 
dinal,  lateral,  and  cross-coupling  types  and  were  analyzed  as  to  impor¬ 
tance.  A  submarine  trainer  that  will  simulate  a  particular  submarine 
can  use  the  following  approximations.  Keep  the  coefficients  in  table  18. 

TABLE  18.  FINAL  NON-ZERO  COEFFICIENTS 


Axial 

Lateral 

Normal 

Roll 

Pitch 

law 

XUD 

YPD 

ZQD 

KRD 

MQD 

NRD 

XDSDS 

YRD 

ZUD 

KPD 

MUD 

NFD 

XDBDB 

TVD 

ZQ 

KVD 

M3 

NVD 

XDRDR 

IR 

zv 

KP 

m 

NR 

All 

IP 

ZSTR 

KV 

KSTR 

HP 

A12 

IV 

ZDS 

MAWQ 

NV 

A13 

YDR 

ZDB 

•  MWAW 

NVAV 

A21 

YVAV 

ZWAW 

KDS 

NPQ 

A22 

zw 

MDB 

A23 

ZRR 

MRR 

A31 

ZVR 

MTR 

A32 

MTV 

A33 

MRP 

IX,  U,  IZ,  ZB,  XG,  ZQ  NDR  is  reduced  by  20JG,  KVAV  is 

reduced  by  $0% 

The  coefficients  in  table  19  are  set  to  zero. 


TABLE  19.  FINAL  ZERO  COEFFICIENTS 


Axial 

Lateral 

Normal 

Roll 

Pitch 

law 

XQQ 

IPAP 

ZPP 

KQR 

MPP 

NQR 

XRR 

IPQ 

ZRP 

KPQ 

MQAQ 

NRAR 

XRP 

TQR 

ZAQDS 

KPAP 

MTP 

KWR 

xuu 

TIQ 

ZVAQ 

KR 

MAQDS 

NWP 

XTV 

INR 

ZAW 

KVQ 

MiV 

NVQ 

xw 

ZARDR 

ZW 

KWP 

MW 

NARDR 

XWE 

W 

ZQS 

KWR 

MQE 

NAVR 
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TABLE  19.  PINAL  ZERO  COEFFICIENTS  (cont.) 


Axil  1 

Lateral 

Normal 

Roll 

Pitch 

Yaw 

XWWE 

rSTR 

ZWE 

KSTR 

MWE 

NSTR 

XBRDRE 

YVW 

ZWAVffi 

KVW 

MWAWE 

NVW 

XDSDSE 

IRE 

ZDSE 

KDR 

MDSE 

NRE 

.m 

XVR 

IVE 

YVAVE 

YDRE 

ZVP 

KSTRE 

NVE 

NVAVE 

NDRE 

nr.  nz. 

SB. 

YB,  IG 

2.  COEFFICIENTS  -  BASIC  EQUATIONS 

The  coefficient  list  and  programing  tor  a  simplified  basic  trainer 
is  identical  to  that  above  except  for  the  following  changes. 

ZWAW,  MWAW,  MAWQ  are  set  to  zero. 

ZG  is  set  to  a  constant  value 
ZW  is  raisad  by  7i># 

MW  is  lowered  by  11# 
sinC  =  .971® 
cos©  =  1.02  -  .261  0 
sin0  =  ,9710 

cos0  =  1.02  -.281  0 

X  and  X  are  not  computed,  so  is  not  necessary. 
Therefore,  ip  is  not  integrated. 

Actually,  eliminating  ZWAW,  MWAW,  and  MAWQ  results  in  such  a  small 
reduction  in  the  total  program  complexity  that  it  would  seem  advisable 
to  keep  then  as  in  the  modified  program,  unless  such  a  reduction  would 
allow  a  smaller  core. 

3.  KINEMATIC  RELATIONS 

The  Euler  angles  should  be  computed  in  the  following  order. 


O  —  £  cos  <J>  —  V  SJN  <f> 

•  _  2  SIN  <j>-h  V  COS  <t> 

r  ~~  cos  0 

9 

<J>  =  p  +  y/  stN  & 

h.  INTEGRATION  METHODS  AND  TIMING 

The  best  integration  method  for  submarine  simulators  is  2nd  order 
Adams  which  has  the  difference  equation  given  below. 

83 


NAVTRADEVCEN  68-C-00$0-l 


* 


It  is  both  simple  and  stable  over  the  entire  operating  range  of 
computation.  The  value  of  the  integration  parameter  H  should  be  0.5 
second  for  an  accurate  modified  trainer  and  1  second  for  a  basic  maneu¬ 
vering  trainer. 

$.  NEAR-SURFACE  COEFFICIENTS 

The  use  of  near-surface  coefficients  in  a  submarine  simulator  is 
not  necessary.  For  particularly  accurate  operation  close  to  the  surface, 
MSTR  and  ZSTR  should  be  computed  as  a  function  of  depth  and  speed,  but 
keeping  the  deep  coefficients  otherwise  does  not  change  the  near-surface 
handling  qualities  to  any  great  degree. 

A  hydrodynamic  analysis  was  made  which  resulted  in  the  determin¬ 
ation  of  the  magnitude  of  the  coefficients  used  for  vave  effects  on  the 
demonstration  submarine. 

D.  MODIFIED  AND  BASIC  PROGRAMING 

One  program  vas  written  for  both  the  modified  and  basic  programs. 
This  program  was  written  in  two  parts.  The  first  part  takes  the  data 
deck  containing  all  coefficients  used  with  the  full  simulation  program, 
and  produces  another  data  deck  containing  the  reduced  coefficients. 

Thir  reduction  is  done  on  the  same  computer  as  the  one  used  for  the 
research  program. 

The  second  part  is  the  actual  simulation  program.  It  can  be  com¬ 
piled  on  any  small  general-purpose  digital  computer  that  has  8K  of  core 
and  a  basic  FORTRAN  compiler.  Usually,  the  latter  requires  the  former. 
The  input  to  this  program  is  the  CCNTR  subroutine  and  the  output  is 
twelve  variables  in  a  teletypewriter  format  for  demonstration  purposes. 

If  the  proper  I/O  subroutines  and  channels  were  available,  this  program 
could  ba  run  in  real-time. 

E.  WAVE  GENERATION 

An  investigation  was  made  into  the  simulation  of  waves  in  a  digi¬ 
tal  computer.  Waves  are  treated  on  a  statistical  basis  with  a  particular 
power  spectra  for  each  sea  state.  A  literature  search  was  made  that 
resulted  In  two  different  mathematical  models  and  a  number  of  different 
spectra  in  current  use.  A  bibliography  is  given  on  the  currant  state- 
of-the-art. 

A  particular  spectra  was  chosen  and  two  math  models  were  developed. 
Digital  programing  techniques  were  applied  to  these  two  models  and  a 
sample  sea  state  was  produced  and  checked  by  means  of  a  power  spectral 
density  run. 

F.  SHALL  COMPUTER  UTILIZATION 

A  small  jgeneral-purpose  16-hit  computer  can  be  uaed  in  real-time 
for  the  submarine  simulation  problem.  The  computer  would  have  to  have 
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8K  of  core  and  an  add  tine  of  approximately  2j#seo,  There  are  at  least 
six  lov-cost  computers  now  on  the  market  in  +M«  class* 

G,  VERIFICATION 

Finally,  a  complete  verification  program  was  performed  in  which 
the  full  range  of  submarine  maneuvering  operations  was  covered  for  the 
basic,  modified,  and  full-set  coefficients.  They  were  in  very  dose 
agreement,  witn  tho  modified  set  in  particular  having  almost  no  differ¬ 
ences  from  the  full  set. 
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SECTION  V 
DISCUSSION 


A.  HISTORY  AND  REVIEW  OF  THE  STATE-OF-THE-ART 

The  art  of  submarine  simulation  is  older  than  that  of  aircraft 
Mmulation,  This  is  a  natural  consequence  of  the  shipbuilding  industry 
which  has  used  models  in  towing  tanks  for  hundreds  of  years.  The  fact 
that  ship  performance  could  be  expressed  by  a  mathematical  model  has 
been  a  well-known  fact  to  the  Theoretical  workers  in  the  field  .ever 
since  Froude  and  Reynolds  worked  out  the  details  in  mid-nineteenth  oentuxy. 

During  VW  II,  the  science  of  feedback  control  theory  caught  up  with 
the  mathematical  structure  that  had  been  built  over  the  years.  This 
resulted  in  the  iskania  diving  trainer  of  which  four  are  still  in  use 
at  the  Submarine  School  in  New  London,  Conn.  The  computer  is  a  mechanical 
one  that  uees  compressed  air,  and  the  servoes  are  controlled  by  a  jet 
pipe  unit.  This  is  a  type  of  feedback  controller  that  consists  of  two 
jets  of  air  blowing  on  each  side  of  a  flapper  vane.  Che  is  the  controller 
and  the  other  is  the  follow-up,  Thio  idea  has  come  back  into  vogue  in 
today’s  fluidic  amplifiers. 

These  trainers  were  followed  in  the  early  193'0’s  by  the  21B20  built 
by  the  Electric  Boat  Co.  They  started  out  with  "Lamb’s  hydrodynamics” 
and  a  paper  by  the  Society  of  Naval  Architects  and  Marine  Engineers 
called  "Nomenclature  for  Treating  the  Motion  of  a  Submerged  Body  Through 
a  Fluid*,  April,  1956,  and  ended  up  with  an  analog  computer  that  had 
210  amplifiers  and  hi  multipliers.  The  trainer  simulated  four  different 
classes  of  submarines  with  two  different  motion  platforms.  However, 
only  one  platform  could  be  used  at  a  times  and  to  change  either  platforms 
or  classes  took  a  full  day  and  included  changing  cables,  panels,  and 
tuaming  junctions.  This  trainer  was  later  updated  to  the  21B20A. 

The  21B2QA  was  followed  by  the  21B56,  which  was  a  Basic  High  Speed 
Control  Trainer.  The  "Training  Effectiveness  as  a  Function  of  Trainer 
Complexity"  study  by  J.  Newton,  Electric  Boat  Division,  Groton, 
was  conducted  between  the  delivery  of  these  two  devices.  Device 
wae  the  first  Navy  trainer  to  reflect  the  findings  of  this  study.  The 
computer  size  was  reduced  by  a  factor  of  ten.  The  simulation  in  this 
trainer  was  simplified  to  the  degree  that  all  that  can  be  said  for  it 
is  that  the  indicators  move  when  the  controls  are  pushed. 

This  23B56  was  updated  to  the  21B56A  with  the  addition  of  a  ballast 
control  panel.  The  weight  and  moment:  computations  were  added  to  the 
hydrodynamic  computer  and  a  complete  simulation  was  designed.  The 
final  result  was  that  about  2$  times  as  many  amplifiers  were  added  to 
increase  the  simulation  capability. 

These  simplified  equations  of  motion  were  also  used  in  the  21A37 
Fleet  Ballistic  Missile  Trainer.  This  is  the  largest  training  device 


Conn., 
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at  New  London  and  the  calculations  for  the  position  and  velocity  of  the 
submarine  form  a  very  small  part  of  the  total  simulation. 

In  spite  of  the  fact  that  submarine  maneuvering  forms  a  small  part 
of  the  training  usefulness  of  the  21A37,  it  was  thought  to  be  necessary 
to  upgrade  the  flight  system  with  the  addition  of  casualty  flooding. 

The  upgraded  version  was  the  21A37/3.  The  main  changes  were  in  the 
simulation  for  ballast  tanks,  engine  power  for  backing  down,  and  calculation 
of  the  pitch  trigonometric  functions.  The  result  was  a  more  complex  but 
more  realistic  simulation. 

In  1967  Hydrosystems,  Inc,,  replaced  the  analog  computer  associated 
with  the  21B20A  with  a  digital  computer.  This  was  followed  by  the  21CJ> 
trainer  which  was  a  digital  trainer  from  the  start.  The  21B20A  update 
uses  a  DDP121  computer  with  a  16K  2ii-bit  core.  At  the  present  time,  it 
is  running  both  platforms  at  once  on  a  side-by-side  basis  in  the  same 
computer.  An  interesting  contrast  can  be  made  between  analog  and  digital 
equipment  with  this  trainer.  The  DDP12li  computer  sits  in  the  middle  of 
a  very  large  room  that  was  once  filled  with  analog  equipment,  Ihe  pro¬ 
graming  of  this  trainer  includes  all  of  the  NSRDC  coefficients. 

The  point  of  the  above  history  is  that  the  simulation  of  submarines 
started  out  on  a  very  complex  level.  As  a  reaction  to  this,  the  next 
generation  of  simulators  was  tco  simple  and  had  to  be  upgraded.  The 
simulation  is  beginning  to  move  back  to  the  complex  side  and  it  is  hoped 
that  this  Submarine  Equations  Study  will  strike  a  happy  medium. 

There  have  been  as  many  sats  of  equations  of  motion  as  there  have 
bee*',  submarine  trainers.  However,  a  close  examination  of  the  dynamics 
reports  seems  to  indicate  that  they  are  pretty  much  the  same  except  for 
the  degree  of  complexity  felt  necessary  by  the  system  analyst  who  actually 
programs  the  simulator. 

One  difference  that  did  come  to  light  was  two  schools  of  thought  as 
to  how  an  equation  Bhould  be  normalized.  The  coefficients  of  the  equa¬ 
tions  of  motion  are  reduced  to  non-dimensional  fora  by  dividing  through 
the  calculated  value  by  the  fundamental  units  necessary  to  remove  any 
dimension.  For  example, 

u.  -  k  sf 

where  is  forward  acceleration  and  £5  is  the  sternplanes.  This  is 
reduced  to 

-m«i=  |  SLlu-z  xSsgs  si 
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If  tho  ^fimsnsions  are  checked  on  both  3idos  of  this  new  equation 
with  in  radians,  it  will  be  seen  that  Xfi^Ss  <&»®nsionless. 

The  difficulty  is  with  £*■  term.  This  is  the  length  of  the  submarine 
squared.  Some  authorities  use  the  Characteristic  Area  "A"  which  is  the 
cross-sectional  area  of  the  submarine  and  which  has  the  same  dimensions 
as  This  procedure  comes  from  "Nomenclature  for  Treating  tho  Motion 

of  a  Submerged  Body  Through  a  Fluid",  SNAME  Technical  and  Research  Bulletin 
No,  1-5  and  is  perfectly  correct  if  the  coefficients  were  reduced  using 
this  parameter.  The  difficulty  occurs  when  the  reduction  is  made  using 
and  the  simulator  recalculates  the  original  value  using  "A"  without  ad¬ 
justing  for  the  differences  in  the  two  characteristics. 

It  Is  this  type  of  conflict  that  the  publication  of  the  Naval  Ship 
Research  and  Development  Center  publication  "Standard  Equations  of  Motion 
for  Submarine  Simulation ",  Report  2.510  will  help  eliminate.  For  the 
first  time,  a  set  of  equations  of  motion  are  available  that  can  be  con¬ 
sidered  as  a  standard,  and  this  should  ease  the  burden  on  the  system 
designer  to  a  great  degree. 

B.  DETERMINATION  OF  SIMULATION  REQUIREMENTS 

In  order  to  determine  the  simulation  requirements  for  a  submarine 
simulator,  background  material  that  would  be  useful  to  a  trainer  designer 
must  be  collected.  Investigation  includes  a  review  of  the  state-of-the- 
art,  interviews  with  submarine  personnel  and  trainer  users,  establishment 
of  data  sources,  determination  of  the  validity  of  the  data  establishing 
the  desired  ranges  and  accuracy  requirements,  and  reduction  of  this  mass 
of  material  into  performance  requirements  that  can  be  used  to  determine 
the  simplification  to  be  applied  to  the  equations  of  motion.  In  many 
cases,  the  information  gained  cannot  be  applied  directly  to  establishing 
a  performance  characteristic  but  must  be  used  when  making  a  judgment  of 
the  final  output. 

The  first  source  of  information  was  interviews  with  submarine  trainer 
users.  They  could  be  divided  into  two  classes;  namely,  those  who  were 
submariners  and  those  who  were  interested  in  the  utilisation  of  submarine 
trainers.  These  interviews  could  be  very  frustrating  because  many  men 
who  have  served  on  submarines  know  what  happens  during  a  particular  maneu¬ 
ver,  but  they  cannot  remember  how  it  happens.  For  example,  a  snap-roll 
occurs  at  the  start  of  a  tuna,  brat  it  was  difficult  to  find  out  if  this 
occurred  to  the  outside  or  the  inside  of  the  turn.  Apparently,  different 
classes  of  submarines  can  do  either.  In  general,  the  men  in  chargo  of 
training  at  the  Submarine  School  were  the  most  knowledgeable  and  could 
explain  what  should  happen. 

The  consensus  of  opinion  was  that  the  Navy  is  satisfied  with  the 
training  values  of  the  simulators  now  on  hand,  but  that  some  of  than 
are  too  simple  for  complete  training.  One  interesting  comment  was  that 
one  of  the  trainers  equipped  with  the  full  set  of  NSRDC  coefficients 
wrs  the  beet  simulator  compared  to  the  actual  submarine  that  they  had 
used.  However,  another  trainer  using  the  same  equations  of  motion  was 
not  n«j  good  and  did  not  represent  a  particular  submarine  as  well  as  the 
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th8  ether.  This  shows  that  the  NSRDC  coefficients  and  equation  cf  motion 
can  do  an  excellent  job  over  the  full  operating  range,  but  it  is  possible 
to  have  a  wrong  value  get  in  along  the  way. 

Control  in  the  longitudinal  channel  should  be  very  precise.  When 
a  300  foot  submarine  dives  at  a  35  degree  angle,  which  is  not  unusual, 
there  is  a  180  foot  difference  in  depth  between  the  bow  and  stern.  There¬ 
fore,  in  shallow  water,  extreme  care  has  to  be  taken  when  changing  depth. 
Most  maneuvering  training  involves  depth  changes  and  a  good  planesman 
can  come  onto  a  preordered  depth  with  no  overshoot  at  all.  Even  a  few 
feet  over  is  sloppy  handling.  Changes  in  course  are  not  as  much  of  a  prob¬ 
lem,  so  simulation  can  be  considerably  more  lenient  in  this  area.  There¬ 
fore,  the  philosophy  was  to  simulate  the  longitudinal  channel  as  closely 
as  possible,  but  to  allow  considerable  freedom  on  the  response  of  the 
lateral  channel  as  long  as  the  magnitudes  were  simulated  to  an  approximate 
degree. 

The  interviews  disclosed  that  there  were  a  number  of  real-world 
effects  that  did  not  transfer  into  simulator  training  through  the  equa¬ 
tions  of  motion.  This  included  such  things  as  the  tendency  on  the  part 
of  one  submarine  to  dive  when  making  turns  to  the  left  but  not  to  the 
right,  and  the  variances  between  different  submarines  of  the  same  class. 
However,  these  do  not  have  to  be  simulated  as  long  as  the  overall  response 
of  the  trainer  seems  to  be  valid  to  the  trainee. 

Until  the  present  time,  submarine  maneuvering  trainers  were  mainly 
used  for  training  at  deep  submergence,  but  near-surface  operation  is 
becoming  more  and  more  important.  Generally  the  older  trainers  had  a 
simple  sine  wave  generator  for  surface  effects,  but  this  is  not  good 
enough  any  more.  Sine  wave  generators  are  limited  because  of  their 
periodic  nature.  A  statistical  wave  generator  is  needed  so  that  an  occa¬ 
sional  wave  will  be  much  larger  than  the  others.  This  is  to  give  training 
to  prevent  the  submarine  from  broaching  when  at  periscope  depth,  which  is 
an  unpardonable  sin  to  a  submarine  captain.  The  Submarine  School  has 
only  one  simulator  with  such  a  wave  generator  in  it,  and  not  enough 
training  has  been  given  in  order  to  determine  how  effective  it  is. 

One  other  factor  was  emphasized  in  discussing  the  use  of  submarine 
simulators.  The  purpose  of  a  maneuvering  trainer  is  different  from  that 
of  an  aircraft  simulator.  In  an  aircraft,  the  pilot,  or  pilots,  is 
the  only  one  who  is  taking  action  to  control  the  vehicle.  The  purpose 
of  the  trainer  is  to  have  the  pilot  respond  instinctively  because  the 
situation  can  change  faster  than  there  is  time  to  think  about  it.  On 
the  other  hand,  a  submarine  trainer  attempts  to  train  the  officers  to 
react  in  a  team  effort.  At  least  five  to  seven  people  are  required 
to  control  the  path  of  the  submarine  and  each  one  has  to  work  as  part 
of  the  team  under  any  conditions  they  may  encounter.  As  the  amount  of 
reserve  buoyancy  on  the  newer  submarines  is  decreased,  the  actions  by 
an  individual  become  more  and  more  critical  while  the  amount  of  time 
to  respond  in  an  emergency  becomes  less.  Since  the  required  response 
to  the  same  casualty  in  two  different  classes  can  be  exactly  the  opposite, 
the  simulation  fidelity  in  an  advanced  trainer  has  to  be  valid  even  at 
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the  extreme  ends  of  the  operating  range  of  the  subaarine, 

the  second  source  for  information  to  determine  the  simulation  re¬ 
quirements  for  a  trainer  that  simulates  a  particular  submarine  is  the 
literature  available  on  that  submarine.  The  user  of  this  study  should 
be  thoroughly  familiar  vith  the  NSRDC  publications  and  the  handbooks 
that  are  available.  This  is  covered  in  part  D  of  this  section.  They 
not  only  cover  the  numerical  values  for  the  submarine,  but  they  also 
give  a  physical  description  of  any  unusual  maneuvers  that  might  be  en¬ 
countered  and  -that  would  need  consideration  for  inclusion  in  the  final 
simulation. 

A  problem  that  has  always  been  associated  with  any  simulator  is 
that  of  the  range  and  accuracy  of  the  variables  in  the  equations  of 
notion.  This  is  because  an  analog  computer  has  a  certain  definite 
resolution  and  noise  level  based  on  the  computing  voltage  used  for  full- 
scale  of  any  variable.  This  is  especially  true  under  conditions  of  level 
flight  because  many  of  the  variables,  and  thus  the  voltage,  are  close 
to  aero  and  the  computer  is  trying  to  close  the  loop  and  compute  response 
to  a  signal  that  is  equivalent  to  the  noise  present.  If  the  range  of  a 
variable  is  reduced  to  improve  the  resolution,  then  saturation  may  occur 
when  normal  maneuvers  are  undertaken. 

This  is  not  a  problem  in  a  digital  computer,  since  only  numbers  are 
being  operated  on  and  the  actions  of  the  machine  are  exactly  repeatable. 
The  accuracy  and  resolution  of  the  machine  is  limited  by  the  number  of 
bits  in  the  conputer  word.  This  is  usually  twenty-four,  even  in  a  six¬ 
teen-bit  cosputer  where  two  words  are  used,  and  so  a  change  in  the  least 
significant  bit  is  equal  to  6  x  10" a  part  of  the  total  resolution  of  each 
variable.  Therefore,  a  continuous  error  of  one  bit  at  each  integration 
cycle  et  one  second  real-time  intervals  in  computing  an  addition  would 
result  in  one-tenth  of  one  percent  error  over  fifteen  minutes  operation. 
Range  is  separate  from  resolution  in  floating  point  notation  because  an 
exponent  is  computed  along  with  the  fractional  part  of  the  word.  This 
can  be  up  to  lliO  orders  of  magnitude  even  in  a  email  computer. 

The  actual  accuracy  of  the  system  as  a  whole  is  limited  by  the  I/O 
equipment  which  was  not  covered  as  part  of  this  study.  The  output  con¬ 
version  should  be  as  accurate  as  the  trainees1  : instruments  and  as  accurate 
as  the  motion  platform  feedback  pick-ups.  The  Input  accuracy  should  be 
good  enough  so  that  an  accurate  conversion  is  made  on  the  order  of  one- 
tenth  of  one  percent  when  the  controls  are  being  moved  at  their  fastest 
rate  by  the  operator. 

The  range  of  all  variables  would  have  to  be  known  if  fixed-point 
coding  were  uaed,  Casualty  runs  would  have  to  be  made  using  the  supplied 
program  to  determine  the  maximum  ranges  under  all  conditions.  This  is 
not  necessary  if  FORTRAN  programing  is  used. 


C.  VERIFICATION  REQUIREMENTS 


There  is  a  need  for  verification  at  three  different  points  in  the 
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inres tigation  of  a  set  of  coefficients  to  be  used  in  a  submarine  simulator. 
It  should  be  determined  that  the  coefficients  match  the  submarine  being 
simulated,  that  individual  changes  made  in  the  coefficients  do  not  affect 
the  simulation  beyond  a  certain  degree,  and  that  the  complete  reduced 
simulation  does  the  same  job  as  effectively  as  the  full  set  of  coefficients. 
The  handling  qualities  of  the  simulator  should  be  the  same  as  the  submarine, 
but  the  output  to  a  particular  set  of  inputs  is  what  is  actually  measured. 
Because  the  research  program  is  limited  to  non-real  time,  measurements 
over  a  wide  spectrum  of  operation  have  to  serve  in  place  of  pilot  evaluation 
which  is  usually  used. 

Discussion  of  the  methods  of  satisfying  the  first  two  requirements 
have  been  covered  in  the  sections  of  this  study  on  impulse  testing  and 
reduction  of  coefficients.  The  procedures  of  the  third  requirement  have 
been  covered  in  the  section  on  Verification  under  Methods  and  Procedures. 

The  results  are  given  in  Advanced  Submarine  Systems  Equations,  NTDC 
Report  Number  68-C-0050-2. 

Generally,  a  concise  definition  of  how  two  sets  of  outputs  are  to 
be  compared  cannot  be  written  unt.\l  the  outputs  have  been  generated. 

However,  some  rules  were  followed  that  apply  to  most  of  the  verification 
runs  that  were  made.  These  were  a  result  of  discussions  with  Goodyear 
Aerospace's  Human  Factors  Department.  It  was  decided  not  to  apply  perfect 
phase  criteria  to  the  comparison  plots.  This  is  because  the  man  in  the 
loop  cannot  determine  the  passage  of  time  to  a  precise  degree.  Phase 
differences  show  up  either  as  a  lead  or  lag  on  the  response  plots.  What 
is  important  is  that  steady-state  values  remain  the  same  and  that  the 
magnitude  of  easily  recognizable  points  be  the  same.  For  example,  the 
snap-roll  angle  at  the  start  of  a  turn  should  be  simulated  accurately. 

In  some  cases,  the  comparison  curves  have  no  definite  points  to  compare. 

This  is  generally  true  when  the  damping  ratio  is  high,  and  in  these  cases 
the  response  was  compared  at  the  end  of  two  minutes.  Hie  difference  was 
held  to  five  percent  change  of  the  change  between  that  point  and  the  value 
at  t«0.  The  five  percent  figure  was  considered  as  below  the  threshold  of 
observation  on  the  part  of  the  planesman,  and  the  two  minutes  is  a  human 
factors  estimate  of  how  long  an  operator  in  a  closed  loop  can  remember 
the  past  actions. 

A  number  of  tests  vers  ran  as  part  of  the  major  verification  pro¬ 
cedure  incorporated  into  the  contract  for  this  study.  Volumes  III  and 
IV  contain  only  the  output  as  a  result  of  these  tests  and  the  total 
number  of  graphs  produced  was  over  550.  As  it  turned  out,  this  extensive 
testing  proved  to  be  of  benefit  to  the  study  because  it  showed  up  sevoral 
weaknesses  in  the  coefficient  reduction  that  was  used  in  the  first  mathe¬ 
matical  model.  However,  for  a  test  guide  to  be  used  with  a  simulator, 
the  output  would  not  have  to  be  nearly  as  extensive  as  was  required  in 
this  case.  The  various  submarine  tests  are  discussed  from  the  point  of 
view  of  producing  a  test  guide  for  a  training  simulator  in  the  field. 

Hie  tests  themselves  are  described  in  the  section  on  Verifications. 

Meander  -  In  this  case,  only  stemplane  angle  and  pitch  should  be 
recorded.  Depth  will  be  measured  during  the  overshoot  runs.  Instead  of 
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calculating  the  damping  ratio,  the  criterion  should  be  to  match  the  amp¬ 
litude  peaks  to  within  five  percent  of  the  amplitude  of  the  first  peak 
and  the  zero  crossings  should  bo  within  five  seconds  of  the  values  of  the 
standard.  The  five  different  speeds  are  a  suitable  number,  but  one  of 
them  should  be  at  the  point  of  maximum  damping  ratio,  which  is  approxi¬ 
mately  7  knots  in  the  demonstration  submarine. 

Vertical  Overshoot  -  These  tests  determine  the  handling  qualities  in 
the  longitudinal  channel.  They  required  forty  runs  to  cover  the  range  of 
operation.  If  the  NSRDC  coefficients  are  being  used,  a  valid  check  can  be 
mads  at  one  high  speed  and  one  low  speed.  The  output  seats  to  be  linear 
with  the  driving  plane  angle,  so  one  angle  should  be  sufficient.  Pitch 
execute  angles  should  be  five  and  twenty  degrees.  This  reduces  the  veri¬ 
fication  tests  to  four,  or  one- tenth  the  number  above.  The  only  time  that 
a  great  number  of  overshoot  tests  should  be  performed  is  when  a  guide  is 
to  be  prepared  for  training  use.  The  planesman  does  not  know  how  big  a 
change  is  required  for  any  depth  change,  particularly  on  the  basic  trainer 
which  does  not  represent  any  specific  submarine.  If  the  *  '  era  tor  is 
supplied  with  a  chart  of  "Overshoot  Pitch  Angle  versus  Approach  Speed" 
and  "Total  Depth  Change  versus  Maximum  Pitch  Angle"  for  the  basic  trainer, 
he  can  predict  what  plane  angle  changes  are  required  for  a  specified 
depth  change.  However,  for  test  guide  purposes,  the  four  curves  should 
be  enough  if  they  are  matched  on  the  same  basis  as  the  meander  curves. 

Steady  Submerged  Turns  -  These  verification  tests  required  that 
eleven  parameters  be  recorded  at  three  rudder  angles  and  five  different 
speeds  for  a  total  of  165  charts.  Much  of  this  information  is  not  used 
by  the  system  programsr  to  check  the  operation  of  the  simulator  and  the 
following  is  suggested.  The  rudder  angles  should  be  15°  and  35°  and  the 
speeds  should  be  $  and  25  knots.  The  parameters  recorded  as  a  function  of 
time  are*  rudder  angle,  stemplane  angle,  speed,  heading,  roll,  and  depth. 
Another  requirement  is  that  these  turns  should  be  made  under  autopilot 
control,  using  the  stem  planes  to  close  the  loop  in  keeping  a  constant 
depth.  The  point  checks  are  snap-roll  angle,  steady-state  roll  angLe, 
change  in  speed  at  steady-state  conditions,  and  heading  change  at  the  end 
of  a  specified  period  of  time. 

The  reason  for  using  autopilot  oontrol  and  recording  stemplane' 
angle  and  depth  is  that  although  these  do  not  change  very  much  during  a 
steady-state  turn,  their  shape  over  time  is  very  important  for  checking 
ail  six  axes  in  the  simulation.  Open  loop  control  will  not  serve  because 
the  longitudinal  channel  is  frozen  and  the  effects  of  lateral  motions  an 
elevation  ijs  not  apparent.  These  effects  can  be  quite  large  with  a  small 
change  in  some  of  the  coefficients. 

Acceleration/Deceleration  -  Two  tests  are  given  in  this  study, 
although  only  one  was  required.  The  test  given  in  the  specifications 
requires  several  minutes  of  computer  time  to  run  and  it  is  felt  that 
the  maximal  acceleration/deceleration  test  will  verify  the  same  operation 
in  a  shorter  time.  The  submarine  is  at  zero  speed  and  a  command  speed  of 
maximum  is  called  for.  Wien  this  speed  is  reached,  a  command  speed  of 
aero  is  called  and  the  program  runs  for  the  same  length  of  time  again. 
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The  curve  of  speed  versus  time  should  match  to  a  close  degree  to  insure 
good  thrust  simulation. 

Longitudinal  Trim  -  This  tests  for  conditions  that  should  be  verified 
under  meander  and  overshoot. 

In  this  discussion,  it  is  assumed  that  the  NSRDC  coefficients  are 
correct,  and  the  approved  criteria  is  determined  by  the  output  of  the 
programs  produced  by  this  study.  In  fact,  such  runs  on  a  general-purpooe 
digital  computer  should  be  made  part  of  the  specifications  for  any  new 
submarine  trainer.  Not  only  will  a  test  guide  be  produced,  but  also  a 
means  of  producing  data  will  be  provided  to  assist  in  maintenance  and 
trouble-shooting  during  the  life  of  the  simulator, 

D.  NSRDC  DATA  SOURCES 


The  Naval  Ship  Research  and  Development  Canter  (NSRDC)  has  a  complete 
data  package  on  about  25  different  existing  submarine  designs.  Most  of 
this  data  carries  the  classification  of  “CONFIDENTIAL*  so  that  before  it 
can  be  furnished  to  •'.n  outside  contractor  the  proper  clearances  and  need- 
to-know  has  to  be  established  through  normal  channels.  Cnee  this  is  done, 
a  complete  set  of  data  can  be  requested  from  NSRDC,  The  programs  devel¬ 
oped  in  this  study  use  the  NSRDC  mathematical  model  as  a  base,  and  there¬ 
fore  such  a  data  package  is  essential  before  simulating  any  particular 
submarine.  This  section  describes  what  is  available  in  the  set  using  the 
demonstration  submarine  as  an  example.  The  name  and  number  has  been  re¬ 
placed  with  NAME  and  #  in  the  titles  of  reports.  The  titles  are  the  same 
except  for  this  so  any  report  can  be  requested  by  inserting  the  proper 
name  and  number  in  the  request  to  NSRDC  or  the  Naval  Training  Device 
Center  Technical  Library  at  Orlando,  Florida,  They  can  also  be  received 
through  the  Defense  Documentation  Center  in  Virginia  with  the  necessary 
clearances. 

The  essential  report  for  the  simulation  of  the  equations  of  motion 
is  the  Naval  Ship  Research  and  Development  Center  Report  2610;  “Standard 
Equations  of  Motion  for  Subsarine  Simulation*.  This  unclassified  report, 
which  was  described  above,  contains  the  equations  of  motion  directly 
applicable  to  submarine  simulation  and  soma  description  of  their  back¬ 
ground,  how  tc  use  them,  and  how  data  is  generated  in  order  to  determine 
the  coefficients  for  a  particular  submarine. 

The  equations  of  motion  listed  in  this  report  are  applicable  to 
the  rigid  body  motions  of  any  submarine  or  submerged  vehicle.  They  are 
writ-can  to  use  non-dimensional  coefficients,  and  it  is  those  coefficients 
that  are  needed  to  simulate  realistically  a  speoific  submarine.  A  complete 
set  of  input  data  consists  of  the  dimensions  of  the  submarine  being  simu¬ 
lated,  the  values  of  the  inertia  terms,  and  values  for  the  hydrodynamic 
coefficients. 

The  numerical  values  can  be  requested  from  NSRDC  under  the  title 
"Hydrodynamic  Coefficients  and  Other  Terms  for  Equations  of  Motion  for 
#(U)“,  They  consist  of  three  pages  of  preprinted  forms  containing  spaces 
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for  all  the  ooeffioients  listed  in  NSRDC  Report  2 $10  along  with  the  dimen¬ 
sions,  buoyancy  and  inertial  terras.  They  do  not  contain  values  for  what 
is  called  systems  simulation  in  a  complete  submarine  simulator.  The 
addition  terms  nesded  if  a  real-time  closed  loop  simulation  is  to  be 
attested  are  such  values  as  thrust  time  delays,  oontrcl  surface  rstes, 
and  time  histories  of  the  ballast  tanks  during  flowing  and  venting,  With 
this  list  of  numerical  values  and  the  equations  in  Report  2510,  it  is 
possible  to  simulate  any  time  outputs,  as  a  result  of  any  specified  input 
conditions,  on  a  digital  computer. 

The  following  reports  are  not  needed  if  the  equations  were  programed 
without  error  and  the  coefficients  are  correct  and  accurately  represent 
the  submarine  being  simulated.  Since  this  is  usually  not  the  case,  it  is 
necessary  to  hare  information  against  a  sot  of  input  conditions.  Hydro¬ 
dynamic  coefficients  can  be  calculated  generally  from  two  different  sources: 
towing  tank  tests  and  full  scale  sea  trials,  A  report  generally  covers 
each  of  these  areas, 

"A  Full-Scale  Evaluation  of  the  Handling  Qualities  of  the  NAME  Class 
Submarine  (#)  (U)«  gives  the  results  of  the  full  scale  sea  trails,  A 
number  of  standard  tests  have  been  developed  for  testing  the  handling 
qualities  of  a  submarine.  They  cover  submerged  performance  in  both  the 
vertical  plane  and  the  horizontal  plane.  This  report  givos  the  specifica¬ 
tions  of  the  submarine,  describes  the  procedures  used  to  conduct  the  tests 
in  each  plane,  gives  samples  of  outputs  to  typical  test  runs,  and  summarizes 
the  results  over  the  range  of  operation  of  the  submarine.  Since  the  input 
conditions  are  specified,  the  full  scale  tests  can  be  repeated  with  the 
computer  program  and.  compared  with  tabulated  results  to  see  if  the  outputs 
of  the  simulation  are  accurate. 

"Model  Investigation  of  the  Stability  and  Control  Characteristics  of 
HAMS  (#)  at  Deep  Submergence  (U)"  contains  the  results  of  model  testing 
and  a  list  of  the  coefficients  developed  from  these  tests.  They  should  be 
the  same  as  the  list  obtained  from  NSRDC  except  for  the  fact  that  work  is 
always  being  performed  to  refine  the  numerical  values.  This  report 
approaches  the  subject  of  response  on  a  frequency  basis  rather  than  on  a 
time  basis  as  in  the  full  soale  trails.  The  outputs  of  the  model  tests 
are  used  as  solutions  to  a  set  of  differential  equations  and  the  various 
roots  of  the  characteristic  equation  are  used  to  determine  the  coefficients. 
These  roots  can  be  used  to  determine  the  characteristic  frequency  and 
damping  ratios  as  a  function  of  speed.  The  simulation  program  produces 
the  same  information  in  response  to  impulse  testing  and  thus  a  comparison 
can  be  made.  This  report' also  discusses  control  effectiveness,  neutral 
angles  and  critical  speed.  It  is  essential  that  these  be  correct  for  the 
simulation  because  they  determine  steady-state  operation  and  the  trainee 
will  notice  errors  in  this  area  before  anything  else. 

The  numerical  results  of  the  coefficients  given  by  the  Wo  above 
can  be  compared  in  "Full  Scale  Evaluation  of  Static  Stability  and  Control 
Derivatives  of  NAME  (#)  (U)«.  This  report  gives  plots  of  the  model  curve 
against  the  full  scale  data  points. 
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If  any  work  is  to  be  performed  in  the  area  of  near-surface  operation, 
■toe  report  "Simulator  Studies  of  Depthkeeping  Ability  of  the  NAME  (#) 

Under  a  Seaway  (U)"  is  helpful.  This  report  contains  the  changes  in  the 
coefficients  nocessary  to  actually  simulate  operation  undor  a  free  surface, 
and  submarine  responses  to  waves  while  trying  to  hold  an  ordered  depth. 

This  publication  also  contains  the  constants  for  a  NSRDC  autopilot  that 
proved  helpful  to  the  present  study, 

"Submerged  Turning  and  Maneuvering  Characteristics  of  the  #  Suomarine 
from  Free  Running  Model  Tests"  supports  the  figures  given  in  the  full 
scale  reports  in  the  horizontal  plane.  Heading  changes,  loss  of  speed, 
and  snap  roll  angles  are  given  for  the  model, 

"Resistance  and  Propulsion  Characteristics  of  the  (#)  Submarine  as 
Predicted  from  Tests  of  Model  and  NAME  (#)"  and  "Machinery  Performance 
Trials  (U)"  are  the  model  testing  and  full  scale  reports  respectively 
that  relate  the  command  speed  used  by  the  NSRDC  equations  of  motion  with 
the  revolution  per  minute  (RPM’s)  seen  by  the  engines.  They  can  also  be 
used  to  determine  control  surface  deflection  rates. 

Finally,  there  are  a  number  of  special  reports  covering  recovery 
procedures  from  casualty  conditions  which  vary  from  submarine  to  submarine 
depending  on  the  amount  of  work  that  has  been  done  in  this  area. 

The  one  area  that  has  not  been  covered  in  these  reports  is  that  of 
tank  simulation.  In  order  to  know  the  shift  in  center  of  gravity  and 
weight  change  due  to  the  ballast  tanks,  it  is  necessary  to  know  their 
location  and  capacity.  This  information  can  be  found  in  the  NAVSHIPS 
■Ship  Information  Books"  and  "Training  Aid  Booklets".  They  cover  hull 
arrangements,  controls  systems,  piping,  electrical  systems,  fire  control 
and  the  other  systems  on  the  actual  submarine. 

E.  INCORPORATION  OF  CASUALTIES 


Training  for  casualties  is  an  important  part  of  a  maneuvering  trainer. 
This  study  did  not  examine  response  to  casualties  as  such,  but  provisions 
were  made  to  incorporate  casualties  into  the  computer  programing.  The 
information  on  typical  casualties  comes  from  a  study  performed  by  Goodyear 
Aerospace  for  NTDC  titled  "Submarine  Casualty  Control  Training". 

Typical  casualties  are  ship  command  and  control  failures,  flooding, 
and  air  bank  failures. 

Ship  command  and  control  casualties  are  those  concerned  with  the 
diving  planes  and  rudder.  Such  failures  include  failure  of  the  planes  or 
rudder  either  at  a  fixed  position  or  after  moving  to  the  end  of  their 
travel.  Sternplans  angle  DS,  sailplane  angle  DB,  and  rudder  angle  DR, 
are  located  in  COMMON  of  the  FORTRAN  program.  This  means  that  these  var¬ 
iables  can  be  controlled  from  any  subroutine  that  has  an  identical  common 
statement.  At  the  present  time,  these  variables  are  scheduled  from  the 
CONTR  subroutine;  but  in  a  submarine  trainer,  they  would  be  calculated 
from  the  Input/Output  subroutine.  It  would  be  a  simple  matter  to  have  an 
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additional  input  to  the  I/O  subroutine,  representing  an  instructor’s 
failure,  to  interrupt  the  normal  control  channel. 

Flooding  can  be  incorporated  into  the  tank  simulation  system  as  an 
additional  set  of  tanks  whose  fill  starting  time  is  controlable  by  the 
instructor.  The  program  can  be  as  simple  or  as  complex  as  desired,  but 
the  final  results  affect  the  flight  section  of  the  programing  through  the 
terms  shifting  the  center  of  gravity  and  through  the  mass  and  weight 
computation.  Again,  these  variablos  are  located  in  the  COMMON  area  so 
that  reprograming  of  the  present  subroutine  is  not  required. 

Air  bank  failures  are  incorporated  into  the  normal  and  emergency 
main  ballast  tank  blow  systems,  variable  ballast  tank  systems,  and  trim 
and  drain  systems.  Since  simulation  is  not  provided  for  these  systems 
as  part  of  this  study,  these  casualties  would  have  to  be  programed  when 
the  system  simulation  was  done. 

F.  USES  AND  LIMITATIONS  OF  THE  SUBMARINE  EQUATION  STUDY 

1,  USEFULNESS 

The  usefulness  of  this  study  has  bean  indicated  by  the  topics  covered 
up  to  this  point,  but  it  is  worthwhile  to  restate  them  here. 

a.  Coefficient  Cheoks  -  A  set  of  coefficients  can  be  checked 
against  any  available  data  on  the  real  submarine.  If  new  data  is  received, 
further  cheoks  can  be  made  easily.  No  simulation  is  required. 

b.  Research'  in  Reduction  and  Casualties  -  The  programs  can  be 
used  to  investigate  the  effect  of  further  reductions  and  casualties. 

This  work  could  be  extended  with  the  aid  of  a  real-time  computing  system 
to  put  the  man  into  the  loop  in  real-time. 

c.  Test  Guides  -  A  complete  set  of  data  points  can  be  run  off 
for  any  simulator  for  which  the  coefficients  are  known.  Any  range  of 
operation  can  be  reproduced  in  a  very  short  period  of  time  once  the  coef¬ 
ficients  have  been  justified.  Such  data  outputs  are  useful  in  initial 
checking  of  the  programing  for  the  simulator  and  in  maintenance  and  trouble¬ 
shooting  in  the  field.  Both  static  and  dynamic  checks  can  be  made. 

d.  Compact  Programing  -  When  programing  is  done  for  a  submarine 
simulator,  as  many  constants  as  possible  are  combined  into  a  single  term 
to  save  on  the  amount  of  computation  required  in  the  compact  mathematical 
model.  These  programs  will  assist  in  this  task  by  calculating  the  new 
coefficients  required. 

2,  LIMITATIONS 

The  biggest  limitation  in  this  study  is  that  operation  at  the  extreme 
ends  of  the  operating  range  was  not  covered.  This  is  usually  due  to 
casualties,  such  as  flooding,  when  attitudes  are  encountered  that  would 
never  be  entered  under  normal  operation.  However,  NSRDC  states  that  the 
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equations  of  motion  include  effects  due  to  “backing"  on  propellers  and 
emergency  blowing  of  ballast  tanks  even  to  such  extremes  as  an  emergency 
recovery  from  a  hard-climb  st emplane  jam,  including  the  buoyant  ascent 
after  blowing  of  bsllast.  It  is  felt  that  the  modified  equations  of 
motion  will  respond  effectively  to  these  conditions,  even  though  they 
were  not  tested  in  this  range,  but  that  the  basic  set  probably  would  not. 
There  are  also  a  number  of  subroutines  that  would  have  to  be  written  in 
order  to  make  tests  such  as  emergency  recoveries.  They  are  listed  below. 

a.  Tanks  -  The  tank  blow  system  is  not  simulated.  However,  the 
effect  of  any  tank  arrangement,  no-c  including  partially  filled  tank  slosh, 
can  be  computed  through  Ye  ,  and  Z$» 

b.  Controls  -  Hie  time  history  for  rudder  and  plane  changes  are 
not  included  except  that  schedules  for  all  of  the  standard  submarine  test 
maneuvers  are  provided. 

c.  Failures  -  Failures,  such  as  plane  jams,  are  not  included. 
However,  many  of  the  hardeners  and  jams  of  the  planes  can  be  simulated  by 
uce  ctf  the  proper  constants  in  the  CONTR  subroutine  of  the  main  simulation 
program.  Air  Bank  end  flooding  failures  would  also  have  to  be  added. 
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SECTION  VI 


CONCLUSIONS 


The  "Study,  Advanced  Submarine  Systems  Equations"  was  a  success. 
The  "Standard  Equations  of  Motion  for  Submarine  Simulation"  NS ROC 
Report  2510  was  used  as  a  basis  for  producing  a  submarine  simulation 
program  that  could  be  utilized  in  a  small  general-purpose  digital 
computer. 
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NOTATION 


Symbol 

Dimensionless  Form 

Definition 

1  &i» 

Seta  of  propallor  thrust  constants 

B 

B1  =  ^ 

*P^U> 

Buoyancy  force,  positive  upward 

CB 

Center  of  buoyancy  of  submarine 

CG 

Ix 

I  '  = 
x  *P*6 

Center  of  mass  of  submarine 

I 

X 

Moment  of  inertia  of  submarine  about  x  axis 

I 

y 

ly 

ly  ~  ipiB 

Moment  of  inertia  of  submarine  aboui  y  axis 

V 

l 

*  Z 

i  •=  lz 
*  if*6 

Moment  of  inertia  of  submarine  about  z  axis 

T 

*xy 

I  '  r  Ixy 

xy  ipt8 

Product  of  inertia  about  xy  axis 

I 

y* 

I  «  =  Jx2_ 
y*  ipi6 

Product  of  inertia  about  yz  axes 

j  ,  _ 

zx  ip*6 

Product  of  inertia  about  zx  axes 

K 

K'  -  K 

K  'SPu1 

Hydrodynamic  moment  component  about  x 
axis  (rolling  moment) 

K* 

K'-  K* 

*  ip*3uz 

Rolling  moment  when  body  angle  (ry,  jS)  and 
control  surface  angles  are  zero 

K, 

K*  '  =  — *22.. 

.  *n  |p*3uz 

Coefficient  used  in  representing  K*  as  a 
function  of  (rj-1) 

K 

P 

K  '  =  — -^P— 

p  ipt4u 

First  order  coefficient  used  in  representing 

K  as  a  function  of  p 

K. 

.  P 

y  «  - 
»  '  »*• 

Coefficient  used  in  representing  K  as  a  function 
of  -  p 

KpiPt 

.  Kptp| 

4 

Second  order  coefficient  used  in  representing 

K  as  a  function  of  p 

Kpq 

K  '  r=  --PS 

pq  ip*B 

Coefficient  used  in  representing  K  as  a  function 
of  the  product  pq 
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Kqr 

K  ■=-?%. 
q*  fcpi6 

Coefficient  used  in  representing  K  as  a 
function  of  the  product  qr 

Kr 

K  <  - 

r  “ipfU 

First  order  coefficient  used  in  representing 

K  as  a  function  of  r 

Kf 

K  ,  -  -iSi 

Kf  “St5- 

Coefficient  used  in  representing  K  as  a  . 
function  of  f 

Kv 

K  •- 

v  "  ipPTT 

First  order  coefficient  used  in  representing 

K  as  a  function  of  v 

Kv 

»  K.'  s  ,  v 

v  hi* 

Coefficient  used  in  representing  K  as  a 
function  of  v 

Kv|v| 

y  ,  _  KvM 

Kv|vi  -f&r 

Second  order  coefficient  used  in  representing 

K  as  a  function  of  v 

Kvq 

K  ..-*SL 
vq 

Coefficient  used  in  representing  K  as  a  function 
of  the  product  vq 

K 

VW 

K  •  -  Kvw 

VW 

Coefficient  used  in  representing  K  as  a  function 
of  the  product  vw 

Kwp 

Kwd 

K  '  =  -- P 

wp  ipt* 

Coefficient  used  in  representing  K  as  a  function 
of  the  product  wp 

Kwr 

vr  .  -  ^wr 
W*.'*P*< 

Coefficient  used  in  representing  K  as  a  function 
of  the  product  wr 

K. 

6r 

„  . _  K6r 

fcpi3U2 

First  order  coefficient  used  in  representing 

K  as  a  function  of  6r 

l 

*'  =  1 

Overall  length  of  submarine 

m 

_.i  _  «  m 
m  '  .3 

tpt3 

Mass  of  submarine,  including  water  in  free- 
flooding  spaces 

M 

M 

"*p43U2 

Hydrodynamic  moment  component  about  y  axis 
(pitching  moment) 

M* 

M* 

M*  "tp4SU* 

Pitching  moment  when  body  angles  (or,  S )  and 
control  surface  angles  are  zero 

M 

PP 

M  *=^$- 
pp  ipir 

Second  order  coefficient  used  in  representing 

M  as  a  function  of  n.  First  ordeT  coefficient  is 

u 

zero. 

M 

q 

First  order  coefficient  used  in  representing 

M  as  a  function  of  q 

Mq>) 

m  •  =  J^aiL. 
q v  *pl*u 

First  order  coefficient  used  in  representing 
M^Cs  a  function  of  (fj-1) 

M. 

q 

Mi 

M.’  =  r-rf- 

q  5P<-6 

Coefficient  used  in  representing  M  as  a 
function  of  q 
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* 


M=|q| 

MqU! 

-  Ms! 

Second  order  coefficient  used  in  representing 

M  as  a  function  of  q 

“W»« 

M.  .  ;.3sl5* 
Ul68 

Coefficient  used  in  representing  M,  as  a 
function  q  * 

M 

*P 

M  '  = 
rp 

wrp 

ipt6 

Coefficient  used  in  representing  M  as  a 
function  of  the  product  rp 

M 

rr 

M 

vp 

M  '  = 
rr 

%'  = 

Mrr 

Second  order  coefficient  used  in  repress  .  g 

M  as  a  function  of  r.  First  order  coefficient 
is"  zero 

Coefficient  used  in  representing  M  as  a 
function  of  the  product  vp 

•» 

Myp 

ipt4 

Mvr 

Mvr'  = 

Myr 

ipt4 

Coefficient  used  in  representing  M  as  a 
function  of  the  product  vr 

Mvv 

M  '  = 
vv 

MVV 

Second  order  coefficient  used  in  representing 

M  as  a  function  of  v 

M 

w 

M  1  = 

w 

Mw 

iptAu 

First  order  coefficient  used  in  representing 

M  as  a  function  of  w 

M 

.  WT) 

MW: 

=  —  xn.. 

ipru 

First  order  coefficient  used  xn  representing 

Mw  as  a  function  of  H-l) 

M. 

w 

Mw’  = 

Mw 

Coefficient  used  in  representing  M  as  a  function 
_  of  w 

M  ■  | 

M 

M.  , 

jwjq 

M  •-  M!w| 

mM  iptau 

w  ,  _  M|w  U 

mM<i  £p*4 

First  order  coefficient  used  in  representing  M 
as  a  function  of  w;  equal  to  zero  for  symmetrical 
function 

Coefficient  used  in  representing  Mq  as  a  function 
of  w 

Mw|w| 

0 

U  I 
w  1  w 

,  |w  1 

!  -  *pt*  ■ 

Second  order  coefficient  used  in  representing 

M  as  a  function  of  w 

M  .  , 
w|w|TJ 

M  . 
w|  w 

First  order  coefficient  used  in  representing 

M  i  ,  as  a  function  of  (77- 1 ) 

Mww 

M  • 
ww 

Mww 

Second  order  coefficient  used  in  representing 

M  as  a  function  of  w;  equal  to  zero  for  sym« 
metrical  function 

M6b 

•  M6b*  = 

Mjb 

ipFu2 

First  order  coefficient  used  in  representing 

M  as  a  function  of  6^ 

M6. 

M6s 

First  order  coefficient  used  in  representing 

M  as  a  function  of  6, 

\ 

Jpi3U2 

m4.„ 

M63T}' 

Mftsn 

First  order  coefficient  used  in  representing 

M$5  as  a  function  of  (77- 2) 
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N 

x„  N 

N  -^Fu2 

Hydrodynamic  moment  component  about  z 
axis  (yawing  moment) 

N* 

N* 

N*'=-££FiP 

Yawing  moment  when  body  angles  (ot,  /?)  and 
control  surface  angles  are  zero 

N 

P 

Np 

p  JpFiT 

First  order  coefficient  used  in  representing  N 
as  a  function  of  p 

N. 

P 

N- 

N.'  =  T-fi— 

P  ipl° 

Coefficient  used  in  representing  N  as  a  function 
of  p 

N 

pq 

N  '  -  -^E£L 

pq  ipie 

Coefficient  used  in  representing  N  as  a  function 
of  the  product  pq 

N 

V 

*  Nqr 

V  tP-t6 

Coefficient  used  in  representing  N  as  a  function 
of  the  product  qr 

N 

r 

N 

N  '= 

r  ip44u 

First  order  coefficient  used  in  representing  N 
as  a  function  of  r 

N 

TV 

n  1  =  3a_ 

rt}  4pt,*U 

First  order  coefficient  used  in  representing 

Nr  as  a  function  of  (rj- 1) 

N. 

r 

N* 

N. 1  =  -—V- 
r  ipt6 

Coefficient  used  in  representing  N  as  a  function 
of  r 

N  .  . 

N  Nr,r| 

Nrlr|  “ 

Second  order  coefficient  used  in  representing 

N  as  a  function  of  r 

^|r| fir 

N  '  ,  _  N|rl6r 
|r|6r  Jpt«U 

Coefficient  used  in  representing  N  as  a 
function  of  r  r 

N 

V 

Nv 

N  '  =  v 

v  $pi3U 

First  order  coefficient  used  in  representing  N 
as  a  function  of  v 

N 

VTJ 

Nvn 

N  1  s  — 12L 

V7J  7£p43U 

First  order  coefficient  used  in  representing  Nv 
as  a  function  of  (fj-1) 

N. 

V 

M* 

N. '  =  * 

v  ip*4 

Cocfficv-  used  in  representing  N  as  a 
functic  .  v 

N 

vq 

n  -  =il^a_ 
yq  ipt* 

Coefficient  used  in  representing  N  as  a  function 
of  the  product  vq 

N,  , 

1  vl  r 

*  ,_Nlv|r 

Vlr  - 

Coefficient  used  in  representing  Nr  as  a 
function  of  v 

N  .  . 
v|v| 

XT  ,  _  Nv|vf 

NyM  ^pt3 

Second  order  coefficient  used  in  representing 

N  as  a  function  of  v 

Nv|v|»j 

XT  ,  nv|v!tj 

vlvh  ipt3 

Fitjst  order  coefficient  used  in  representing 
Nvjv,  as  a  function  of  (tj-1) 

103 


NAVTRADEVGEN  68-C-0Q50-1 


,  _  ^vw 


Coefficient  used  in  representing  N  as  a  function 
of  the  product  vw 


N  '  =  t-t^- 

wp 


Coefficient  used  in  representing  N  ae  a  function 
of  the  product  wp 


N6r'  =  U* 

*  N  ■  =  3j3L. 

tot) 


Coefficient  used  in  representing  N  as  a  function 
of  the  product  wr 

First  order  coefficient  used  in  representing  N 
as  a -function  of  6r 

First  order  coefficient  used  in  representing 
N^r  as  a  function  of  (tj-1) 

Angular  velocity  component  about  y  axis 
relative  to  fluid  (roll) 


p  ’I F 


Angular  acceleration  component  about 
relative  to  fluid 


Angular  velocity  component  about  y  axis  relative 
to  fluid  (pitch) 


-I  _ 

r  --jj- 


l.  .  H* 


Angular  acceleration  component  about  y  axis 
relative  to  fluid 

Angular  velocity  component  about  z  axis 
relative  to  fluid  (yaw) 


Angular  acceleration  component  about  z  axis 
relative  to  fluid 


U'-Hi 


Linear  velocity  of  origin  of  body  axes  relative 
to  fluid 


-  u 
U 


Component  of  U  in  direction  of  the  x  axis 


ft,  _  ut 


Time  rate  of  change  of  u  in  direction  of  the 
x  axis 


u  '  =  — 
c  U 


Command  speed:  steady  value  of  ahead  speed 
component  u  for  a  given  propeller  rpm  uhen 
body  angles  (a,  8)  and  control  surface  angles 
are  zero.  Sign  changes  with  propeller  reversal 


Component  of  U  in  direction  of  the  y  axis 


=  iL 

u* 


Time  rate  of  change  of  v  in  direction  of  the 
y  axis 


M ATf'ro  a nrirr t»t  £q  n  r\^rv\  i 

4WM  iJLUUMfl  UCM  UU**U*OVpU*i 


w 

,  w 

w*  = - 

U 

Component  of  U  in  direction  of  the  z  axis 

w 

.  ,  wt 

U1 

Time  rate  of  change  of  w  in  direction  of  the 
z  axis 

W 

W-  -  W 

W  -iw’iF 

Weight,  including  water  in  free  flooding  spaces 

X 

4 

Longitudinal  body  axis;  also  the  coordinate  of  a 
point  relative  to  the  origin  of  body  axes 

XB 

A* 

ii 

> 

» 

The  x  coordinate  of  CB 

XG 

The  x  coordinate  of  CG 

*o 

oU 

X  1 

II 

'o 

X 

X  coordinate  of  the  displacement  of  CG  relative 
to  the  origin  of  a  set  of  fixed  axes 

X 

X-  -  x 
"  *p42U2 

Hydrodynamic  force  component  along  x  axis 
(longitudinal,  or  axial,  force) 

x  <  = 
qq  £p44 

Second  order  coefficient  used  in  representing 

X  as  a  function  of  q.  First  order  coefficient 

is  zero 

X 

rp 

V  .  _  ■  xrp 
rP  *p44 

Coefficient  used  in  representing  X  as  a  function 
of  the  product  rp 

X 

rr 

X.r 

X  •  =  -j-Ii- 
rr  $p44 

Second  order  coefficient  used  in  representing 

X  as  a  function  of  r.  First  order  coefficient  is 

zero 

X. 

u 

X£ 

Xu*=i5F- 

Coefficient  used  in  representing  X  as  a  function 
of  u 

Xuu 

x  '  =  -£uu. 
uu  *p42 

Second  order  coefficient  used  in  representing 

X  as  a  function  of  u  in  the  non-propelled  case. 
First  order  coefficient  is  zero 

Xvr 

Xvr 

X  •  —  v* 
vr  -£p*3 

Coefficient  used  in  representing  X  as  a  function 
of  the  product  vr 

Xw 

v  t  _  Xyv 
w 

Second  order  coeffient  used  in  representing  X 
as  a  function  of  v,  First  order  coefficient  is  zero 

% 

XVvr> 

JW = 

First  order  coefficient  used  in  representing  Xyv 
as  a  function  of  (tj-I) 

Xwq 

r  *  _  xwq 

wq  “  i p43 

Coefficient  used  in  representing  X  as  a  function 
of  the  product  wq 
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Xww 

v  ,  _  xww 

AWW  ^pt2 

Second  order  coefficient  used  in  representing 

X  as  a  function  of  w.  First  order  coefficient  is 

zero 

x 

UWtJ 

x  <  -  Xwwfl 

wwt)  ^p<,i 

First  order  coefficient  used  in  representing  Xww 
as  a  function  of  (77- 1) 

X6b6b 

v  ,  X6b6b 

6b6b  "jpt2U2 

Second  order  coefficient  used  in  representing  X 
as  a  function  of  6^,  First  order  coefficient 
is  zero 

X6r6r 

v  i  _  X6r6r 

«A6r6r  '  £pt2U2 

X  '  =  -Jeter] 

or6r?j  Ifpt^U1 

Second  order  coefficient  used  in  representing 

X  as  a  function  of  6r.  First  order  coefficient  is 
zero 

1 

X6r6rj) 

First  order  coefficient  used  in  representing 
X6r6r  as  a  ^unct*on  °*  (*7~1) 

X8s6s 

v  .  X3s6s 

6*6b  _  Ipi2 Ul 

Second  order  coefficient  used  in  representing  X 
as  a  function  of  6S.  First  order  coefficient  is 

zero 

X6s6»n 

v  .  _  x5r68T? 

6s  6*77  Jp^2U2 

First  order  coefficient  used  in  representing 
X6s6s  aS  a  *unct*on  °*  (tj-  1) 

y 

y. 

Lateral  body  axis;  also  the  coordinate  of  a 
point  relative  to  the  origin  of  body  axes 

yB 

V  ’-VB 

yB  -1“ 

The  y  coordinate  of  CB 

yG 

,  yC 
yG  =  T~ 

The  y  coordinate  of  CG 

yo 

A- 

II 

> 

A  coordinate  of  the  displacement  of  CG  relative 
to  the  origin  of  a  set  of  fixed  axes 

Y 

Y'  -  Y 

Hydrodynamic  force  component  along  y  axis 
(lateral  force} 

Y, 

Y*'  =  Y  - 
*pi2  U2 

Lateral  force  when  body  angles  (a,  /3)  and  control 
surface  angles  are  zero 

Y 

P 

YP 

Yp  ~*p43U 

First  order  coefficient  used  in  representing 

Y  as  a  function  of  p 

Y. 

P 

Coefficient  used  in  representing  Y  as  a  function 
of  p 

Yplpf 

Y  ,  , '  =  Yp  M 

p!p> 

Second  order  coefficient  used  in  representing 

Y  as  a  function  of  p 
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Y 

pq 

v  •  =  -iP.S- 

pq 

Coefficient  used  in  representing  Y  as  a  function 
of  the  product  pq 

Yqr 

Ynr 

V  * 

Coefficient  used  in  representing  Y  as  a  function 
of  the  product  qr 

Yr 

Yr  "ipt3  U 

First  order  coefficient  used  in  representing  Y 
asa  function  of  r 

Yrrj 

Y  <  -  YXT? 

yxtj 

First  order  coefficient  used  in  representing 

Yr  as  a  function  of  (?J-1) 

Y* 

r 

Yj 

Yr  '  4p^ 

Coefficient  used  in  representing  Y  as  a  function 
of  r 

Y|rUr 

Coefficient  used  in  representing  Y.  as  a 
function  of  r 

Y 

V 

y..  yv 
v  $pt*u 

First  order  coefficient  used  in  representing 

Y  as  a  function  of  v 

Yv»7 

y  •  =  -2yn. 
vv  ip^2u 

First  order  coefficient  used  in  representing 

Yv  as  a  function  of  (rj-1) 

Y. 

V 

y  •  -  Y* 

v  *pA3 

Coefficient  used  in  representing  Y  as  a 
function  of  v 

Yvo 

V s  # 

Coefficient  used  in  representing  Y  as  a  function 
of  the  product  vq 

Yvl  r  | 

V  .  -  T*M 

VM  ipt3 

Coefficient  used  in  representing  Yv  as  a  function 
of  r 

Yv|  vl 

Yy|vl#  ip? 

Second  order  coefficient  used  in  representing 

Y  as  a  function  of  v 

Yvlv|jji 

„  ,  _  Y vl v|tj 

v!v!t?  tp-t* 

First  order  coefficient  used  in  representing 
Yv|vj  as  a  function  of  (tj-1) 

Yvw 

Y 

Y  »  -  w 
™  "ipi2 

Coefficient  used  in  representing  Y  as  a 
function  of  the  product  vw 

Y 

wp 

Y 

y  •  =  -  yp 
wp  ipt3 

Coefiuient  used  in  representing  Y  as  a 
function  of  the  product  wp 

Ywr 

Y 

.  Y  '  =  wr 

**  $p*3 

Coefficient  used  in  representing  Y  as  a 
function  of  the  product  wr 

Y5r 

y  Y6r 

&r  $pi*U* 

First  order  coefficient  used  in  representing 

Y  a£  a  function  of  6r 

Y6r»j 

Y  .  _  Y  6rn 
«rJ?  ~£p42U2 

First  order  coefficient  used  in  representing 

Ygr  as  a  function  of  (tj-1) 
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* 


£ 

n|«o  ' 

II 

Normal  body  axis;  also  the  coordinate  of  a 
point  relative  to  the  origin  of  body  axes 

lB 

a  B. 

2b  -  r 

The  z  coordinate  of  CB 

ZG 

*c  -  1 

The  z  coordinate  of  CG 

zo 

,  ,  -  *o 
*o  "  T" 

A  coordinate  of  the  displacement  of  CG 
relative  to  the  origin  of  a  set  of  fixed  axes 

Z 

ft 

Z'  =  2 

ipt2u2 

Hydrodynamic  force  component  along  z 
axis  (normal  force) 

Z* 

z* 

Z*  *peja* 

Normal  force  when  body  angles  (a,  j3)  and 
control  surface  angles  are  sero 

z 

pp 

zpp  ~*p?* 

Second  order  coefficient  used  in  representing 
2  as  a  function  of  p.  First  order  coefficient 
is  zero 

z 

q 

2  '  = -^3 — 

3  ipt3U 

First  order  coefficient  used  in  representing 

Z  as  a  function  of  q 

Zqq 

2  '  =  .ZS3.. 

qn  $p-t3u 

First  order  coefficient  used  in  representing 
Zq  as  a  function  of  (17- 1) 

z4 

Z  i 

zi,=i^ 

Coefficient  used  in  representing  Z  as  a 
function  of  q 

Zlql6s 

zlq|6.'  = 

Coefficient  used  in  representing  Z  as  a 

function  of  q 

z 

*p 

Z  •  =  3*- 
rP  Jp4« 

Coefficient  used  in  representing  Z  as  a 
function  of  the  product  rp 

z 

XX 

♦  Z  V*. 

*y  i  —  ** 

**  $pt« 

Second  order  coefficient  used  in  representing 
Z  as  a  function  of  r.  First  order  coefficient 
is  zero 

Z 

w 

_  zw 

zw  "|**u 

First  order  coefficient  used  in  representing 

Z  as  a  function  of  w 

z 

WT? 

2  <  =  JEwn 
«*?  ipFu 

First  order  coefficient  used  in  representing 
Zw  as  a  function  of  (7)- 1) 

Z  . 

MV 

7  .  _  Zw 

w  -j^r 

Coefficient  used  in  representing  Z  is  a 
function  of  w 

Z,  1 

M 

7  i  _  Z!wl 

ZM  * iPi*u 

First  order  coefficient  used  in  representing 

Z  as  a  function  of  w;  equal  to  zero  for  sym¬ 
metrical  function 

Zwjq| 

2  ,  .•szfbdai 
wlqi  ioi3 

Coefficient  used  in  representing  Z  as  a 
function  of  q  w 
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The  goal  of  this  study  was  to  determine  the  optimum  sets  of  equations  to  be  used 
with  two  general  types  of  submarine  trainers.  The  equations  of  motion  used  were 
those  developed  by  the  Naval  Ship  Research  and  Development  Center.  They  were  re¬ 
duced  into  two  sets  for  use  in  training  simulators,  one  for  s  basic  submerged  con¬ 
trol  maneuvering  trainer  and  one  for  an  advanced  submerged  control  simulator  to 
provide  highly  realistic  ship  control  training  through  the  full  range  of  normal 
and  casualty  conditions. 

The  report  outlines  a  general-pur po3e  digital  computer  program,  following  the 
NSRDC  standard  equations,  written  in  FORTRAN.  Integration  methods  for  digital 
simulators  are  discussed.  A  number  of  programs  for  testing  the  degree  of  simu¬ 
lation  of  a  digital  simulator  program  are  given. 

Use  of  these  programs  as  applied  to  submarine  simulation  is  shown  with  two  sets 
of  equations  which  eliminate  73  of  the  139  coefficients  used  in  the  original 
NSRDC  equations.  The  mathematical  model  for  a  submarine  simulator  using  a  very 
small  computer  is  developed.  A  discussion  of  near-surface  operation  and  wave- 
generation  is  followed  by  the  general  requirements  for  determination  of  the 
simulation  requirements  for  training  a3  opposed  to  research. 
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